24{
25 TCLAP::CmdLine cmd(
26 "Takes two DEMs located at the exact same spatial position (but at "
27 "different elevation) and calculates n raster DEMs located at "
28 "equidistant intervals between them (i.e. for n=1, one new raster "
29 "located precisely in the middle will be created).\n\n"
30 "OpenGeoSys-6 software, version " +
32 ".\n"
33 "Copyright (c) 2012-2025, OpenGeoSys Community "
34 "(http://www.opengeosys.org)",
36 TCLAP::ValueArg<std::size_t> number_arg(
37 "n", "number",
38 "number of rasters to be calculated, "
39 "(min = 0)",
40 false, 1, "int");
41 cmd.add(number_arg);
42 TCLAP::ValueArg<std::string> output_arg("o", "output-file",
43 "Output (.asc). Raster output file",
44 true, "", "OUTPUT_FILE");
45 cmd.add(output_arg);
46 TCLAP::ValueArg<std::string> input2_arg(
47 "", "file2", "Input (.asc). Second DEM-raster input file", true, "",
48 "INPUT_FILE");
49 cmd.add(input2_arg);
50 TCLAP::ValueArg<std::string> input1_arg(
51 "", "file1", "Input (.asc) First DEM-raster input file", true, "",
52 "INPUT_FILE");
53 cmd.add(input1_arg);
55 cmd.add(log_level_arg);
56
57 cmd.parse(argc, argv);
58
61
62 std::unique_ptr<GeoLib::Raster> dem1(
64 std::unique_ptr<GeoLib::Raster> dem2(
66
67 if (dem1 == nullptr || dem2 == nullptr)
68 {
69 return 1;
70 }
71
74
75 bool errors_found(false);
77 {
78 ERR(
"Origin x-coordinate is not the same in both raster files.\n");
79 errors_found = true;
80 }
82 {
83 ERR(
"Origin y-coordinate is not the same in both raster files.\n");
84 errors_found = true;
85 }
87 {
88 ERR(
"Cellsize is not the same in both raster files.\n");
89 errors_found = true;
90 }
92 {
93 ERR(
"Raster width is not the same in both raster files.\n");
94 errors_found = true;
95 }
97 {
98 ERR(
"Raster height is not the same in both raster files.\n");
99 errors_found = true;
100 }
101
102 if (errors_found)
103 {
104 return 2;
105 }
106
107 std::size_t const n = number_arg.getValue();
108 std::vector<std::vector<double>> raster;
109 for (std::size_t i = 0; i < n; ++i)
110 {
111 std::vector<double> r;
113 raster.push_back(r);
114 }
115
116 auto it2 = dem2->begin();
117 for (auto it1 = dem1->begin(); it1 != dem1->end(); ++it1)
118 {
119 if (it2 == dem2->end())
120 {
121 ERR(
"Error: File 2 is shorter than File 1.");
122 return 3;
123 }
125 {
126 for (std::size_t i = 0; i < n; ++i)
127 {
128 raster[i].push_back(h1.
no_data);
129 }
130 }
131 else
132 {
133 double const min = std::min(*it1, *it2);
134 double const max = std::max(*it1, *it2);
135 double const step = (max - min) / static_cast<double>(n + 1);
136 for (std::size_t i = 0; i < n; ++i)
137 {
138 raster[i].push_back(max - ((i + 1) * step));
139 }
140 }
141 it2++;
142 }
143 if (it2 != dem2->end())
144 {
145 ERR(
"Error: File 1 is shorter than File 2.");
146 return 3;
147 }
148
149 std::string const filename = output_arg.getValue();
150 for (std::size_t i = 0; i < n; ++i)
151 {
154
157 r, basename + std::to_string(i) + ext);
158 INFO(
"Layer {:d} written.", i + 1);
159 }
160 return EXIT_SUCCESS;
161}
void INFO(fmt::format_string< Args... > fmt, Args &&... args)
void ERR(fmt::format_string< Args... > fmt, Args &&... args)
static void writeRasterAsASC(GeoLib::Raster const &raster, std::string const &file_name)
Writes an Esri asc-file.
static GeoLib::Raster * readRaster(std::string const &fname)
Class Raster is used for managing raster data.
TCLAP::ValueArg< std::string > makeLogLevelArg()
void initOGSLogger(std::string const &log_level)
std::string getFileExtension(const std::string &path)
std::string dropFileExtension(std::string const &filename)
GITINFOLIB_EXPORT const std::string ogs_version