C#用GDAL创建SHP文件,升成的PRJ文件在ARCGIS中不识别的问题

在用GDAL创建SHP文件时,用SpatialReference类创建的坐标系统,在升成时会有一个PRJ文件,但是该文件在用ARCGIS不会被识别,原因为ARCGIS版本不一样所识别的文件也不一样)采用另外一种方法,直接升成PRJ文件,该文件必须要按照ASCII码升成,否则ARCGIS不识别
思路,用坐标取带号,该方法只对带上带号的2000高期吕克坐标系有用。

int dh = Convert.ToInt16(zb[0].X.ToString().Substring(0, 2));
string wkt = 
                                "PROJCS[\"CGCS2000 / 3-degree Gauss-Kruger zone"+dh.ToString() +"\"," +
                                "GEOGCS[\"China Geodetic Coordinate System 2000\"," +
                                "DATUM[\"China_2000\"," +
                                "SPHEROID[\"CGCS2000\",6378137.0,298.257222101]]," + 
                                "PRIMEM[\"Greenwich\",0.0]," +
                                "UNIT[\"Degree\",0.0174532925199433]]," +
                                "PROJECTION[\"Gauss_Kruger\"]," +
                                "PARAMETER[\"False_Easting\","+dh.ToString()+"500000.0]," +
                                "PARAMETER[\"False_Northing\",0.0]," +
                                "PARAMETER[\"Central_Meridian\","+(dh*3).ToString("F1")+"]," +
                                "PARAMETER[\"Scale_Factor\",1.0]," +
                                "PARAMETER[\"Latitude_of_Origin\",0.0]," +
                                "UNIT[\"Metre\",1.0]]";
 string prjPath = zdPath + @"\ZD.prj";
 using (StreamWriter write = new StreamWriter(prjPath, false, Encoding.ASCII))
     {
         write.WriteLine(wkt);
     }

最重要是升成PRJ该文件必须是ASCII码升成,否则ARCGIS不会识别,以上代码在ARCGIS10.6中实测通过。

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
要解决GDAL读取Shapefile文件后输出文乱码的问题,可以通过设置合适的编码来处理。下面是一个示例代码,演示如何使用GDAL读取Shapefile文件正确处理文字符: ```python from osgeo import ogr # 打开Shapefile dataset = ogr.Open('input_shapefile.shp', 0) layer = dataset.GetLayer() # 获取属性表编码 layer_encoding = layer.GetLayerDefn().GetFieldDefn(0).GetEncoding() # 设置输出编码 output_encoding = 'utf-8' # 遍历要素并输出属性值 for feature in layer: for i in range(feature.GetFieldCount()): field_value = feature.GetField(i) # 解码属性值为Unicode字符串 field_value_unicode = field_value.decode(layer_encoding) # 编码为指定输出编码的字符串 field_value_encoded = field_value_unicode.encode(output_encoding) # 输出属性值 print(field_value_encoded) # 关闭数据集 dataset = None ``` 在上述代码,我们首先通过`GetLayerDefn()`方法获取属性表的编码,并将其保存在`layer_encoding`变量。然后,我们将输出编码设置为`utf-8`,可以根据需要修改为其他编码。接下来,我们遍历要素并使用`decode()`方法将属性值解码为Unicode字符串,然后使用`encode()`方法将其编码为指定的输出编码字符串。最后,我们使用`print()`语句输出属性值。 请将`input_shapefile.shp`替换为您要读取的Shapefile文件路径。如果您的Shapefile文件使用了其他编码,请在`layer.GetLayerDefn().GetFieldDefn(0).GetEncoding()`指定正确的编码。 希望这可以解决您的问题!如果您有任何进一步的疑问,请随时提问。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值