GPS文件处理
GPS数据导入计算机后，有多种格式，易于处理、易于识别的就是TXT和CSV格式，在本文中我处理的是TXT格式。
建立Shapefile和处理TXT数据都比较简单，后面我就只贴代码（代码里有详细注释），不讲原理了。
主代码：
public static void CreateShpfromTXT(IMapControl3 map_control,progressFrm progressfrm)
{
ProgressBar progressBar1 = progressfrm.progressBar1;
OpenFileDialog openDia = new OpenFileDialog();
openDia.Title = "打开GPS文件(TXT).......";
openDia.Filter = "文本文件(*.TXT)|*.txt";
openDia.RestoreDirectory = true;
openDia.Multiselect = true;
if (openDia.ShowDialog() == DialogResult.OK)
{
string[] arrayFileName = openDia.FileNames;
string[] arraySafeName = openDia.SafeFileNames;
for (int i = arrayFileName.Length-1; i >=0; i--)
{
//根据用户选择的多个TXT文件 依次生成Shp文件
string filefullname = arrayFileName[i];
string filesafename = arraySafeName[i];
string directoryPath = filefullname.Remove(filefullname.LastIndexOf(filesafename), filesafename.Length);
string fileName = filesafename.Remove(filesafename.LastIndexOf(".TXT"), 4);
string shapeName = fileName + ".shp";
foreach (FileInfo file in new DirectoryInfo(directoryPath).GetFiles())
{
if (file.Name == shapeName)
{
file.Delete();//如果该Shp文件已经存在，则删除它以及与它相关的文件
(new FileInfo(directoryPath + fileName + ".shx")).Delete();
(new FileInfo(directoryPath + fileName + ".prj")).Delete();
(new FileInfo(directoryPath + fileName + ".dbf")).Delete();
}
}
#region 开始创建Shp
IFeatureClass newfeatureclass = CreatNewShapeFile(directoryPath, shapeName);
if (newfeatureclass != null)
{
string[] strs = File.ReadAllLines(filefullname);
int hasline = strs.Count();//获取该TXT文件有多少个点
progressBar1.Maximum = hasline;
progressBar1.Value = 0;
progressfrm.TipText ="";
progressfrm.TipText = string.Format("正在处理GPS文件：{0},请稍等……", filefullname);
progressfrm.Refresh();
System.IO.FileStream filestreamTXT = new System.IO.FileStream(filefullname, FileMode.Open);
StreamReader streamreader = new StreamReader(filestreamTXT);
streamreader.BaseStream.Seek(0, SeekOrigin.Begin);
streamreader.ReadLine(); streamreader.ReadLine(); streamreader.ReadLine();
while (streamreader.Peek() > 0)
{
string strLine = streamreader.ReadLine() + " ";
List<string> filedArray = findWord(strLine);
if (filedArray != null)
{
List<string> filedArrayUnion = GetUnionfieldList(filedArray);
IFeature newfeature = newfeatureclass.CreateFeature();
IPoint newpoint = new PointClass();
newpoint.X = Convert.ToDouble(filedArrayUnion[2]);
newpoint.Y = Convert.ToDouble(filedArrayUnion[1]);
newpoint.Z = Convert.ToDouble(filedArrayUnion[3]);
newfeature.Shape = newpoint as IGeometry;
newfeature.Store();
progressBar1.Value++;
}
}
strs = null;
}
#endregion
map_control.AddShapeFile(directoryPath, shapeName);
}
}
}
创建一个空的Shapefile：
public static IFeatureClass CreatNewShapeFile(string directoryPath, string shapeName)
{
IAoInitialize aoinitialize = new AoInitializeClass();
aoinitialize.Initialize(esriLicenseProductCode.esriLicenseProductCodeEngine);
//步骤：打开工作空间；建立字段集；建立字段；IFeatureWorkspace.CreatFeatureClass
//关键在于几何字段的IGeometryDefEdit，另外还有esriFieldType.esriFieldTypeGeometry，esriFeatureType.esriFTSimple
const string strShapeFieldName = "shape";
IWorkspaceFactory pWSF = new ShapefileWorkspaceFactoryClass();
IFeatureWorkspace pWS = (IFeatureWorkspace)pWSF.OpenFromFile(directoryPath, 0);
// directoryPath 父文件夹
//设置字段集
IFields pFields = new FieldsClass();
IFieldsEdit pFieldsEdit = (IFieldsEdit)pFields;
//设置字段
IField pFieldgeo = new FieldClass();
IFieldEdit pFieldEditgeo = (IFieldEdit)pFieldgeo;
//创建类型为几何类型的字段，添加Shape字段，OID是不用添加的
pFieldEditgeo.Name_2 = "Shape";//Name代表只读，Name_2代表只写
pFieldEditgeo.Type_2 = esriFieldType.esriFieldTypeGeometry;
//为esriFieldTypeGeometry类型的字段创建几何定义，包括类型和空间参照
IGeometryDef pGeoDef = new GeometryDefClass(); //The geometry definition for the field if IsGeometry is TRUE.
IGeometryDefEdit pGeoDefEdit = (IGeometryDefEdit)pGeoDef;
pGeoDefEdit.GeometryType_2 = esriGeometryType.esriGeometryPoint;
// pGeoDefEdit.SpatialReference_2 = new UnknownCoordinateSystemClass() ;
// pFieldEditgeo.GeometryDef_2 = pGeoDef;
ISpatialReferenceFactory3 spatialReferenceFactory = new SpatialReferenceEnvironmentClass();
ISpatialReference spatialReference = spatialReferenceFactory.CreateSpatialReference((int)esriSRGeoCSType.esriSRGeoCS_WGS1984);
ISpatialReferenceInfo spatialReferenceInfo = spatialReference as ISpatialReferenceInfo;
ISpatialReference spatial = spatialReferenceInfo as ISpatialReference;
pGeoDefEdit.SpatialReference_2 = spatial;
pFieldEditgeo.GeometryDef_2 = pGeoDef;
pFieldsEdit.AddField(pFieldgeo);
//创建shapefile
IFeatureClass newfeatureclass = pWS.CreateFeatureClass(shapeName, pFields, null, null, esriFeatureType.esriFTSimple, "Shape", "");
return newfeatureclass;
}
