/// <summary>
/// Projection name is written as suffix to the filename in case of reprojection.
/// Let's try to find file with the same name, but correct suffix for projection.
/// It's assumed here that projection for the layer passed doesn't match the project one.
/// </summary>
public static bool SeekSubstituteFile(LayerSource layer, MapWinGIS.GeoProjection targetProjection, out LayerSource newLayer)
{
newLayer = null;
string testName = CoordinateTransformation.FilenameWithProjectionSuffix(layer.Filename, layer.Projection, targetProjection);
if (File.Exists(testName))
{
LayerSource layerTest = new LayerSource(testName);
if (layerTest.Type != layer.Type)
{
return(false);
}
FileInfo f1 = new FileInfo(layer.Filename);
FileInfo f2 = new FileInfo(testName);
// the size of .shp files must be exactly the same
bool equalSize = !(layerTest.Type == LayerSourceType.Shapefile && f1.Length != f2.Length);
if (layerTest.Projection.get_IsSameExt(targetProjection, layerTest.Extents, 10) && equalSize)
{
newLayer = layerTest;
return(true);
}
}
return(false);
}
/// <summary>
/// 测试一个类型是由拓展形成的图层的单层投影
/// </summary>
/// <param name="filename">文件名</param>
/// <param name="newName">输出新的文件名</param>
/// <returns>返回测试结果</returns>
public TestingResult TestLayer(string filename, out string newName)
{
newName = filename;
LayerSource layer = new LayerSource(filename, this as MapWinGIS.ICallback); //打开指定文件的图层
if (layer.Type == LayerSourceType.Undefined) //图层类型是没有定义的
{
string message = layer.GetErrorMessage();
if (message == "")
{
message = "Unspecified error";
}
MessageBox.Show("Invalid datasource: " + message.ToLower() + Environment.NewLine + filename, m_mapWin.ApplicationInfo.ApplicationName, MessageBoxButtons.OK, MessageBoxIcon.Error);
return(TestingResult.Error);
}
else
{
LayerSource newLayer = null;
TestingResult result = this.TestLayer(layer, out newLayer);
if (result == TestingResult.Substituted)
{
newName = newLayer.Filename;
newLayer.Close();
}
layer.Close();
return(result);
}
}
/// <summary>
/// Runs the operation
/// </summary>
private void btnOk_Click(object sender, EventArgs e)
{
CoordinateSystem cs = this.ProjectionTreeView1.SelectedCoordinateSystem;
if (cs == null)
{
MessageBox.Show("No projection is selected", m_mapWin.ApplicationInfo.ApplicationName, MessageBoxButtons.OK, MessageBoxIcon.Information);
}
else
{
if (this.LayersControl1.Filenames.Count() == 0)
{
MessageBox.Show("No files are selected", m_mapWin.ApplicationInfo.ApplicationName,
MessageBoxButtons.OK, MessageBoxIcon.Information);
}
else
{
MapWinGIS.GeoProjection projection = new MapWinGIS.GeoProjection();
if (projection.ImportFromEPSG(cs.Code))
{
frmTesterReport report = new frmTesterReport();
int count = 0; // number of successfully processed files
foreach (string name in this.LayersControl1.Filenames)
{
LayerSource layer = new LayerSource(name);
string projName = layer.Projection != null ? layer.Projection.Name : "";
if (layer.Type != LayerSourceType.Undefined)
{
layer.Projection = projection;
count++;
}
else
{
report.AddFile(name, projName, ProjectionOperaion.Skipped, "");
}
}
if (count > 0)
{
MessageBox.Show(string.Format("The projection was successfully assigned to the files: {0}", count), m_mapWin.ApplicationInfo.ApplicationName,
MessageBoxButtons.OK, MessageBoxIcon.Information);
}
if (report.MismatchedCount > 0)
{
report.ShowReport(projection, "The following files were not processed:", ReportType.Assignment);
}
this.LayersControl1.UpdateProjections();
}
else
{
MessageBox.Show("Failed to initialize the selected projection", m_mapWin.ApplicationInfo.ApplicationName,
MessageBoxButtons.OK, MessageBoxIcon.Information);
}
}
}
}
/// <summary>
/// Reprojects image
/// </summary>
public static TestingResult Reproject(MapWinGIS.Image image, out MapWinGIS.Image imageNew, MapWinGIS.GeoProjection projection, frmTesterReport report)
{
MapWinGIS.GeoProjection projImage = new MapWinGIS.GeoProjection();
projImage.ImportFromProj4(image.GetProjection());
string sourcePrj = image.GetProjection();
string targetPrj = projection.ExportToProj4();
string origFilename = image.Filename;
MapWinGIS.ICallback callback = image.GlobalCallback;
imageNew = null;
LayerSource obj = new LayerSource(image);
LayerSource objNew = null;
if (CoordinateTransformation.SeekSubstituteFile(obj, projection, out objNew))
{
imageNew = objNew.Image;
return(TestingResult.Substituted);
}
string newFilename = CoordinateTransformation.FilenameWithProjectionSuffix(origFilename, projImage, projection);
newFilename = CoordinateTransformation.GetSafeNewName(newFilename);
// setting callback
if (report != null)
{
if (!report.Visible)
{
report.InitProgress(projection);
}
report.ShowFilename(image.Filename);
}
MapWinGIS.GeoProcess.SpatialReference.ProjectImage(sourcePrj, targetPrj, origFilename, newFilename, report);
if (report != null)
{
report.ClearFilename();
}
imageNew = new MapWinGIS.Image();
if (imageNew.Open(newFilename, MapWinGIS.ImageType.USE_FILE_EXTENSION, false, callback))
{
return(TestingResult.Ok);
}
else
{
imageNew = null;
return(TestingResult.Error);
}
}
/// <summary>
/// Reprojects a grid
/// </summary>
/// <param name="grid">Grid object to reproject</param>
/// <param name="inPlace">Whether reprojected file should replace the initial one</param>
/// <returns>True on success and false otherwise</returns>
public static TestingResult Reproject(MapWinGIS.Grid grid, out MapWinGIS.Grid gridNew, MapWinGIS.GeoProjection projection, frmTesterReport report)
{
string sourcePrj = grid.Header.GeoProjection.ExportToProj4();
string targetPrj = projection.ExportToProj4();
string origFilename = grid.Filename;
MapWinGIS.ICallback callback = grid.GlobalCallback;
gridNew = null;
LayerSource obj = new LayerSource(grid);
LayerSource objNew = null;
if (CoordinateTransformation.SeekSubstituteFile(obj, projection, out objNew))
{
gridNew = objNew.Grid;
return(TestingResult.Substituted);
}
string newFilename = CoordinateTransformation.FilenameWithProjectionSuffix(origFilename, grid.Header.GeoProjection, projection);
newFilename = CoordinateTransformation.GetSafeNewName(newFilename);
// setting callback
if (report != null)
{
if (!report.Visible)
{
report.InitProgress(projection);
}
report.ShowFilename(grid.Filename);
}
bool result = MapWinGIS.GeoProcess.SpatialReference.ProjectGrid(ref sourcePrj, ref targetPrj, ref origFilename, ref newFilename, true, report);
if (report != null)
{
report.ClearFilename();
}
if (!result)
{
return(TestingResult.Error);
}
// TODO: no need to open it if only a name is supposed to be returned
gridNew = new MapWinGIS.Grid();
gridNew.Open(newFilename, MapWinGIS.GridDataType.UnknownDataType, false, MapWinGIS.GridFileType.UseExtension, callback);
gridNew.AssignNewProjection(projection.ExportToProj4());
return(TestingResult.Ok);
}
/// <summary>
/// Resolve a MaybeTask<IImageProvider> that represents the whole layer stack.
/// </summary>
/// <param name="backgroundImage">A MaybeTask<IImageProvider> resulting in the image to use for the background layer.</param>
/// <param name="backgroundImageSize">Optionally, the size of the background image. If omitted, the natural size of the source will be used, but the size will not available in the LayerLinkingContext for lazy resolvers.</param>
/// <param name="renderSizeHint">Optionally, the size of image that the layer stack will be rendered to. Note that this is just a hint, the resulting <see cref="IImageProvider" /> should also be rendered at this size by the user.</param>
/// <returns>An IImageProvider that represents the whole layer stack.</returns>
public MaybeTask <IImageProvider> ToImageProvider(MaybeTask <IImageProvider> backgroundImage, Size backgroundImageSize = default(Size), Size renderSizeHint = default(Size))
{
var backgroundLayer = GetActiveBackgroundLayer(backgroundImage, backgroundImageSize);
var layerContextInvariants = new LayerContext.Invariants(backgroundLayer, backgroundImage, renderSizeHint);
var layerContexts = new LayerContext[Count];
layerContexts[0] = new LayerContext(layerContextInvariants, null, backgroundLayer, 0);
backgroundImage = backgroundLayer.GetImageProvider(layerContexts[0]);
var currentImage = backgroundImage;
for (int layerIndex = 1; layerIndex < layerContexts.Length; ++layerIndex)
{
Debug.WriteLine("Resolving layer " + layerIndex);
var context = new LayerContext(layerContextInvariants, layerContexts[layerIndex - 1].CurrentLayer, this[layerIndex], layerIndex);
context.PreviousImage = currentImage;
currentImage = context.CurrentLayer.GetImageProvider(context);
Debug.Assert(!currentImage.IsEmpty, "Resolved imageprovider is null.");
// Bind the previous image as the source of the current one.
currentImage = LayerSource.Bind(currentImage, context.PreviousImage);
// If blending is used, wrap the current image in a BlendEffect and bind the source to it as well.
var blendedCurrentImage = LayerBlend.Apply(context, currentImage);
if (!blendedCurrentImage.IsEmpty)
{
blendedCurrentImage = LayerSource.Bind(blendedCurrentImage, context.PreviousImage);
currentImage = blendedCurrentImage;
}
layerContexts[layerIndex] = context;
}
Debug.WriteLine("Resolved final image provider.");
return(currentImage);
}
/// <summary>
/// Reprojects layer of undefined type
/// </summary>
/// <param name="filename">Filename of layer to reproject. A changed name will be returned when new file was created</param>
/// <param name="projection">New projection</param>
/// <param name="inPlace">Whether new files should be written</param>
/// <param name="report">A reference to report form</param>
/// <returns>True on success and false otherwise</returns>
public static TestingResult ReprojectLayer(string filename, out string newFilename, MapWinGIS.GeoProjection projection, frmTesterReport callback)
{
newFilename = "";
LayerSource source = new LayerSource(filename, callback);
if (source.Type != LayerSourceType.Undefined)
{
LayerSource sourceNew = null;
TestingResult result = CoordinateTransformation.ReprojectLayer(source, out sourceNew, projection, callback);
if (sourceNew != null)
{
newFilename = sourceNew.Filename;
}
return(result);
}
else
{
return(TestingResult.Error);
}
}
/// <summary>
/// Reprojects layer of undefined type
/// </summary>
/// <param name="filename">Filename of layer to reproject. A changed name will be returned when new file was created</param>
/// <param name="projection">New projection</param>
/// <param name="inPlace">Whether new files should be written</param>
/// <param name="report">A reference to report form</param>
/// <returns>True on success and false otherwise</returns>
public static TestingResult ReprojectLayer(LayerSource layer, out LayerSource newLayer, MapWinGIS.GeoProjection projection, frmTesterReport report)
{
newLayer = null;
TestingResult result = TestingResult.Error;
switch (layer.Type)
{
case LayerSourceType.Shapefile:
MapWinGIS.Shapefile sfNew = null;
result = CoordinateTransformation.Reproject(layer.Shapefile, out sfNew, projection, report);
if (sfNew != null)
{
newLayer = new LayerSource(sfNew);
}
break;
case LayerSourceType.Grid:
MapWinGIS.Grid gridNew = null;
result = CoordinateTransformation.Reproject(layer.Grid, out gridNew, projection, report);
if (gridNew != null)
{
newLayer = new LayerSource(gridNew);
}
break;
case LayerSourceType.Image:
MapWinGIS.Image imageNew = null;
result = CoordinateTransformation.Reproject(layer.Image, out imageNew, projection, report);
if (imageNew != null)
{
newLayer = new LayerSource(imageNew);
}
break;
default:
System.Diagnostics.Debug.Print("Coordinate transformation: unsupported interface");
break;
}
return(result);
}
/// <summary>
/// 测试单层投影
/// </summary>
/// <param name="layer">图层源（Shapefile或grid对象）</param>
/// <param name="newLayer">输出新图层</param>
/// <returns>返回测试结果</returns>
public TestingResult TestLayer(LayerSource layer, out LayerSource newLayer)
{
if (layer == null)
{
throw new ArgumentException("空图层引用被通过");
}
newLayer = null;
MapWinGIS.GeoProjection projectProj = m_mapWin.Project.GeoProjection;
MapWinGIS.GeoProjection layerProj = layer.Projection;
bool isSame = projectProj.get_IsSameExt(layerProj, layer.Extents, 10);
// 让我们看看我们是否有项目的投影的一种方言
if (!isSame && !projectProj.IsEmpty)
{
ProjectionDatabase db = (ProjectionDatabase)m_mapWin.ProjectionDatabase;//投影数据源
if (db != null)
{
CoordinateSystem cs = db.GetCoordinateSystem(projectProj, ProjectionSearchType.Enhanced);//坐标系统
if (cs != null)
{
db.ReadDialects(cs);
foreach (string dialect in cs.Dialects)
{
MapWinGIS.GeoProjection projTemp = new MapWinGIS.GeoProjection();
if (!projTemp.ImportFromAutoDetect(dialect))
{
continue;
}
if (layerProj.get_IsSame(projTemp))
{
isSame = true;
break;
}
}
}
}
}
// 投影中可以包含的文件的后缀名，让我们试着用正确的后缀搜索文件
if (!isSame)
{
if (CoordinateTransformation.SeekSubstituteFile(layer, projectProj, out newLayer))
{
m_report.AddFile(layer.Filename, layer.Projection.Name, ProjectionOperaion.Substituted, newLayer.Filename);
return(TestingResult.Substituted);
}
}
if (!layer.Projection.IsEmpty)
{
if (projectProj.IsEmpty)
{
// 层具有投影，项目没有；分配到投影，不提示用户
// 让我们找个众所周知的投影与EPSG编码
ProjectionDatabase db = m_mapWin.ProjectionDatabase as ProjectionDatabase;
if (db != null)
{
CoordinateSystem cs = db.GetCoordinateSystem(layerProj, ProjectionSearchType.UseDialects);
if (cs != null)
{
MapWinGIS.GeoProjection proj = new MapWinGIS.GeoProjection();
if (proj.ImportFromEPSG(cs.Code))
{
layerProj = proj;
}
}
}
m_mapWin.Project.GeoProjection = layerProj;
return(TestingResult.Ok);
}
else if (isSame)
{
// 相同投影
return(TestingResult.Ok);
}
else
{
// 用户必须被提示
if (!m_usePreviousAnswerMismatch && !m_mapWin.ApplicationInfo.NeverShowProjectionDialog)
{
bool dontShow = false;
bool useForOthers = false;
ArrayList list = new ArrayList();
list.Add("Ignore mismatch");
list.Add("Reproject file");
//list.Add("Skip file");
// PM 2013-05-03:
list.Add("Don't load the layer");
frmProjectionMismatch form = new frmProjectionMismatch((ProjectionDatabase)m_mapWin.ProjectionDatabase);
int choice = form.ShowProjectionMismatch(list, (int)m_mapWin.ApplicationInfo.ProjectionMismatchBehavior,
projectProj, layer.Projection, out useForOthers, out dontShow);
form.Dispose();
if (choice == -1)
{
return(TestingResult.CancelOperation);
}
m_usePreviousAnswerMismatch = useForOthers;
m_mapWin.ApplicationInfo.ProjectionMismatchBehavior = (ProjectionMismatchBehavior)choice;
m_mapWin.ApplicationInfo.NeverShowProjectionDialog = dontShow;
}
MapWinGIS.Interfaces.ProjectionMismatchBehavior behavior = m_mapWin.ApplicationInfo.ProjectionMismatchBehavior;
switch (behavior)
{
case ProjectionMismatchBehavior.Reproject:
TestingResult result = CoordinateTransformation.ReprojectLayer(layer, out newLayer, projectProj, m_report);
if (result == TestingResult.Ok || result == TestingResult.Substituted)
{
ProjectionOperaion oper = result == TestingResult.Ok ? ProjectionOperaion.Reprojected : ProjectionOperaion.Substituted;
string newName = newLayer == null ? "" : newLayer.Filename;
m_report.AddFile(layer.Filename, layer.Projection.Name, oper, newName);
return(newName == layer.Filename ? TestingResult.Ok : TestingResult.Substituted);
}
else
{
m_report.AddFile(layer.Filename, layer.Projection.Name, ProjectionOperaion.FailedToReproject, "");
return(TestingResult.Error);
}
case ProjectionMismatchBehavior.IgnoreMismatch:
m_report.AddFile(layer.Filename, layer.Projection.Name, ProjectionOperaion.MismatchIgnored, "");
return(TestingResult.Ok);
case ProjectionMismatchBehavior.SkipFile:
m_report.AddFile(layer.Filename, layer.Projection.Name, ProjectionOperaion.Skipped, "");
return(TestingResult.SkipFile);
}
}
}
else if (!projectProj.IsEmpty) // 图层投影是空的
{
bool projectProjectionExists = !projectProj.IsEmpty;
// 用户必须被提示
if (!m_usePreviousAnswerAbsence && !m_mapWin.ApplicationInfo.NeverShowProjectionDialog)
{
bool dontShow = false;
bool useForOthers = false;
ArrayList list = new ArrayList();
// 当在投影第一变体应排除
int val = projectProjectionExists ? 0 : 1;
if (projectProjectionExists)
{
// PM 2013-05-03:
//list.Add("Assign projection from project");
list.Add("Use the project's projection");
}
// list.Add("Ignore the absence");
// list.Add("Skip the file");
list.Add("Ignore the missing of projection file");
list.Add("Don't load the layer");
frmProjectionMismatch form = new frmProjectionMismatch((ProjectionDatabase)m_mapWin.ProjectionDatabase);
int choice = form.ShowProjectionAbsence(list, (int)m_mapWin.ApplicationInfo.ProjectionAbsenceBehavior - val, projectProj, out useForOthers, out dontShow);
form.Dispose();
if (choice == -1)
{
return(TestingResult.CancelOperation);
}
choice += val;
m_usePreviousAnswerAbsence = useForOthers;
m_mapWin.ApplicationInfo.ProjectionAbsenceBehavior = (ProjectionAbsenceBehavior)choice;
m_mapWin.ApplicationInfo.NeverShowProjectionDialog = dontShow;
}
// 当在项目没有投影，它不能分配层
ProjectionAbsenceBehavior behavior = m_mapWin.ApplicationInfo.ProjectionAbsenceBehavior;
if (!projectProjectionExists && m_mapWin.ApplicationInfo.ProjectionAbsenceBehavior == ProjectionAbsenceBehavior.AssignFromProject)
{
behavior = ProjectionAbsenceBehavior.IgnoreAbsence;
}
switch (behavior)
{
case ProjectionAbsenceBehavior.AssignFromProject:
m_report.AddFile(layer.Filename, layer.Projection.Name, ProjectionOperaion.Assigned, "");
layer.Projection = projectProj;
return(TestingResult.Ok);
case ProjectionAbsenceBehavior.IgnoreAbsence:
m_report.AddFile(layer.Filename, layer.Projection.Name, ProjectionOperaion.AbsenceIgnored, "");
return(TestingResult.Ok);
case ProjectionAbsenceBehavior.SkipFile:
m_report.AddFile(layer.Filename, layer.Projection.Name, ProjectionOperaion.Skipped, "");
return(TestingResult.SkipFile);
}
}
else
{
// 层没有投影，项目也没有，不在这里
}
System.Diagnostics.Debug.Print("Invalid result in projection tester");
return(TestingResult.Ok);
}
/// <summary>
/// Does the reprojection work
/// </summary>
private void DoReprojection(IEnumerable <string> filenames, MapWinGIS.GeoProjection projection, bool inPlace)
{
frmTesterReport report = new frmTesterReport();
report.InitProgress(projection);
List <string> files = new List <string>();
int count = 0; // number of successfully reprojected shapefiles
foreach (string filename in filenames)
{
LayerSource layer = new LayerSource(filename);
LayerSource layerNew = null;
if (projection.get_IsSame(layer.Projection))
{
report.AddFile(layer.Filename, projection.Name, ProjectionOperaion.SameProjection, "");
files.Add(layer.Filename);
}
else
{
TestingResult result = CoordinateTransformation.ReprojectLayer(layer, out layerNew, projection, report);
if (result == TestingResult.Ok || result == TestingResult.Substituted)
{
ProjectionOperaion oper = result == TestingResult.Ok ? ProjectionOperaion.Reprojected : ProjectionOperaion.Substituted;
string newName = layerNew == null ? "" : layerNew.Filename;
report.AddFile(layer.Filename, layer.Projection.Name, oper, newName);
files.Add(newName == "" ? layer.Filename : newName);
count++;
}
else
{
ProjectionOperaion operation = result == TestingResult.Error ? ProjectionOperaion.FailedToReproject : ProjectionOperaion.Skipped;
report.AddFile(layer.Filename, layer.Projection.Name, ProjectionOperaion.Skipped, "");
}
}
layer.Close();
if (layerNew != null)
{
layerNew.Close();
}
}
report.ShowReport(projection, "Reprojection results:", ReportType.Loading);
IEnumerable <string> names = m_mapWin.Layers.Select(l => l.FileName);
names = files.Except(names);
if (count > 0)
{
if (projection.get_IsSame(m_mapWin.Project.GeoProjection))
{
if (names.Count() > 0)
{
if (MessageBox.Show("Do you want to add layers to the project?", m_mapWin.ApplicationInfo.ApplicationName,
MessageBoxButtons.YesNo, MessageBoxIcon.Question) == DialogResult.Yes)
{
m_mapWin.Layers.StartAddingSession();
foreach (string filename in names)
{
m_mapWin.Layers.Add(filename);
}
m_mapWin.Layers.StopAddingSession();
}
}
else
{
MessageBox.Show("No files to add to the map.",
m_mapWin.ApplicationInfo.ApplicationName, MessageBoxButtons.OK, MessageBoxIcon.Information);
}
}
else
{
MessageBox.Show("Chosen projection is different from the project one. The layers can't be added to map.",
m_mapWin.ApplicationInfo.ApplicationName, MessageBoxButtons.OK, MessageBoxIcon.Information);
}
}
else
{
MessageBox.Show("No files to add to the map.",
m_mapWin.ApplicationInfo.ApplicationName, MessageBoxButtons.OK, MessageBoxIcon.Information);
}
}
/// <summary>
/// Reprojects a shapefile
/// </summary>
/// <param name="grid">Shapefile object to reproject</param>
/// <param name="inPlace">Whether reprojected file should replace the initial one</param>
/// <returns>True on success and false otherwise</returns>
public static TestingResult Reproject(MapWinGIS.Shapefile sfSource, out MapWinGIS.Shapefile sfNew, MapWinGIS.GeoProjection projection, frmTesterReport report)
{
sfNew = null;
string origFilename = sfSource.Filename;
MapWinGIS.Shapefile sf = new MapWinGIS.Shapefile();
int count = 0;
LayerSource obj = new LayerSource(sfSource);
LayerSource objNew = null;
if (CoordinateTransformation.SeekSubstituteFile(obj, projection, out objNew))
{
sfNew = objNew.Shapefile;
return(TestingResult.Substituted);
}
string newFilename = CoordinateTransformation.FilenameWithProjectionSuffix(origFilename, sfSource.GeoProjection, projection);
newFilename = CoordinateTransformation.GetSafeNewName(newFilename);
// settings callback
MapWinGIS.ICallback callback = sfSource.GlobalCallback;
if (report != null)
{
sfSource.GlobalCallback = report;
if (!report.Visible)
{
report.InitProgress(projection);
}
report.ShowFilename(sfSource.Filename);
}
// doing the job
sf = sfSource.Reproject(projection, ref count);
// restoring callback
if (report != null)
{
sfSource.GlobalCallback = callback;
report.ClearFilename();
}
if (sf != null && count == sfSource.NumShapes)
{
sf.GlobalCallback = sfSource.GlobalCallback;
bool result = sf.SaveAs(newFilename, null); // it doesn't close the editing mode
if (!result)
{
System.Diagnostics.Debug.Print("Error while saving reprojected shapefile: " + sf.get_ErrorMsg(sf.LastErrorCode));
}
sfNew = sf;
return(TestingResult.Ok);
}
//sf.Close();
return(TestingResult.Error);
}