def test_polar_image(self): ds = get_path('Mars_MGS_MOLA_ClrShade_MAP2_90.0N0.0_POLA.tif') testname = os.path.join(self.path, 'polartest.svg') s = bar.ScaleBar.from_image(ds, outputname = testname) ref, test = getsvg(get_path('mars_polar.svg'), testname) self.assertEqual(ref, test) os.remove(testname)
def test_lambert_image(self): ds = get_path('Lunar_LRO_LOLA_Shade_MAP2_90.0N20.0_LAMB.tif') testname = os.path.join(self.path, 'lambtest.svg') s = bar.ScaleBar.from_image(ds, outputname=testname) ref, test = getsvg(get_path('lunar_lamb.svg'), testname) self.assertEqual(ref, test) os.remove(testname)
def test_polar_image(self): ds = get_path('Mars_MGS_MOLA_ClrShade_MAP2_90.0N0.0_POLA.tif') testname = os.path.join(self.path, 'polartest.svg') s = bar.ScaleBar.from_image(ds, outputname=testname) ref, test = getsvg(get_path('mars_polar.svg'), testname) self.assertEqual(ref, test) os.remove(testname)
def getsvg(im1, im2): reference_merc = get_path(im1) with open(reference_merc, 'r') as f: ref = f.readlines() with open(im2, 'r') as f: test = f.readlines() return ref, test
def test_mercator_wktstring(self): wkt = """PROJCS["Mercator_MARS", GEOGCS["GCS_MARS", DATUM["MARS", SPHEROID["MARS",3396190,169.8944472236118]], PRIMEM["Reference_Meridian",0], UNIT["Degree",0.017453292519943295]], PROJECTION["Mercator_1SP"], PARAMETER["central_meridian",0], PARAMETER["false_easting",0], PARAMETER["false_northing",0], UNIT["Meter",1], PARAMETER["latitude_of_origin",0.0]]""" testname = os.path.join(self.path,'mars_merc_wkt.svg') s = bar.ScaleBar.from_projstring(wkt, (0, 0, 180,65), outputname=testname) ref, test = getsvg(get_path('mars_merc_wkt.svg'), testname) self.assertEqual(ref, test) os.remove(testname)
def test_polar_wktstring(self): wkt = """PROJCS["Mars_South_Pole_Stereographic", GEOGCS["Mars 2000", DATUM["D_Mars_2000", SPHEROID["Mars_2000_IAU_IAG",3396190.0,169.89444722361179]], PRIMEM["Greenwich",0], UNIT["Decimal_Degree",0.0174532925199433]], PROJECTION["Stereographic"], PARAMETER["False_Easting",0], PARAMETER["False_Northing",0], PARAMETER["Central_Meridian",0], PARAMETER["Scale_Factor",1], PARAMETER["Latitude_Of_Origin",-90], UNIT["Meter",1]]""" testname = os.path.join(self.path,'mars_polar_wkt.svg') s = bar.ScaleBar.from_projstring(wkt, (0, -90, 180,-40), outputname=testname, cliplat=-90) ref, test = getsvg(get_path('mars_polar_wkt.svg'), testname) self.assertEqual(ref, test) os.remove(testname)
def test_mercator_wktstring(self): wkt = """PROJCS["Mercator_MARS", GEOGCS["GCS_MARS", DATUM["MARS", SPHEROID["MARS",3396190,169.8944472236118]], PRIMEM["Reference_Meridian",0], UNIT["Degree",0.017453292519943295]], PROJECTION["Mercator_1SP"], PARAMETER["central_meridian",0], PARAMETER["false_easting",0], PARAMETER["false_northing",0], UNIT["Meter",1], PARAMETER["latitude_of_origin",0.0]]""" testname = os.path.join(self.path, 'mars_merc_wkt.svg') s = bar.ScaleBar.from_projstring(wkt, (0, 0, 180, 65), outputname=testname) ref, test = getsvg(get_path('mars_merc_wkt.svg'), testname) self.assertEqual(ref, test) os.remove(testname)
def test_lambert_wktstring(self): wkt = """PROJCS["Moon_Lambert_Conformal_Conic_AUTO", GEOGCS["Moon 2000", DATUM["D_Moon_2000", SPHEROID["Moon_2000_IAU_IAG",1737400.0,0.0]], PRIMEM["Greenwich",0], UNIT["Decimal_Degree",0.0174532925199433]], PROJECTION["Lambert_Conformal_Conic"], PARAMETER["False_Easting",0], PARAMETER["False_Northing",0], PARAMETER["Central_Meridian",0], PARAMETER["Standard_Parallel_1",43], PARAMETER["Standard_Parallel_2",73], PARAMETER["Latitude_Of_Origin",0], UNIT["Meter",1]]""" testname = os.path.join(self.path, 'lunar_lamb_wkt.svg') s = bar.ScaleBar.from_projstring(wkt, (33.25, 36, 60, 78.11), outputname=testname) ref, test = getsvg(get_path('lunar_lamb_wkt.svg'), testname) self.assertEqual(ref, test) os.remove(testname)
def test_polar_wktstring(self): wkt = """PROJCS["Mars_South_Pole_Stereographic", GEOGCS["Mars 2000", DATUM["D_Mars_2000", SPHEROID["Mars_2000_IAU_IAG",3396190.0,169.89444722361179]], PRIMEM["Greenwich",0], UNIT["Decimal_Degree",0.0174532925199433]], PROJECTION["Stereographic"], PARAMETER["False_Easting",0], PARAMETER["False_Northing",0], PARAMETER["Central_Meridian",0], PARAMETER["Scale_Factor",1], PARAMETER["Latitude_Of_Origin",-90], UNIT["Meter",1]]""" testname = os.path.join(self.path, 'mars_polar_wkt.svg') s = bar.ScaleBar.from_projstring(wkt, (0, -90, 180, -40), outputname=testname, cliplat=-90) ref, test = getsvg(get_path('mars_polar_wkt.svg'), testname) self.assertEqual(ref, test) os.remove(testname)
# -*- coding: utf-8 -*- import sys import numpy as np import matplotlib.pylab as plt from scipy.interpolate import interp1d from scalebar.examples import get_path from scalebar.fileio import gdalio from scalebar.utils import util from scalebar.bar import bar if __name__ == '__main__': #Read file and get projection information ds = get_path('Mars_MGS_MOLA_ClrShade_MAP2_0.0N0.0_MERC.tif') #ds = get_path('Mars_MGS_MOLA_ClrShade_MAP2_90.0N0.0_POLA.tif') #ds = get_path('Lunar_LRO_LOLA_Shade_MAP2_90.0N20.0_LAMB.tif') S = bar.ScaleBar.from_image(ds, outputname='mars_merc.svg', lon_major_ticks=[25,50,75])
def setUp(self): self.ds = gdalio.GeoDataSet(get_path('Mars_MGS_MOLA_ClrShade_MAP2_90.0N0.0_POLA.tif'))
def setUp(self): self.ds = gdalio.GeoDataSet(get_path('Lunar_LRO_LOLA_Shade_MAP2_90.0N20.0_LAMB.tif'))