/// <summary> /// Function to identify rendered and source dataset pixel values for one or more /// raster, mosaic and image service layers. /// </summary> /// <param name="mapPoint">Map point to identify pixel values.</param> public static async void CustomRasterIdentify(MapPoint mapPoint) { // Create a new list of popup pages. var popupContents = new List <PopupContent>(); // Create an empty string that will represent what goes into a popup page. string identifiedVal = ""; // Create the popup pages to show. await QueuedTask.Run(() => { // Check if there is an active map view. if (MapView.Active != null) { // Get the active map view. var mapView = MapView.Active; // Get the list of selected layers. IReadOnlyList <Layer> selectedLayerList = MapView.Active.GetSelectedLayers(); if (selectedLayerList.Count == 0) { // If no layers are selected fill the popup page with the appropriate message. // Note: you can use html tags to format the text. identifiedVal += "<p>No Layers selected. Please select one or more Raster, Mosaic or Image Service layers.</p>"; // Add the popup page to the list of pages. popupContents.Add(new PopupContent(identifiedVal, "Custom Raster Identify")); } else { // Iterate over the list of selected layers. foreach (Layer currentSelectedLayer in selectedLayerList) { #region Get a basic raster layer from the selected layer. BasicRasterLayer currentRasterLayer = null; if (currentSelectedLayer is MosaicLayer) { // If the current selected layer is a mosaic layer, MosaicLayer mosaicLayer = currentSelectedLayer as MosaicLayer; // Get the image sub-layer from the mosaic layer. This is a basic raster layer. currentRasterLayer = mosaicLayer.GetImageLayer() as BasicRasterLayer; } else if (currentSelectedLayer is BasicRasterLayer) { // If the current selected layer is a raster layer or image service layer, // both are already basic raster layers. currentRasterLayer = currentSelectedLayer as BasicRasterLayer; } else { // If the current selected layer is neither a mosaic nor a raster or image service layer, // fill the popup page with the appropriate message. identifiedVal += "<p>Selected layer is not a raster layer. Please select one or more Raster, Mosaic or Image Service layers.</p>"; // Add the popup page to the list of pages. popupContents.Add(new PopupContent(identifiedVal, "Custom Raster Identify for: " + currentSelectedLayer.Name)); continue; } #endregion #region Get the pixel value for the rendered raster. // Add the label for the rendered pixel value. identifiedVal += "<b>Rendered Pixel value: </b>"; // Get the raster from the current selected raster layer. Raster raster = currentRasterLayer.GetRaster(); // If the map spatial reference is different from the spatial reference of the raster, // set the map spatial reference on the raster. This will ensure the map points are // correctly reprojected to image points. if (mapView.Map.SpatialReference.Name != raster.GetSpatialReference().Name) { raster.SetSpatialReference(mapView.Map.SpatialReference); } // Convert the map point to be identified into an image point. Tuple <int, int> imagePoint = raster.MapToPixel(mapPoint.X, mapPoint.Y); if ((int)imagePoint.Item1 < 0 || (int)imagePoint.Item1 > raster.GetWidth() || (int)imagePoint.Item2 < 0 || (int)imagePoint.Item2 > raster.GetHeight()) { // If the point is outside the image, fill the pixel value with the appropriate message. identifiedVal += "Point is not within image. \n"; } else { // If the point is within the image. Iterate over the bands in the raster. for (int band = 0; band < raster.GetBandCount(); band++) { // Get the pixel value based on the band, column and row and add the // formatted pixel value to the popup page. if (band == 0) { identifiedVal += raster.GetPixelValue(band, (int)imagePoint.Item1, (int)imagePoint.Item2).ToString(); } else { identifiedVal += ", " + raster.GetPixelValue(band, (int)imagePoint.Item1, (int)imagePoint.Item2).ToString(); } } identifiedVal += "."; } // Add the rendered pixel value to the popup page contents. string htmlContent = "<p>" + identifiedVal + "</p>"; #endregion #region Get the pixel value for the source raster dataset. // Add the label for the source dataset pixel value. identifiedVal = "<b>Dataset Pixel value: </b>"; // Get the basic raster dataset from the raster. BasicRasterDataset basicRasterDataset = raster.GetRasterDataset(); if (!(basicRasterDataset is RasterDataset)) { // If the dataset is not a raster dataset, fill the pixel value with the appropriate message. identifiedVal += "Selected layer is not a Raster Layer. Please select one or more Raster, Mosaic or Image Service layers."; htmlContent += "<p>" + identifiedVal + "</p>"; popupContents.Add(new PopupContent(htmlContent, "Custom Raster Identify for " + currentSelectedLayer.Name)); } // Get the raster dataset. RasterDataset rasterDataset = basicRasterDataset as RasterDataset; // Create a full raster from the raster dataset. raster = rasterDataset.CreateFullRaster(); // If the map spatial reference is different from the spatial reference of the raster, // Set the map spatial reference on the raster. This will ensure the map points are // correctly reprojected to image points. if (mapView.Map.SpatialReference.Name != raster.GetSpatialReference().Name) { raster.SetSpatialReference(mapView.Map.SpatialReference); } // Convert the map point to be identified to an image point. imagePoint = raster.MapToPixel(mapPoint.X, mapPoint.Y); if ((int)imagePoint.Item1 < 0 || (int)imagePoint.Item1 > raster.GetWidth() || (int)imagePoint.Item2 < 0 || (int)imagePoint.Item2 > raster.GetHeight()) { // If the point is outside the image, fill the pixel value with the appropriate message. identifiedVal += "Point is not within image. \n"; } else { // If the point is within the image. Iterate over the bands in the raster. for (int band = 0; band < raster.GetBandCount(); band++) { // Get the pixel value based on the band, column and row and add the // formatted pixel value to the popup page. if (band == 0) { identifiedVal += raster.GetPixelValue(band, (int)imagePoint.Item1, (int)imagePoint.Item2).ToString(); } else { identifiedVal += ", " + raster.GetPixelValue(band, (int)imagePoint.Item1, (int)imagePoint.Item2).ToString(); } } identifiedVal += "."; } // Add the source dataset pixel value to the popup page contents. htmlContent += "<p>" + identifiedVal + "</p>"; #endregion // Add the popup page to the list of pages. popupContents.Add(new PopupContent(htmlContent, "Custom Raster Identify for " + currentSelectedLayer.Name)); // Reset identifiedVal = ""; } } } }); // Show custom popup with the list of popup pages created above. MapView.Active.ShowCustomPopup(popupContents); }
/// <summary> /// Mask raster pixels based on the rectangle given and save the output in the /// current project folder. /// </summary> /// <param name="geometry">Rectangle to use to mask raster pixels.</param> public static async void MaskRaster(Geometry geometry) { try { // Check if there is an active map view. if (MapView.Active != null) { // Get the active map view. var mapView = MapView.Active; // Get the list of selected layers. IReadOnlyList <Layer> selectedLayerList = mapView.GetSelectedLayers(); if (selectedLayerList.Count == 0) { // If no layers are selected show a message box with the appropriate message. MessageBox.Show("No Layers selected. Please select one Raster layer."); } else { // Get the most recently selected layer. Layer firstSelectedLayer = mapView.GetSelectedLayers().First(); if (firstSelectedLayer is RasterLayer) { // Working with rasters requires the MCT. await QueuedTask.Run(() => { #region Get the raster dataset from the currently selected layer // Get the raster layer from the selected layer. RasterLayer currentRasterLayer = firstSelectedLayer as RasterLayer; // Get the raster from the current selected raster layer. Raster inputRaster = currentRasterLayer.GetRaster(); // Get the basic raster dataset from the raster. BasicRasterDataset basicRasterDataset = inputRaster.GetRasterDataset(); if (!(basicRasterDataset is RasterDataset)) { // If the dataset is not a raster dataset, show a message box with the appropriate message. MessageBox.Show("No Raster Layers selected. Please select one Raster layer."); return; } // Get the input raster dataset from the basic raster dataset. RasterDataset rasterDataset = basicRasterDataset as RasterDataset; #endregion #region Save a copy of the raster dataset in the project folder and open it // Create a full raster from the input raster dataset. inputRaster = rasterDataset.CreateFullRaster(); // Setup the paths and name of the output file and folder inside the project folder. string ouputFolderName = "MaskedOuput"; string outputFolder = Path.Combine(Project.Current.HomeFolderPath, ouputFolderName);; string outputName = "MaskedRaster.tif"; // Delete the output directory if it exists and create it. // Note: You will need write access to the project directory for this sample to work. if (Directory.Exists(outputFolder)) { Directory.Delete(outputFolder, true); } Directory.CreateDirectory(outputFolder); // Create a new file system connection path to open raster datasets using the output folder path. FileSystemConnectionPath outputConnectionPath = new FileSystemConnectionPath( new System.Uri(outputFolder), FileSystemDatastoreType.Raster); // Create a new file system data store for the connection path created above. FileSystemDatastore outputFileSytemDataStore = new FileSystemDatastore(outputConnectionPath); // Create a new raster storage definition. RasterStorageDef rasterStorageDef = new RasterStorageDef(); // Set the pyramid level to 0 meaning no pyramids will be calculated. This is required // because we are going to change the pixels after we save the raster dataset and if the // pyramids are calculated prior to that, the pyramids will be incorrect and will have to // be recalculated. rasterStorageDef.SetPyramidLevel(0); // Save a copy of the raster using the file system data store and the raster storage definition. inputRaster.SaveAs(outputName, outputFileSytemDataStore, "TIFF", rasterStorageDef); // Open the raster dataset you just saved. rasterDataset = OpenRasterDataset(outputFolder, outputName); // Create a full raster from it so we can modify pixels. Raster outputRaster = rasterDataset.CreateFullRaster(); #endregion #region Get the Min/Max Row/Column to mask // If the map spatial reference is different from the spatial reference of the input raster, // set the map spatial reference on the input raster. This will ensure the map points are // correctly reprojected to image points. if (mapView.Map.SpatialReference.Name != inputRaster.GetSpatialReference().Name) { inputRaster.SetSpatialReference(mapView.Map.SpatialReference); } // Use the MapToPixel method of the input raster to get the row and column values for the // points of the rectangle. Tuple <int, int> tlcTuple = inputRaster.MapToPixel(geometry.Extent.XMin, geometry.Extent.YMin); Tuple <int, int> lrcTuple = inputRaster.MapToPixel(geometry.Extent.XMax, geometry.Extent.YMax); int minCol = (int)tlcTuple.Item1; int minRow = (int)tlcTuple.Item2; int maxCol = (int)lrcTuple.Item1; int maxRow = (int)lrcTuple.Item2; // Ensure the min's are less than the max's. if (maxCol < minCol) { int temp = maxCol; maxCol = minCol; minCol = temp; } if (maxRow < minRow) { int temp = maxRow; maxRow = minRow; minRow = temp; } // Ensure the mins and maxs are within the raster. minCol = (minCol < 0) ? 0 : minCol; minRow = (minRow < 0) ? 0 : minRow; maxCol = (maxCol > outputRaster.GetWidth()) ? outputRaster.GetWidth() : maxCol; maxRow = (maxRow > outputRaster.GetHeight()) ? outputRaster.GetHeight() : maxRow; #endregion #region Mask pixels based on the rectangle drawn by the user // Calculate the width and height of the pixel block to create. int pbWidth = maxCol - minCol; int pbHeight = maxRow - minRow; // Check to see if the output raster can be edited. if (!outputRaster.CanEdit()) { // If not, show a message box with the appropriate message. MessageBox.Show("Cannot edit raster :("); return; } // Create a new pixel block from the output raster of the height and width calculated above. PixelBlock currentPixelBlock = outputRaster.CreatePixelBlock(pbWidth, pbHeight); // Iterate over the bands of the output raster. for (int plane = 0; plane < currentPixelBlock.GetPlaneCount(); plane++) { // For each band, clear the pixel block. currentPixelBlock.Clear(plane); //Array noDataMask = currentPixelBlock.GetNoDataMask(plane, true); //for (int i = 0; i < noDataMask.GetLength(0); i++) // noDataMask.SetValue(Convert.ToByte(0), i); //currentPixelBlock.SetNoDataMask(plane, noDataMask); } // Write the cleared pixel block to the output raster dataset. outputRaster.Write(minCol, minRow, currentPixelBlock); // Refresh the properties of the output raster dataset. outputRaster.Refresh(); #endregion // Create a new layer from the masked raster dataset and add it to the map. LayerFactory.Instance.CreateLayer(new Uri(Path.Combine(outputFolder, outputName)), mapView.Map); // Disable the layer representing the original raster dataset. firstSelectedLayer.SetVisibility(false); }); } else { // If the selected layer is not a raster layer show a message box with the appropriate message. MessageBox.Show("No Raster layers selected. Please select one Raster layer."); } } } } catch (Exception exc) { MessageBox.Show("Exception caught in MaskRaster: " + exc.Message); } }