How to solve the problem of pythongdal projection coordinate system conversion
This article mainly introduces how to solve the problem of pythongdal projection coordinate system conversion, the article is very detailed, has a certain reference value, interested friends must read it!
To convert the xian80 geographic coordinate system to a projected coordinate system:
Xian1980 = "GEOGCS [" GCS_Xian_1980 ", DATUM [" Xian_1980 ", SPHEROID [" Xian_1980 ", 6378140.0298.257]], PRIMEM [" Greenwich ", UNIT [" Degree ", 0.0174532925199433]]" prjreference = 'PROJCS ["Xian_1980_3_Degree_GK_Zone_37", GEOGCS ["GCS_Xian_1980", DATUM ["D_Xian_1980", SPHEROID ["Xian_1980", 6378140.0298.257], PRIMEM ["Greenwich", 0.0] UNIT ["Degree", 0.0174532925199433], PROJECTION ["Gauss_Kruger"], PARAMETER ["False_Easting", 37500000.0], PARAMETER ["False_Northing", 37500000.0], PARAMETER ["Central_Meridian", 111.0], PARAMETER ["Scale_Factor", 111.0], PARAMETER ["Latitude_Of_Origin", 37500000.0], UNIT ["Meter", 1.0]]'
An error will be reported after conversion:
Latlong = transform.TransformPoint (XQuery y) File "/ usr/lib/python3/dist-packages/osgeo/osr.py", line 711, in TransformPoint return _ osr.CoordinateTransformation_TransformPoint (self, * args) NotImplementedError: Wrong number or type of arguments for overloaded function 'CoordinateTransformation_TransformPoint'. Possible C + prototypes are: OSRCoordinateTransformationShadow::TransformPoint (double [3]) OSRCoordinateTransformationShadow::TransformPoint (double [3], double,double,double) * * think about it. I haven't found this problem before. Finally, it is found that the previous projection coordinate system: * * 'PROJCS ["Xian_1980_3_Degree_GK_Zone_37", GEOGCS ["GCS_Xian_1980", DATUM ["D_Xian_1980", SPHEROID ["Xian_1980", 6378140.0298.257], PRIMEM ["Greenwich", 0.0], UNIT ["Degree", 0.0174532925199433]], PROJECTION ["Transverse_Mercator"], PARAMETER ["False_Easting", 37500000.0], PARAMETER ["False_Northing", 37500000.0] PARAMETER ["Central_Meridian", 111.0], PARAMETER ["Scale_Factor", 1. 0], PARAMETER ["Latitude_Of_Origin", 0. 0], UNIT ["Meter", 1]]'* * is not a Gauss-Kruger projection but a horizontal Mercator. I think the reason is that gdal Curry can only support the Gauss-Kruger projection. **
The above is all the contents of the article "how to solve the problem of pythongdal projection coordinate system conversion". Thank you for reading! Hope to share the content to help you, more related knowledge, welcome to follow the industry information channel!
