I have a data set in one projection, and a shapefile in another, and I am trying to combine the two. I am using Python. I apologize for posting a "help me transform coordinates" question on this site - there are so many.
I assume the way to do this conversion is with proj4. I am unable to find/understand the proj4 string needed to do this conversion.
The shapefile.prj contains this:
PROJCS["North_Pole_Stereographic",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Stereographic"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",0.0],PARAMETER["Scale_Factor",1.0],PARAMETER["Latitude_Of_Origin",90.0],UNIT["Meter",1.0]]Is there a simple way to find the proj4 code for this?
And I have a netCDF file that I can plot correctly with the following Python code:
map = Basemap(llcrnrlon=lon[0,0], llcrnrlat=lat[0,0], urcrnrlon=lon[560,300], urcrnrlat=lat[560,300], #rsphere=(6378137.00,6356752.3142),\ resolution='l', area_thresh=1E6, projection='lcc', #lat_1=90, lon_0=0) lat_1=71.0, lon_0=-39.0)I would like to keep the Basemap (python) projection, and convert the shapefile contents to that projection.
Edit
Further research and I found this post: Shapefile PRJ to PostGIS SRID lookup table? with a helpful bit of code. Running that and I find the following:
Proj4 is: +proj=stere +lat_0=90 +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defsEPSG is: NoneEdit
I found out I can determine the basemap proj string with:
map.proj4stringI'm now using the advice here http://all-geo.org/volcan01010/2012/11/change-coordinates-with-pyproj/ to try to convert from one to the other...
أكثر...
I assume the way to do this conversion is with proj4. I am unable to find/understand the proj4 string needed to do this conversion.
The shapefile.prj contains this:
PROJCS["North_Pole_Stereographic",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Stereographic"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",0.0],PARAMETER["Scale_Factor",1.0],PARAMETER["Latitude_Of_Origin",90.0],UNIT["Meter",1.0]]Is there a simple way to find the proj4 code for this?
And I have a netCDF file that I can plot correctly with the following Python code:
map = Basemap(llcrnrlon=lon[0,0], llcrnrlat=lat[0,0], urcrnrlon=lon[560,300], urcrnrlat=lat[560,300], #rsphere=(6378137.00,6356752.3142),\ resolution='l', area_thresh=1E6, projection='lcc', #lat_1=90, lon_0=0) lat_1=71.0, lon_0=-39.0)I would like to keep the Basemap (python) projection, and convert the shapefile contents to that projection.
Edit
Further research and I found this post: Shapefile PRJ to PostGIS SRID lookup table? with a helpful bit of code. Running that and I find the following:
Proj4 is: +proj=stere +lat_0=90 +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defsEPSG is: NoneEdit
I found out I can determine the basemap proj string with:
map.proj4stringI'm now using the advice here http://all-geo.org/volcan01010/2012/11/change-coordinates-with-pyproj/ to try to convert from one to the other...
أكثر...