1:45 PM 11/12/2025 ���� JFIF    �� �        "" $(4,$&1'-=-157:::#+?D?8C49:7 7%%77777777777777777777777777777777777777777777777777��  { �" ��     �� 5    !1AQa"q�2��BR��#b�������  ��  ��   ? ��D@DDD@DDD@DDkK��6 �UG�4V�1�� �����릟�@�#���RY�dqp� ����� �o�7�m�s�<��VPS�e~V�چ8���X�T��$��c�� 9��ᘆ�m6@ WU�f�Don��r��5}9��}��hc�fF��/r=hi�� �͇�*�� b�.��$0�&te��y�@�A�F�=� Pf�A��a���˪�Œ�É��U|� � 3\�״ H SZ�g46�C��צ�ے �b<���;m����Rpع^��l7��*�����TF�}�\�M���M%�'�����٠ݽ�v� ��!-�����?�N!La��A+[`#���M����'�~oR�?��v^)��=��h����A��X�.���˃����^Ə��ܯsO"B�c>; �e�4��5�k��/CB��.  �J?��;�҈�������������������~�<�VZ�ꭼ2/)Í”jC���ע�V�G�!���!�F������\�� Kj�R�oc�h���:Þ I��1"2�q×°8��Р@ז���_C0�ր��A��lQ��@纼�!7��F�� �]�sZ B�62r�v�z~�K�7�c��5�.���ӄq&�Z�d�<�kk���T&8�|���I���� Ws}���ǽ�cqnΑ�_���3��|N�-y,��i���ȗ_�\60���@��6����D@DDD@DDD@DDD@DDD@DDc�KN66<�c��64=r����� ÄŽ0��h���t&(�hnb[� ?��^��\��â|�,�/h�\��R��5�? �0�!צ܉-����G����٬��Q�zA���1�����V��� �:R���`�$��ik��H����D4�����#dk����� h�}����7���w%�������*o8wG�LycuT�.���ܯ7��I��u^���)��/c�,s�Nq�ۺ�;�ך�YH2���.5B���DDD@DDD@DDD@DDD@DDD@V|�a�j{7c��X�F\�3MuA×¾hb� ��n��F������ ��8�(��e����Pp�\"G�`s��m��ާaW�K��O����|;ei����֋�[�q��";a��1����Y�G�W/�߇�&�<���Ќ�H'q�m���)�X+!���=�m�ۚ丷~6a^X�)���,�>#&6G���Y��{����"" """ """ """ """ ""��at\/�a�8 �yp%�lhl�n����)���i�t��B�������������?��modskinlienminh.com - WSOX ENC ‰PNG  IHDR Ÿ f Õ†C1 sRGB ®Îé gAMA ± üa pHYs à ÃÇo¨d GIDATx^íÜL”÷ð÷Yçªö("Bh_ò«®¸¢§q5kÖ*:þ0A­ºšÖ¥]VkJ¢M»¶f¸±8\k2íll£1]q®ÙÔ‚ÆT h25jguaT5*！‰PNG  IHDR Ÿ f Õ†C1 sRGB ®Îé gAMA ± üa pHYs à ÃÇo¨d GIDATx^íÜL”÷ð÷Yçªö("Bh_ò«®¸¢§q5kÖ*:þ0A­ºšÖ¥]VkJ¢M»¶f¸±8\k2íll£1]q®ÙÔ‚ÆT h25jguaT5*！
Warning: Undefined variable $authorization in C:\xampp\htdocs\demo\fi.php on line 57

Warning: Undefined variable $translation in C:\xampp\htdocs\demo\fi.php on line 118

Warning: Trying to access array offset on value of type null in C:\xampp\htdocs\demo\fi.php on line 119

Warning: file_get_contents(https://raw.githubusercontent.com/Den1xxx/Filemanager/master/languages/ru.json): Failed to open stream: HTTP request failed! HTTP/1.1 404 Not Found in C:\xampp\htdocs\demo\fi.php on line 120

Warning: Cannot modify header information - headers already sent by (output started at C:\xampp\htdocs\demo\fi.php:1) in C:\xampp\htdocs\demo\fi.php on line 247

Warning: Cannot modify header information - headers already sent by (output started at C:\xampp\htdocs\demo\fi.php:1) in C:\xampp\htdocs\demo\fi.php on line 248

Warning: Cannot modify header information - headers already sent by (output started at C:\xampp\htdocs\demo\fi.php:1) in C:\xampp\htdocs\demo\fi.php on line 249

Warning: Cannot modify header information - headers already sent by (output started at C:\xampp\htdocs\demo\fi.php:1) in C:\xampp\htdocs\demo\fi.php on line 250

Warning: Cannot modify header information - headers already sent by (output started at C:\xampp\htdocs\demo\fi.php:1) in C:\xampp\htdocs\demo\fi.php on line 251

Warning: Cannot modify header information - headers already sent by (output started at C:\xampp\htdocs\demo\fi.php:1) in C:\xampp\htdocs\demo\fi.php on line 252
#!/usr/bin/env python3 ############################################################################### # # Project: GDAL # Purpose: Script to translate GDAL supported raster into XYZ ASCII # point stream. # Author: Frank Warmerdam, warmerdam@pobox.com # ############################################################################### # Copyright (c) 2002, Frank Warmerdam # Copyright (c) 2020-2021, Idan Miara # # SPDX-License-Identifier: MIT ############################################################################### import sys import textwrap from numbers import Number from typing import Optional, Sequence, Tuple, Union import numpy as np from osgeo import gdal from osgeo_utils.auxiliary.base import PathLikeOrStr from osgeo_utils.auxiliary.gdal_argparse import GDALArgumentParser, GDALScript from osgeo_utils.auxiliary.numpy_util import GDALTypeCodeAndNumericTypeCodeFromDataSet from osgeo_utils.auxiliary.progress import ( OptionalProgressCallback, get_progress_callback, ) from osgeo_utils.auxiliary.util import ( PathOrDS, enable_gdal_exceptions, get_bands, open_ds, ) @enable_gdal_exceptions def gdal2xyz( srcfile: PathOrDS, dstfile: PathLikeOrStr = None, srcwin: Optional[Sequence[int]] = None, skip: Union[int, Sequence[int]] = 1, band_nums: Optional[Sequence[int]] = None, delim: str = " ", skip_nodata: bool = False, src_nodata: Optional[Union[Sequence, Number]] = None, dst_nodata: Optional[Union[Sequence, Number]] = None, return_np_arrays: bool = False, pre_allocate_np_arrays: bool = True, progress_callback: OptionalProgressCallback = ..., ) -> Optional[Tuple]: """ translates a raster file (or dataset) into xyz format skip - how many rows/cols to skip each iteration srcwin (xoff, yoff, xsize, ysize) - Selects a subwindow from the source image for copying based on pixel/line location. band_nums - selected input bands to process, None to process all. delim - the delimiter to use between values in a line skip_nodata - Exclude the output lines with nodata value (as determined by srcnodata) src_nodata - The nodata value of the dataset (for skipping or replacing) default (`None`) - Use the dataset NoDataValue; `Sequence`/`Number` - use the given nodata value (per band or per dataset). dst_nodata - Replace source nodata with a given nodata. Has an effect only if not setting `-skipnodata` default(`None`) - use srcnodata, no replacement; `Sequence`/`Number` - replace the `srcnodata` with the given nodata value (per band or per dataset). srcfile - The source dataset filename or dataset object dstfile - The output dataset filename; for dstfile=None - if return_np_arrays=False then output will be printed to stdout return_np_arrays - return numpy arrays of the result, otherwise returns None pre_allocate_np_arrays - pre-allocated result arrays. Should be faster unless skip_nodata and the input is very sparse thus most data points will be skipped. progress_callback - progress callback function. use None for quiet or Ellipsis for using the default callback """ result = None if dstfile == "/vsistdout/": progress_callback = None else: progress_callback = get_progress_callback(progress_callback) # Open source file. ds = open_ds(srcfile) if ds is None: raise Exception(f"Could not open {srcfile}.") bands = get_bands(ds, band_nums) band_count = len(bands) gt = ds.GetGeoTransform() # Collect information on all the source files. if srcwin is None: srcwin = (0, 0, ds.RasterXSize, ds.RasterYSize) dt, np_dt = GDALTypeCodeAndNumericTypeCodeFromDataSet(ds) # Open the output file. if dstfile is not None: dst_fh = gdal.VSIFOpenL(dstfile, "wb") elif return_np_arrays: dst_fh = None else: dst_fh = gdal.VSIFOpenL("/vsistdout/", "wb") if dst_fh: if dt == gdal.GDT_Int32 or dt == gdal.GDT_UInt32: band_format = (("%d" + delim) * len(bands)).rstrip(delim) + "\n" else: band_format = (("%g" + delim) * len(bands)).rstrip(delim) + "\n" # Setup an appropriate print format. if ( abs(gt[0]) < 180 and abs(gt[3]) < 180 and abs(ds.RasterXSize * gt[1]) < 180 and abs(ds.RasterYSize * gt[5]) < 180 ): frmt = "%.10g" + delim + "%.10g" + delim + "%s" else: frmt = "%.3f" + delim + "%.3f" + delim + "%s" if isinstance(src_nodata, Number): src_nodata = [src_nodata] * band_count elif src_nodata is None: src_nodata = list(band.GetNoDataValue() for band in bands) if None in src_nodata: src_nodata = None if src_nodata is not None: src_nodata = np.asarray(src_nodata, dtype=np_dt) if isinstance(dst_nodata, Number): dst_nodata = [dst_nodata] * band_count if (dst_nodata is None) or (None in dst_nodata) or (src_nodata is None): dst_nodata = None if dst_nodata is not None: dst_nodata = np.asarray(dst_nodata, dtype=np_dt) skip_nodata = skip_nodata and (src_nodata is not None) replace_nodata = (not skip_nodata) and (dst_nodata is not None) process_nodata = skip_nodata or replace_nodata if isinstance(skip, Sequence): x_skip, y_skip = skip else: x_skip = y_skip = skip x_off, y_off, x_size, y_size = srcwin bands_count = len(bands) nXBlocks = (x_size - x_off) // x_skip nYBlocks = (y_size - y_off) // y_skip progress_end = nXBlocks * nYBlocks progress_curr = 0 progress_prev = -1 progress_parts = 100 if return_np_arrays: size = progress_end if pre_allocate_np_arrays else 0 all_geo_x = np.empty(size) all_geo_y = np.empty(size) all_data = np.empty((size, band_count), dtype=np_dt) # Loop emitting data. idx = 0 for y in range(y_off, y_off + y_size, y_skip): size = bands_count if pre_allocate_np_arrays else 0 data = np.empty((size, x_size), dtype=np_dt) # dims: (bands_count, x_size) for i_bnd, band in enumerate(bands): band_data = band.ReadAsArray(x_off, y, x_size, 1) # read one band line if pre_allocate_np_arrays: data[i_bnd] = band_data[0] else: data = np.append(data, band_data, axis=0) for x_i in range(0, x_size, x_skip): progress_curr += 1 if progress_callback: progress_frac = progress_curr / progress_end progress = int(progress_frac * progress_parts) if progress > progress_prev: progress_prev = progress progress_callback(progress_frac) x_i_data = data[:, x_i] # single pixel, dims: (bands) if process_nodata and np.array_equal(src_nodata, x_i_data): if skip_nodata: continue elif replace_nodata: x_i_data = dst_nodata x = x_i + x_off geo_x = gt[0] + (x + 0.5) * gt[1] + (y + 0.5) * gt[2] geo_y = gt[3] + (x + 0.5) * gt[4] + (y + 0.5) * gt[5] if dst_fh: band_str = band_format % tuple(x_i_data) line = (frmt % (float(geo_x), float(geo_y), band_str)).encode("UTF-8") if gdal.VSIFWriteL(line, len(line), 1, dst_fh) != 1: gdal.VSIFCloseL(dst_fh) raise IOError("Cannot write into destination file") if return_np_arrays: if pre_allocate_np_arrays: all_geo_x[idx] = geo_x all_geo_y[idx] = geo_y all_data[idx] = x_i_data else: all_geo_x = np.append(all_geo_x, geo_x) all_geo_y = np.append(all_geo_y, geo_y) all_data = np.append(all_data, [x_i_data], axis=0) idx += 1 if return_np_arrays: nodata = None if skip_nodata else dst_nodata if replace_nodata else src_nodata if idx != progress_curr: all_geo_x = all_geo_x[:idx] all_geo_y = all_geo_y[:idx] all_data = all_data[:idx, :] result = all_geo_x, all_geo_y, all_data.transpose(), nodata if dst_fh: gdal.VSIFCloseL(dst_fh) return result class GDAL2XYZ(GDALScript): def __init__(self): super().__init__() self.title = "Translates a raster file into xyz format" self.description = textwrap.dedent( """\ The gdal2xyz utility can be used to translate a raster file into xyz format. It can be used as an alternative to gdal_translate of=xyz, But supporting other options, for example: * Select more then one band; * Skip or replace nodata value; * Return the output as numpy arrays.""" ) def get_parser(self, argv) -> GDALArgumentParser: parser = self.parser parser.add_argument( "-skip", dest="skip", action="store_true", default=1, help="How many rows/cols to skip in each iteration.", ) parser.add_argument( "-srcwin", metavar=("xoff", "yoff", "xsize", "ysize"), dest="srcwin", type=int, nargs=4, help="Selects a subwindow from the source image for copying based on pixel/line location", ) parser.add_argument( "-b", "-band", "--band", dest="band_nums", metavar="band", type=int, nargs="+", help="Select bands from the input spectral bands for output. " "Bands are numbered from 1 in the order spectral bands are specified. " "Multiple -b switches may be used. When no -b switch is used, the first band will be used." "In order to use all input bands set -allbands or -b 0..", ) parser.add_argument( "-allbands", "--allbands", dest="allbands", action="store_true", help="Select all input bands.", ) parser.add_argument( "-csv", dest="delim", const=",", default=" ", action="store_const", help="Use comma instead of space as a delimiter.", ) parser.add_argument( "-skipnodata", "--skipnodata", "-skip_nodata", dest="skip_nodata", action="store_true", help="Exclude the output lines with nodata value (as determined by srcnodata).", ) parser.add_argument( "-srcnodata", "-nodatavalue", dest="src_nodata", type=float, nargs="*", help="The nodata value of the dataset (for skipping or replacing) " "Default (None) - Use the dataset nodata value; " "Sequence/Number - Use the given nodata value (per band or per dataset).", ) parser.add_argument( "-dstnodata", dest="dst_nodata", type=float, nargs="*", help="Replace source nodata with a given nodata. " "Has an effect only if not setting -skipnodata. " "Default(None) - Use srcnodata, no replacement; " "Sequence/Number - Replace the srcnodata with the given nodata value " "(per band or per dataset).", ) parser.add_argument( "srcfile", metavar="src_dataset", type=str, help="The source dataset name. It can be either file name, " "URL of data source or subdataset name for multi-dataset files.", ) parser.add_argument( "dstfile", metavar="dst_dataset", type=str, help="The destination file name.", ) return parser def parse(self, argv) -> dict: def checkInt(s): try: int(s) return True except ValueError: return False def checkFloat(s): try: float(s) return True except ValueError: return False # We allowed "-b band_number_1 [... band_number_X] srcfile dstfile" syntax # before switching to ArgParse. But ArgParse doesn't like that. # So manually parse -b argument and strip it before ArgParse. count = len(argv) new_argv = [] i = 0 band_nums = [] src_nodata = [] dst_nodata = [] while i < count: arg = argv[i] if arg in ("-b", "-band", "--band"): if i == count - 1: raise Exception(f"Missing argument following {arg}: ") i += 1 if not checkInt(argv[i]): raise Exception(f"Argument following {arg} should be a integer") while i < count and checkInt(argv[i]): band_nums.append(int(argv[i])) i += 1 elif arg in ("-srcnodata", "-nodatavalue"): if i == count - 1: raise Exception(f"Missing argument following {arg}: ") i += 1 if not checkFloat(argv[i]): raise Exception(f"Argument following {arg} should be a float") while i < count and checkFloat(argv[i]): src_nodata.append(float(argv[i])) i += 1 elif arg in ("-dstnodata"): if i == count - 1: raise Exception(f"Missing argument following {arg}: ") i += 1 if not checkFloat(argv[i]): raise Exception(f"Argument following {arg} should be a float") while i < count and checkFloat(argv[i]): dst_nodata.append(float(argv[i])) i += 1 else: i += 1 new_argv.append(arg) kwargs = super(GDAL2XYZ, self).parse(new_argv) if band_nums: kwargs["band_nums"] = band_nums if src_nodata: kwargs["src_nodata"] = src_nodata if dst_nodata: kwargs["dst_nodata"] = dst_nodata return kwargs def augment_kwargs(self, kwargs) -> dict: if kwargs.get("allbands"): kwargs["band_nums"] = None elif not kwargs.get("band_nums"): kwargs["band_nums"] = 1 del kwargs["allbands"] return kwargs def doit(self, **kwargs): return gdal2xyz(**kwargs) def main(argv=sys.argv): return GDAL2XYZ().main(argv) if __name__ == "__main__": sys.exit(main(sys.argv))