KJ135
/
ISCE_INSAR
mirror of http://172.16.0.12:5000/WBJY/ISCE_INSAR.git
4
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

gdal version smarts for lat/lon order

LT1AB
Russell Grew 2021-05-07 08:39:13 +10:00
parent e499aec055
commit 8ba8f32f67
1 changed files with 13 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
contrib/stack/topsStack/plotBursts.py
View File

@ -48,7 +48,7 @@ def main(iargs=None):
inps=cmdLineParse(iargs)
from osgeo import ogr, osr
from osgeo import gdal, ogr, osr
import matplotlib
if inps.shapefile is not None:
matplotlib.use('Agg')
@ -138,11 +138,18 @@ def main(iargs=None):
lateNear = burst.orbit.rdr2geo(t1,r0)
ring = ogr.Geometry(ogr.wkbLinearRing)
ring.AddPoint(earlyNear[0], earlyNear[1])
ring.AddPoint(earlyFar[0], earlyFar[1])
ring.AddPoint(lateFar[0], lateFar[1])
ring.AddPoint(lateNear[0], lateNear[1])
ring.AddPoint(earlyNear[0], earlyNear[1])
if int(gdal.__version__[0]) >= 3:
ring.AddPoint(earlyNear[0], earlyNear[1])
ring.AddPoint(earlyFar[0], earlyFar[1])
ring.AddPoint(lateFar[0], lateFar[1])
ring.AddPoint(lateNear[0], lateNear[1])
ring.AddPoint(earlyNear[0], earlyNear[1])
else:
ring.AddPoint(earlyNear[1], earlyNear[0])
ring.AddPoint(earlyFar[1], earlyFar[0])
ring.AddPoint(lateFar[1], lateFar[0])
ring.AddPoint(lateNear[1], lateNear[0])
ring.AddPoint(earlyNear[1], earlyNear[0])
feature = ogr.Feature(layer.GetLayerDefn())
feature.SetField('Name', 'IW{0}-{1}'.format(swath, ii))