Openw, un, output_image_group.filename, /get_lun This_file = output_base + output_image_group.extension this_tag = "_" + timestamp() + "_" + landtrendr_version()įor i = 0, n_elements(output_image_group)-1 do begin Output_image_group.DATA = "Segment Mean Spectral" Output_image_group.extension = '_segmean.bsq' Output_image_group.DATA = "Source Spectral Stack" Output_image_group.extension = '_source.bsq' Output_image_group.extension = '_segmse.bsq' Output_image_group.DATA = "Segment Statistics" Output_image_group.extension = '_stats.bsq' Output_image_group.DATA = "Fitted Spectral Stack" Output_image_group.n_layers = n_elements(years) Output_image_group.extension = '_fitted.bsq' cloud-mosaicking, then we need to compensate for in a given year are provided with the idea of doing 6/27/08 for fitted image, the output will be the
Output_image_group.DATA = "Trajectory Disturbance/Recovery"
Output_image_group.extension = '_distrec.bsq' Output_image_group.DATA = "Segment Duration" Output_image_group.extension = '_durs.bsq' Output_image_group.DATA = "Segment Magnitude" Output_image_group.n_layers = max_segments Output_image_group.extension = '_mags.bsq' Output_image_group.DATA = "Vertex Spectral Value" Output_image_group.extension = '_vertvals.bsq' Output_image_group.n_layers = max_segments+1 Output_image_group.extension = '_vertyrs.bsq' Output_image_group = replicate(output_image_base, 10) Same number of years as all of the inputs, 3: scaled ratio between disturbance and recovery mags: the magnitude of the segment change between two vertices vertvals: the values of the input band at those years Ok = define_chunks3(subset, pixsize, max_pixels_per_chunk, kernelsize) Zot_img, image_info.image_file, hdr, img, subset=subset, /hdronly Print, "process_cm_biomass_image.pro: Image in If file_exists(image_info.image_file) eq 0 then begin to get the pixel size of the image, assume that all are the same If n_elements(max_pixels_per_chunk) eq 0 then $ SET UP THE OUTPUT FILE AND CHUNK INFORMATION Vertexcountovershoot = diag_info.vertexcountovershootīestmodelproportion = diag_info.bestmodelproportionĮnd else begin if this image has not been set up before, then Recovery_threshold = diag_info.recovery_thresholdĭistweightfactor = diag_info.distweightfactor Pixels_per_chunk = diag_info.pixels_per_chunk Output_image_group = diag_info.output_image_group Print, 'This file has already had processing done on it' process, and needs to just be picked up again see if it has a save file to indicate that it was already in first, check to see if the output image already exists. Print, 'bestmodelproportion cannot exceed 1. Print, 'Vertexcountovershoot cannot exceed 3' Vertexcountovershoot = run_params.vertexcountovershootīestmodelproportion = run_params.bestmodelproportion Recovery_threshold=run_params.recovery_thresholdĭistweightfactor = run_params.distweightfactor Index = run_params.index string with name of indexīackground_val = run_params.background_valĭesawtooth_val = run_sawtooth_valįix_doy_effect = run_params.fix_doy_effect copy over just to make it work with historical code. The filename, directory name, or volume label syntax is incorrect.įunction process_tbcd_chunks, run_params, progressbaryesno I don't know if there is enough background provided here for any of you to help my get past this error, but any help you may be able to provide would be greatly appreciated. I've made it through a good number of modules, but have hit a breaking point.
Idl openw code#
I have virtually no IDL experience, but was told that since the code is already written and has been tested, I should have very few hiccups. I'm attempting to run the LandTrendr model in IDL, and have hit a snag.