FD2ADDTOPO - Revision history

Ak3gseis at 12:53, 10 October 2022

2022-10-10T12:53:35Z

← Older revision		Revision as of 12:53, 10 October 2022
Line 55:		Line 55:
	\|language=Fortran2003		\|language=Fortran2003
	\|add_software= -		\|add_software= -
	\|contact_original=~~[mailto:marcus.boehlich@baw.de~~ M. Boehlich], ~~[mailto:guntram.seiss@baw.de~~ G. Seiß]		\|contact_original=M. Boehlich, G. Seiß
	\|contact_maintenance=[mailto:~~guntram~~.~~seiß~~@baw.de ~~Group~~ PRE]		\|contact_maintenance=[mailto:pre.proghome@baw.de working group PRE]
	\|documentation=$PROGHOME/examples/fd2addtopo/*		\|documentation=$PROGHOME/examples/fd2addtopo/*
	}}		}}

Ak3gseis: Change contact, version, supported EPSG

2022-09-06T09:22:23Z

Change contact, version, supported EPSG

← Older revision		Revision as of 09:22, 6 September 2022
Line 2:		Line 2:
	\|name_de=FD2ADDTOPO		\|name_de=FD2ADDTOPO
	\|name=FD2ADDTOPO		\|name=FD2ADDTOPO
	\|version=~~June 2019~~		\|version=October 2021
	\|version_descr=June 2019		\|version_descr=June 2019
	\|catchwords=		\|catchwords=
Line 45:		Line 45:


	~~Supported~~ EPSG codes ~~are:~~		For supported EPSG codes see [[GEOTRANSFORMER]]!
	* Gauß-Krüger 3 Grad / DHDN: EPSG 31466-31469
	* Gauß-Krüger 3 Grad / STN : EPSG 02398-02399
	* UTM / ETRS89 : EPSG 25831-25833
	* UTM / ED50 : EPSG 23031-23033
	* UTM / WGS84 (North) : EPSG 32601-32660
	* UTM / WGS84 (North) : EPSG 32701-32760
	* UPS / WGS84 : EPSG 32661, 32761
	* RD Amersfoort New : EPSG 28992
	* Geographic ETRS89 : EPSG 4258
	* Geographic WGS84 : EPSG 4326
	* Geographic ED 50 : EPSG 4230
	For the coordinate shift from DHDN/STN to ETRS89 the NTv2 method is used with the grid shift files, which are		For the coordinate shift from DHDN/STN to ETRS89 the NTv2 method is used with the grid shift files, which are
	defined by NTV2GRID (DHDN) and NTV2STN (STN).		defined by NTV2GRID (DHDN) and NTV2STN (STN).
Line 64:		Line 53:
	\|postprocessor=		\|postprocessor=
	[[FD2DEL]], [[FD2TRIM]], [[FDGITTER05]], GISMO		[[FD2DEL]], [[FD2TRIM]], [[FDGITTER05]], GISMO
	\|language=~~Fortran95~~		\|language=Fortran2003
	\|add_software= -		\|add_software= -
	\|contact_original=[mailto:marcus.boehlich@baw.de M. Boehlich], [mailto:guntram.~~seiß~~@baw.de G. Seiß]		\|contact_original=[mailto:marcus.boehlich@baw.de M. Boehlich], [mailto:guntram.seiss@baw.de G. Seiß]
	\|contact_maintenance=[mailto:guntram.seiß@baw.de ~~G. Seiß~~]		\|contact_maintenance=[mailto:guntram.seiß@baw.de Group PRE]
	\|documentation=$PROGHOME/examples/fd2addtopo/*		\|documentation=$PROGHOME/examples/fd2addtopo/*
	}}		}}

imported>Seiss Guntram: Denglish to English

2019-06-18T13:29:54Z

Denglish to English

← Older revision		Revision as of 13:29, 18 June 2019
Line 54:		Line 54:
	* UPS / WGS84 : EPSG 32661, 32761		* UPS / WGS84 : EPSG 32661, 32761
	* RD Amersfoort New : EPSG 28992		* RD Amersfoort New : EPSG 28992
	* ~~Geographisch~~ ETRS89 : EPSG 4258		* Geographic ETRS89 : EPSG 4258
	* ~~Geographisch~~ WGS84 : EPSG 4326		* Geographic WGS84 : EPSG 4326
	* ~~Geographisch~~ ED 50 : EPSG 4230		* Geographic ED 50 : EPSG 4230
	For the coordinate shift from DHDN/STN to ETRS89 the NTv2 method is used with the grid shift files, which are		For the coordinate shift from DHDN/STN to ETRS89 the NTv2 method is used with the grid shift files, which are
	defined by NTV2GRID (DHDN) and NTV2STN (STN).		defined by NTV2GRID (DHDN) and NTV2STN (STN).

imported>Seiss Guntram: decompression of input files on demand

2019-06-18T13:23:48Z

decompression of input files on demand

← Older revision		Revision as of 13:23, 18 June 2019
Line 2:		Line 2:
	\|name_de=FD2ADDTOPO		\|name_de=FD2ADDTOPO
	\|name=FD2ADDTOPO		\|name=FD2ADDTOPO
	\|version=~~December 2013~~		\|version=June 2019
	\|version_descr=~~December 2013~~		\|version_descr=June 2019
	\|catchwords=		\|catchwords=
	finite differences grid (type FIDISOR/FIDIRB)<br />		finite differences grid (type FIDISOR/FIDIRB)<br />
Line 19:		Line 19:
	# input control file [[FD2ADDTOPO.DAT\|fd2addtopo.dat]]<br />It is possible to generate input control files for the procedure GITGENER automatically. The program gen_fd2addtopo_single.hp generates files that control the procedure for single echo sounding files. The program gen_fd2addtopo_cover.hp generates files that control creating a topography out of echo sounding files that cover the area of a reference topography.		# input control file [[FD2ADDTOPO.DAT\|fd2addtopo.dat]]<br />It is possible to generate input control files for the procedure GITGENER automatically. The program gen_fd2addtopo_single.hp generates files that control the procedure for single echo sounding files. The program gen_fd2addtopo_cover.hp generates files that control creating a topography out of echo sounding files that cover the area of a reference topography.
	# (when GITGENER shall be executed or a topography shall be updated with ADDTOP) (example) topography (type [[TOPO.BIN\|topo.bin]]).		# (when GITGENER shall be executed or a topography shall be updated with ADDTOP) (example) topography (type [[TOPO.BIN\|topo.bin]]).
	# (when GITGENER shall be executed) ~~file~~ which contains echo soundings (type [[GEOM.DAT\|geom.dat]] or FORTRAN format '(2(f10.2,1x),f6.2)'. Other formats with ~~seperating~~ blanks are supported too.		# (when GITGENER shall be executed) files which contains echo soundings (type [[GEOM.DAT\|geom.dat]] or FORTRAN format '(2(f10.2,1x),f6.2)'. Other formats with separating blanks are supported too. Packed files are allowed (GZIP, COMPRESS, ZIP)!
	# (when GITPOL or ADDTOP shall run independently from GITGENER) difference topography (type [[TOPO.BIN\|topo.bin]]).		# (when GITPOL or ADDTOP shall run independently from GITGENER) difference topography (type [[TOPO.BIN\|topo.bin]]).
	# (when GITPOL or ADDTOP shall run independently from GITGENER) compatible mask file for difference topography (type topo.bin.log).		# (when GITPOL or ADDTOP shall run independently from GITGENER) compatible mask file for difference topography (type topo.bin.log).
	# (optional) digitized structures (file type [[DIGI.GKK\|digi.gkk]])		# (optional) digitized structures (file type [[DIGI.GKK\|digi.gkk]]). Packed files are allowed (GZIP, COMPRESS, ZIP)!
	\|outputfiles=		\|outputfiles=
	# after GITGENER) difference topography (type [[TOPO.BIN\|topo.bin]])		# after GITGENER) difference topography (type [[TOPO.BIN\|topo.bin]])
Line 42:		Line 42:

	Remark:		Remark:
	Coordinates of raw data files are now transformed to a user defined coordinate reference systen (CRS) if necessary. To support this feature, the user must include a magic comment with the content "C CRS=xxxxx" at the beginning of the file. Here "xxxxx" stands for the EPSG code of the CRS. Additionally the user has to set the environment variable BAWCRS to a valid CRS. All Input files of type topo.bin must be realized in this user CRS.		Coordinates of raw data files are now transformed to a user defined coordinate reference systen (CRS) if necessary. To support this feature, the user '''must''' include a magic comment with the content "C CRS=xxxxx" at the beginning of the file. Here "xxxxx" stands for the EPSG code of the CRS. Additionally the user '''has to''' set the environment variable BAWCRS to a valid CRS. All Input files of type topo.bin must be realized in this user CRS.


	Supported EPSG codes are:		Supported EPSG codes are:
	* ~~Gauss~~-~~Krueger~~ 3 ~~degree~~ / DHDN: EPSG 31466-31469		* Gauß-Krüger 3 Grad / DHDN: EPSG 31466-31469
	* ~~Gauss~~-~~Krueger~~ 3 ~~degree~~ / STN : EPSG 02398-02399		* Gauß-Krüger 3 Grad / STN : EPSG 02398-02399
	* UTM / ETRS89 : EPSG 25831-25833		* UTM / ETRS89 : EPSG 25831-25833
	* UTM / ED50 : EPSG 23031-23033		* UTM / ED50 : EPSG 23031-23033
	* RD Amersfoort New : EPSG 28992		* UTM / WGS84 (North) : EPSG 32601-32660
	* ~~Geographic~~ ETRS89 : EPSG 4258		* UTM / WGS84 (North) : EPSG 32701-32760
	* ~~Geographic~~ WGS84 : EPSG 4326		* UPS / WGS84 : EPSG 32661, 32761
	* ~~Geographic~~ ED 50 : EPSG 4230		* RD Amersfoort New : EPSG 28992
			* Geographisch ETRS89 : EPSG 4258
			* Geographisch WGS84 : EPSG 4326
			* Geographisch ED 50 : EPSG 4230
	For the coordinate shift from DHDN/STN to ETRS89 the NTv2 method is used with the grid shift files, which are		For the coordinate shift from DHDN/STN to ETRS89 the NTv2 method is used with the grid shift files, which are
	defined by NTV2GRID (DHDN) and NTV2STN (STN).		defined by NTV2GRID (DHDN) and NTV2STN (STN).

imported>Seiss Guntram: Some additions/ corrections to the last edit

2013-12-18T16:06:16Z

Some additions/ corrections to the last edit

← Older revision		Revision as of 16:06, 18 December 2013
Line 42:		Line 42:

	Remark:		Remark:
	Coordinates of raw data files are now transformed to a user defined coordinate reference systen (CRS) if necessary. To support this feature, the user must include a magic comment with the content "C CRS=xxxxx" at the beginning of the file. Additionally the user has to set the ~~Environment~~ variable BAWCRS to a valid CRS.		Coordinates of raw data files are now transformed to a user defined coordinate reference systen (CRS) if necessary. To support this feature, the user must include a magic comment with the content "C CRS=xxxxx" at the beginning of the file. Here "xxxxx" stands for the EPSG code of the CRS. Additionally the user has to set the environment variable BAWCRS to a valid CRS. All Input files of type topo.bin must be realized in this user CRS.
	~~Here "xxxxx" stands for the EPSG code~~ of ~~the~~ CRS. Supported EPSG codes are:

			Supported EPSG codes are:
	* Gauss-Krueger 3 degree / DHDN: EPSG 31466-31469		* Gauss-Krueger 3 degree / DHDN: EPSG 31466-31469
	* Gauss-Krueger 3 degree / STN : EPSG 02398-02399		* Gauss-Krueger 3 degree / STN : EPSG 02398-02399
Line 58:		Line 60:
	Digitizer (PC), [[FD2BASIS]], [[FDGITTER05]]		Digitizer (PC), [[FD2BASIS]], [[FDGITTER05]]
	\|postprocessor=		\|postprocessor=
	[[FD2DEL]], [[FD2TRIM]], [[FDGITTER05]], [[GISMO]]		[[FD2DEL]], [[FD2TRIM]], [[FDGITTER05]], GISMO
	\|language=~~Fortran90~~		\|language=Fortran95
	\|add_software= -		\|add_software= -
	\|contact_original=[mailto:marcus.boehlich@baw.de M. Boehlich], [mailto:guntram.seiß@baw.de G. Seiß]		\|contact_original=[mailto:marcus.boehlich@baw.de M. Boehlich], [mailto:guntram.seiß@baw.de G. Seiß]

imported>Seiss Guntram: Support of raw data with different coordinate reference systems

2013-12-18T15:58:19Z

Support of raw data with different coordinate reference systems

← Older revision		Revision as of 15:58, 18 December 2013
Line 2:		Line 2:
	\|name_de=FD2ADDTOPO		\|name_de=FD2ADDTOPO
	\|name=FD2ADDTOPO		\|name=FD2ADDTOPO
	\|version=~~Juli 2010~~		\|version=December 2013
	\|version_descr=~~Juli 2010~~		\|version_descr=December 2013
	\|catchwords=		\|catchwords=
	finite differences grid (type FIDISOR/FIDIRB)<br />		finite differences grid (type FIDISOR/FIDIRB)<br />
Line 40:		Line 40:
	During the second step ( GITPOL ) depth values for further cells will be interpolated from cells, which contain already a real depth. This step eventually has to run twice for situations with a very scattered distribution of the echo soundings. The step can be suppressed by setting the count of interpolations to 0 inside the input control file. By this interpolation cells are searched inside a region of neighbourhood, which contain valid depths. The search is done radial in eight directions (north, northeast, east, southeast, south, southwest, west, northwest). Then a inverse distance weighted average of the found depth values is calculated (simple Shepard interpolation).		During the second step ( GITPOL ) depth values for further cells will be interpolated from cells, which contain already a real depth. This step eventually has to run twice for situations with a very scattered distribution of the echo soundings. The step can be suppressed by setting the count of interpolations to 0 inside the input control file. By this interpolation cells are searched inside a region of neighbourhood, which contain valid depths. The search is done radial in eight directions (north, northeast, east, southeast, south, southwest, west, northwest). Then a inverse distance weighted average of the found depth values is calculated (simple Shepard interpolation).
	During the third step ( ADDTOP ) depth values of the (example) topography are updated where they exist.		During the third step ( ADDTOP ) depth values of the (example) topography are updated where they exist.

			Remark:
			Coordinates of raw data files are now transformed to a user defined coordinate reference systen (CRS) if necessary. To support this feature, the user must include a magic comment with the content "C CRS=xxxxx" at the beginning of the file. Additionally the user has to set the Environment variable BAWCRS to a valid CRS.
			Here "xxxxx" stands for the EPSG code of the CRS. Supported EPSG codes are:
			* Gauss-Krueger 3 degree / DHDN: EPSG 31466-31469
			* Gauss-Krueger 3 degree / STN : EPSG 02398-02399
			* UTM / ETRS89 : EPSG 25831-25833
			* UTM / ED50 : EPSG 23031-23033
			* RD Amersfoort New : EPSG 28992
			* Geographic ETRS89 : EPSG 4258
			* Geographic WGS84 : EPSG 4326
			* Geographic ED 50 : EPSG 4230
			For the coordinate shift from DHDN/STN to ETRS89 the NTv2 method is used with the grid shift files, which are
			defined by NTV2GRID (DHDN) and NTV2STN (STN).

	\|preprocessor=		\|preprocessor=
	Digitizer (PC), [[FD2BASIS]], [[FDGITTER05]]		Digitizer (PC), [[FD2BASIS]], [[FDGITTER05]]
	\|postprocessor=		\|postprocessor=
	[[FD2DEL]], [[FD2TRIM]], [[FDGITTER05]]		[[FD2DEL]], [[FD2TRIM]], [[FDGITTER05]], [[GISMO]]
	\|language=Fortran90		\|language=Fortran90
	\|add_software= -		\|add_software= -

imported>Seiss Guntram: More explanation of the method GITPOL

2010-07-12T09:35:19Z

More explanation of the method GITPOL

← Older revision		Revision as of 09:35, 12 July 2010
Line 2:		Line 2:
	\|name_de=FD2ADDTOPO		\|name_de=FD2ADDTOPO
	\|name=FD2ADDTOPO		\|name=FD2ADDTOPO
	\|version=~~December 1998~~		\|version=Juli 2010
	\|version_descr=~~December 1998~~		\|version_descr=Juli 2010
	\|catchwords=		\|catchwords=
	finite differences grid (type FIDISOR/FIDIRB)<br />		finite differences grid (type FIDISOR/FIDIRB)<br />
Line 38:		Line 38:
	* 1: the weightning of data decreases linear with increasing distance to the center of the cell. data points outside the correlation radius have no influence.		* 1: the weightning of data decreases linear with increasing distance to the center of the cell. data points outside the correlation radius have no influence.
	* the weightning of data decreases exponentially with increasing distance to the center of the cell.		* the weightning of data decreases exponentially with increasing distance to the center of the cell.
	During the second step ( GITPOL ) depth values for further cells will be interpolated from cells, which contain already a real depth. This step eventually has to run twice for situations with a very scattered distribution of the echo soundings. The step can be suppressed by setting the count of interpolations to 0 inside the input control file.		During the second step ( GITPOL ) depth values for further cells will be interpolated from cells, which contain already a real depth. This step eventually has to run twice for situations with a very scattered distribution of the echo soundings. The step can be suppressed by setting the count of interpolations to 0 inside the input control file. By this interpolation cells are searched inside a region of neighbourhood, which contain valid depths. The search is done radial in eight directions (north, northeast, east, southeast, south, southwest, west, northwest). Then a inverse distance weighted average of the found depth values is calculated (simple Shepard interpolation).
	During the third step ( ADDTOP ) depth values of the (example) topography are updated where they exist.		During the third step ( ADDTOP ) depth values of the (example) topography are updated where they exist.
	\|preprocessor=		\|preprocessor=

imported>BAWiki 1: Created page with '{{ProgramDescription |name_de=FD2ADDTOPO |name=FD2ADDTOPO |version=December 1998 |version_descr=December 1998 |catchwords= finite differences grid (type FIDISOR/FIDIRB)
…'

2010-05-10T09:22:55Z

Created page with '{{ProgramDescription |name_de=FD2ADDTOPO |name=FD2ADDTOPO |version=December 1998 |version_descr=December 1998 |catchwords= finite differences grid (type FIDISOR/FIDIRB) …'

New page

{{ProgramDescription
|name_de=FD2ADDTOPO
|name=FD2ADDTOPO
|version=December 1998
|version_descr=December 1998
|catchwords=
finite differences grid (type FIDISOR/FIDIRB) 
update of topography 
echo soundings 
|shortdescription=
The program FD2ADDTOPO creates a difference topography from echo soundings in the format of the given example topography (type FIDISOR/FIDIRB). This difference topography (after optional interpolation) is used to update the grid cells of the original topography with new values through direct replacement. The last difference topography can be written into new ASCII files.
The program contains three components:
# procedure GITGENER creates a new difference topography from echo soundings.
# procedure GITPOL interpolates an already existing or just created difference topography.
# procedure ADDTOP updates an old topography with values of an already existing or just created difference topography.
# procedure STORERAS stores the calculated grid cell center points into new ASCII-files (x,y,z). Two modi are possible: the first stores all data into a single file, the second creates a file for each window of one square kilometer size..
With the help of automatically generated input control and batch files it is possible to estimate the quality of echo soundings. For further information please contact the [mailto:peter.schade@baw.de developer] of this method.
|inputfiles=
# input control file [[FD2ADDTOPO.DAT|fd2addtopo.dat]] It is possible to generate input control files for the procedure GITGENER automatically. The program gen_fd2addtopo_single.hp generates files that control the procedure for single echo sounding files. The program gen_fd2addtopo_cover.hp generates files that control creating a topography out of echo sounding files that cover the area of a reference topography.
# (when GITGENER shall be executed or a topography shall be updated with ADDTOP) (example) topography (type [[TOPO.BIN|topo.bin]]).
# (when GITGENER shall be executed) file which contains echo soundings (type [[GEOM.DAT|geom.dat]] or FORTRAN format '(2(f10.2,1x),f6.2)'. Other formats with seperating blanks are supported too.
# (when GITPOL or ADDTOP shall run independently from GITGENER) difference topography (type [[TOPO.BIN|topo.bin]]).
# (when GITPOL or ADDTOP shall run independently from GITGENER) compatible mask file for difference topography (type topo.bin.log).
# (optional) digitized structures (file type [[DIGI.GKK|digi.gkk]])
|outputfiles=
# after GITGENER) difference topography (type [[TOPO.BIN|topo.bin]])
# (after GITGENER) compatible mask file for difference topography (type topo.bin.log).
# (after GITPOL) interpolated difference topography (type [[TOPO.BIN|topo.bin]])
# (after GITPOL) compatible mask file for interpolated difference topography (type topo.bin.log).
# (after ADDTOP) updated topography (type [[TOPO.BIN|topo.bin]]).
# (after STORERAS) one ore several files of type [[GEOM.DAT|geom.dat]] or any other (xyz) format.
# printer output file fd2addtopo.sdr
|methodology=
In a first step ( GITGENER ) the program FD2ADDTOPO creates a FIDISOR/FIDIRB topography file which contains actual depth values. The user can define how many echo soundings have to exist at least in a grid cell in order to assign a depth value. To cells with less echo soundings a default value is assigned. Furthermore a histogram with the number of echo soundings per cell is evaluated. The printer output file contains the histogram.
Additionally a mask file will be created which contains the information whether a depth value is assigned to a grid cell or not. All cells which are not updated contain the value -99.9 as a dummy depth. The filenames of echo sounding data have to be specified with or without a preceeding sign. The meaning of the sign is as follows: If the sign is "-", the depths are multiplied by -1, otherwise the depths are used as they are given.
An avaraging is done for all data which are located in the same grid cell. The following averaging methods can be used:
* 0: alle data are weighted equally with 1.0
* 1: the weightning of data decreases linear with increasing distance to the center of the cell. data points outside the correlation radius have no influence.
* the weightning of data decreases exponentially with increasing distance to the center of the cell.
During the second step ( GITPOL ) depth values for further cells will be interpolated from cells, which contain already a real depth. This step eventually has to run twice for situations with a very scattered distribution of the echo soundings. The step can be suppressed by setting the count of interpolations to 0 inside the input control file.
During the third step ( ADDTOP ) depth values of the (example) topography are updated where they exist.
|preprocessor=
Digitizer (PC), [[FD2BASIS]], [[FDGITTER05]]
|postprocessor=
[[FD2DEL]], [[FD2TRIM]], [[FDGITTER05]]
|language=Fortran90
|add_software= -
|contact_original=[mailto:marcus.boehlich@baw.de M. Boehlich], [mailto:guntram.seiß@baw.de G. Seiß]
|contact_maintenance=[mailto:guntram.seiß@baw.de G. Seiß]
|documentation=$PROGHOME/examples/fd2addtopo/*
}}