Ben senin için ölçekleme yapacak arcpy API hiçbir şey farkında değilim, ama bunu yapmak için bir işlev yazmak nispeten basit olacaktır.
Aşağıdaki kod 2D özelliklerin ölçeklendirmesini yapar ve M veya Z değerlerini dikkate almaz:
import arcpy
import math
def scale_geom(geom, scale, reference=None):
"""Returns geom scaled to scale %"""
if geom is None: return None
if reference is None:
# we'll use the centroid if no reference point is given
reference = geom.centroid
refgeom = arcpy.PointGeometry(reference)
newparts = []
for pind in range(geom.partCount):
part = geom.getPart(pind)
newpart = []
for ptind in range(part.count):
apnt = part.getObject(ptind)
if apnt is None:
# polygon boundaries and holes are all returned in the same part.
# A null point separates each ring, so just pass it on to
# preserve the holes.
newpart.append(apnt)
continue
bdist = refgeom.distanceTo(apnt)
bpnt = arcpy.Point(reference.X + bdist, reference.Y)
adist = refgeom.distanceTo(bpnt)
cdist = arcpy.PointGeometry(apnt).distanceTo(bpnt)
# Law of Cosines, angle of C given lengths of a, b and c
angle = math.acos((adist**2 + bdist**2 - cdist**2) / (2 * adist * bdist))
scaledist = bdist * scale
# If the point is below the reference point then our angle
# is actually negative
if apnt.Y < reference.Y: angle = angle * -1
# Create a new point that is scaledist from the origin
# along the x axis. Rotate that point the same amount
# as the original then translate it to the reference point
scalex = scaledist * math.cos(angle) + reference.X
scaley = scaledist * math.sin(angle) + reference.Y
newpart.append(arcpy.Point(scalex, scaley))
newparts.append(newpart)
return arcpy.Geometry(geom.type, arcpy.Array(newparts), geom.spatialReference)
Bir geometri nesnesi, bir ölçek faktörü (1 = aynı boyut, 0.5 = yarım boyut, 5 = 5 kat daha büyük, vb.) Ve isteğe bağlı bir referans noktası ile arayabilirsiniz:
scale_geom(some_geom, 1.5)
Hedef özellik sınıfının zaten var olduğunu varsayarak, tüm özellik sınıfını ölçeklendirmek için imleçlerle birlikte bunu kullanın:
incur = arcpy.da.SearchCursor('some_folder/a_fgdb.gdb/orig_fc', ['OID@','SHAPE@'])
outcur = arcpy.da.InsertCursor('some_folder/a_fgdb.gdb/dest_fc', ['SHAPE@'])
for row in incur:
# Scale each feature by 0.5 and insert into dest_fc
outcur.insertRow([scale_geom(row[1], 0.5)])
del incur
del outcur
edit: İşte 0.5 ve 5 kez test geometrinizin yaklaşık bir tahminini kullanan bir örnek:
Ayrıca çok halkalı çokgenlerle (delikler) test edilmiştir!
İstendiği gibi bir açıklama:
scale_geom
tek bir çokgen alır ve ondan bir referans noktasına olan mesafeyi ölçerek (varsayılan olarak, çokgenin centroidi) her bir tepe noktasında dolaşır.
Bu mesafe daha sonra yeni 'ölçeklendirilmiş' tepe noktasını oluşturmak için verilen ölçekle ölçeklendirilir.
Ölçekleme esasen referans noktasından orijinal tepe noktasına kadar ölçeklendirilmiş uzunlukta bir çizgi çizilerek yapılır, çizginin sonu ölçeklenmiş tepe noktası haline gelir.
Açı ve döndürme öğeleri oradadır, çünkü çizginin sonunun konumunu tek bir eksen boyunca hesaplamak ve daha sonra 'yerine' döndürmek daha düzdür.