Latest Contributions by LewisTrotter

Hello, The recent update Pro 3.2 (and GDAL) has broken the ability for GDAL to open WCS urls. We rely on this for some plugins. Running below code in Pro <= 3.1 works perfectly (returns a dataset), but results in ERROR 4 for Pro 3.2. I tried below on the latets version of GDAL outside of Pro and it also works fine. Any advice? Thanks! from osgeo import gdal url = 'https://ows.dea.ga.gov.au/wcs?service=WCS&VERSION=1.0.0&REQUEST=GetCoverage&COVERAGE=ga_ls5t_nbart_gm_cyear_3&MEASUREMENTS=blue,green,red,sdev&TIME=1991-01-01&BBOX=-1491464.3902944827,-3760705.2463235464,-1485657.8995810281,-3753979.4918812234&CRS=EPSG:3577&RESX=30&RESY=30&FORMAT=GeoTIFF' ds = gdal.Open(url, gdal.GA_ReadOnly) print(ds)   ... View more

Hello, I recently upgraded to Pro 3.2 and I'm now experiencing unexpected behaviour with the Register Raster tool in certain circumstances. I have a function hat registers numerous rasters in a loop to a reference. Every so often a raster identical to the reference is put through the tool. When this happens, the registered raster is sent kilometres away - but only when both rasters are identical. Is this expected behaviour? Thanks ... View more

Hello all, I am trying to get some basic code that uses the multiprocessing module to work in a Pro 3.1 toolbox but it just won't play nice. The code:   import time import arcpy import multiprocessing as mp def work_log(work_data): arcpy.AddMessage(" Process %s waiting %s seconds" % (work_data[0], work_data[1])) time.sleep(int(work_data[1])) arcpy.AddMessage(" Process %s Finished." % work_data[0]) def pool_handler(work): p = mp.Pool(2) p.map(work_log, work) if __name__ == '__main__': work = (["A", 5], ["B", 2], ["C", 1], ["D", 3]) pool_handler(work)     If I run this in an .atnx or .pyt I get the same error:      PicklingError: Can't pickle <function work_log at 0x00000134B287D430>: attribute lookup work_log on __main__ failed     I also get 2 new instances of Pro open automatically if I run in a .pyt toolbox, I assume due to lack of __main__. Does anyone have any advice on how to get something like this running in a Pro toolbox (preferably .pyt)? I came across a similar question here https://community.esri.com/t5/python-questions/using-multiprocessing/m-p/1280108 but it has never been answered. Thanks! ... View more

Hello, I have a full Landsat scene that contains the usual bands (RGB NIR etc) as well as a seperate QA cloud mask band with the exact same dimensions and pixel size. I am trying to set all pixels in the full Landsat scene to NoData (-1) wherever a value of 0, 2, or 3 exist in the cloud mask. Are there any efficient approaches to this in the SDK? I came across the Mask Function (not the GP tool) in Pro, but from what I understand these are not available in the SDK. I also do not have access to Spatial Analyst to use Extract By Mask. I've currently implemented a a for loop that basically iterates the cloud mask and whenever the value is found, updates the same pixel in the full Landsat scene, but it seems suboptimal to me. I was looking at the GetNoDataMask and SetNoDataMask methods but as far as I'm aware my raw tiffs do not have these encoded. Cheers ... View more

Hello all, I was wonder if it is possible to NuGet GDAL for use in a .NET Addin module?  I got the recommended GDAL 3.5.0 and GDAL.Native 3.5.0 NuGet libraries via Visual Studio 2022 and everything came in fine. However, when running the required GdalConfiguration.ConfigureGdal() command, I get an exception: "The type initializer for 'OSGeo.GDAL.GdalPINVOKE' threw an exception." I've looked around and there's not a heap of information on it. If I create a new C# project without ArcGIS SDK I can use GDAL without any issues, so I am assuming this is a ArcGIS Pro conflict? Any tips appreciated. ... View more

Hello all, I've been trying to implement a progress bar on a dock pane that increments as a function in a separate class downloads a list of GeoTiffs using HttpClient. I've successfully got this to work using the approach provided in the code sample Overlay3D. However, unlike this example, my list of tasks I want to increment occurs in a separate function outside of the dock pane, and I'm not sure how to pass the information from that function back to the dock pane. Any advice is welcome. I have provided a cut down code example below of my current setup, the comments should show where I am getting stuck:   XAML progress bar, with two bindings (MaxProgressValue and Progress Value) back to dock pane view model.     <!-- Progress bar--> <ProgressBar Name="PrgProgressBar" Grid.Column="0" Margin="0,0,5,0" Height = "10" HorizontalAlignment="Stretch" Minimum="1" Maximum="{Binding Path=MaxProgressValue, Mode=OneWay}" Value="{Binding Path=ProgressValue, Mode=OneWay}" IsIndeterminate="False"> </ProgressBar>     ViewModel related to progress bar, with a button to call a function to download geotiffs.     // Getter, Setter for progress bar private double _maxProgressValue = 100; public double MaxProgressValue { get { return _maxProgressValue; } set { SetProperty(ref _maxProgressValue, value, () => MaxProgressValue); } } // Update progress bar method (based on Overlay3D) private int _iProgressValue = -1; private int _iMaxProgressValue = -1; private void ProgressUpdate(int iProgressValue, int iMaxProgressValue) { if (System.Windows.Application.Current.Dispatcher.CheckAccess()) { if (_iMaxProgressValue != iMaxProgressValue) { MaxProgressValue = iMaxProgressValue; } if (_iProgressValue != iProgressValue) { ProgressValue = iProgressValue; } } else { ProApp.Current.Dispatcher.BeginInvoke(DispatcherPriority.Background, (Action)(() => { if (_iMaxProgressValue != iMaxProgressValue) { MaxProgressValue = iMaxProgressValue; } if (_iProgressValue != iProgressValue) { ProgressValue = iProgressValue; } })); } } // Button on dock panel to run download function public ICommand CmdRun { get { return new RelayCommand(async () => { // Excluded various things // ... // Initialise data structure STAC that contains urls and functions Stac stac = new Stac(...) // Reset progressor value and maximum value ProgressUpdate(1, stac.Result.Features.Count); // Start queued task and download geotifs await QueuedTask.Run(() => { // Begin download tifs loop // How do I monitor progession? // Pass a progressor variable via IProgress? await stac.DownloadFeaturesAsync(assets, outputFolder); }); } }      Stac class for DownloadFeatureAsync.     public async Task DownloadFeaturesAsync(List<string> assetNames, string outputFolder) { // Some uninteresting work... // ... // Iterate dates in STAC query and download tif to folder foreach (string date in dates) { // Open a HTTP client var client = new HttpClient(); // Create a list of downloader tasks int i = 0; // test var tasks = new List<Task>(); foreach (var item in items) { string date = item.Key; string url = item.Value[0]; string filename = item.Value[1]; tasks.Add(Task.Run(async () => { // Do some things... // Query URL and ensure success var response = await client.GetAsync(url); response.EnsureSuccessStatusCode(); // Download file to output folder using (FileStream fs = new FileStream(filepath, FileMode.CreateNew)) { await response.Content.CopyToAsync(fs); } // Increment progress bar... how...?! // ... })); } // Run all download tasks await Task.WhenAll(tasks); }     ... View more