| 140 | } |
| 141 | |
| 142 | void ReprojectionFilter::createTransform(const SpatialReference& srsSRS) |
| 143 | { |
| 144 | if (m_inferInputSRS) |
| 145 | { |
| 146 | m_inSRS = srsSRS; |
| 147 | if (m_inSRS.empty()) |
| 148 | throwError("source data has no spatial reference and " |
| 149 | "none is specified with the 'in_srs' option."); |
| 150 | } |
| 151 | |
| 152 | |
| 153 | // If either vector is empty, GDAL's default ordering is used. |
| 154 | if (m_inAxisOrdering.size() || m_outAxisOrdering.size()) |
| 155 | { |
| 156 | |
| 157 | m_transform.reset(new SrsTransform(m_inSRS, |
| 158 | m_inAxisOrdering, |
| 159 | m_outSRS, |
| 160 | m_outAxisOrdering)); |
| 161 | } else { |
| 162 | m_transform.reset(new SrsTransform(m_inSRS, m_outSRS)); |
| 163 | } |
| 164 | |
| 165 | |
| 166 | if (!pdal::Utils::compare_approx(m_inCoordEpochArg, 0.0f, 0.00f)) |
| 167 | { |
| 168 | m_transform->setSrcEpoch(m_inCoordEpochArg); |
| 169 | } |
| 170 | |
| 171 | if (!pdal::Utils::compare_approx(m_outCoordEpochArg, 0.0f, 0.00f)) |
| 172 | { |
| 173 | m_transform->setDstEpoch(m_outCoordEpochArg); |
| 174 | } |
| 175 | |
| 176 | } |
| 177 | |
| 178 | |
| 179 | PointViewSet ReprojectionFilter::run(PointViewPtr view) |
nothing calls this directly
no test coverage detected