Mercurial > hg > orthanc-stone
annotate Framework/Loaders/DicomStructureSetLoader.cpp @ 841:266e2b0b9abc
better error reporting in DicomStructureSetLoader + fixed POST request logic
in WebAssemblyOracle + support for LookupTableTextureSceneLayer in OpenGL (NOT
using shaders!) (2 new files) + a few small non-functional changes
| author | Benjamin Golinvaux <bgo@osimis.io> |
|---|---|
| date | Tue, 11 Jun 2019 15:41:21 +0200 |
| parents | 28f99af358fa |
| children | 2b245953b44b |
| rev | line source |
|---|---|
| 815 | 1 /** |
| 2 * Stone of Orthanc | |
| 3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics | |
| 4 * Department, University Hospital of Liege, Belgium | |
| 5 * Copyright (C) 2017-2019 Osimis S.A., Belgium | |
| 6 * | |
| 7 * This program is free software: you can redistribute it and/or | |
| 8 * modify it under the terms of the GNU Affero General Public License | |
| 9 * as published by the Free Software Foundation, either version 3 of | |
| 10 * the License, or (at your option) any later version. | |
| 11 * | |
| 12 * This program is distributed in the hope that it will be useful, but | |
| 13 * WITHOUT ANY WARRANTY; without even the implied warranty of | |
| 14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
| 15 * Affero General Public License for more details. | |
| 16 * | |
| 17 * You should have received a copy of the GNU Affero General Public License | |
| 18 * along with this program. If not, see <http://www.gnu.org/licenses/>. | |
| 19 **/ | |
| 20 | |
| 21 | |
| 22 #include "DicomStructureSetLoader.h" | |
| 23 | |
| 24 #include "../Scene2D/PolylineSceneLayer.h" | |
| 25 #include "../Toolbox/GeometryToolbox.h" | |
| 26 | |
| 27 namespace OrthancStone | |
| 28 { | |
| 29 class DicomStructureSetLoader::AddReferencedInstance : public LoaderStateMachine::State | |
| 30 { | |
| 31 private: | |
| 32 std::string instanceId_; | |
| 33 | |
| 34 public: | |
| 35 AddReferencedInstance(DicomStructureSetLoader& that, | |
| 36 const std::string& instanceId) : | |
| 37 State(that), | |
| 38 instanceId_(instanceId) | |
| 39 { | |
| 40 } | |
| 41 | |
| 42 virtual void Handle(const OrthancRestApiCommand::SuccessMessage& message) | |
| 43 { | |
| 44 Json::Value tags; | |
| 45 message.ParseJsonBody(tags); | |
| 46 | |
| 47 Orthanc::DicomMap dicom; | |
| 48 dicom.FromDicomAsJson(tags); | |
| 49 | |
| 50 DicomStructureSetLoader& loader = GetLoader<DicomStructureSetLoader>(); | |
| 51 loader.content_->AddReferencedSlice(dicom); | |
| 52 | |
| 53 loader.countProcessedInstances_ ++; | |
| 54 assert(loader.countProcessedInstances_ <= loader.countReferencedInstances_); | |
| 55 | |
| 56 if (loader.countProcessedInstances_ == loader.countReferencedInstances_) | |
| 57 { | |
| 58 // All the referenced instances have been loaded, finalize the RT-STRUCT | |
| 59 loader.content_->CheckReferencedSlices(); | |
| 60 loader.revision_++; | |
| 61 } | |
| 62 } | |
| 63 }; | |
| 64 | |
| 65 | |
| 66 // State that converts a "SOP Instance UID" to an Orthanc identifier | |
| 67 class DicomStructureSetLoader::LookupInstance : public LoaderStateMachine::State | |
| 68 { | |
| 69 private: | |
| 70 std::string sopInstanceUid_; | |
| 71 | |
| 72 public: | |
| 73 LookupInstance(DicomStructureSetLoader& that, | |
| 74 const std::string& sopInstanceUid) : | |
| 75 State(that), | |
| 76 sopInstanceUid_(sopInstanceUid) | |
| 77 { | |
| 78 } | |
| 79 | |
| 80 virtual void Handle(const OrthancRestApiCommand::SuccessMessage& message) | |
| 81 { | |
| 82 DicomStructureSetLoader& loader = GetLoader<DicomStructureSetLoader>(); | |
| 83 | |
| 84 Json::Value lookup; | |
| 85 message.ParseJsonBody(lookup); | |
| 86 | |
| 87 if (lookup.type() != Json::arrayValue || | |
| 88 lookup.size() != 1 || | |
| 89 !lookup[0].isMember("Type") || | |
| 90 !lookup[0].isMember("Path") || | |
| 91 lookup[0]["Type"].type() != Json::stringValue || | |
| 92 lookup[0]["ID"].type() != Json::stringValue || | |
| 93 lookup[0]["Type"].asString() != "Instance") | |
| 94 { | |
|
841
266e2b0b9abc
better error reporting in DicomStructureSetLoader + fixed POST request logic
Benjamin Golinvaux <bgo@osimis.io>
parents:
828
diff
changeset
|
95 std::stringstream msg; |
|
266e2b0b9abc
better error reporting in DicomStructureSetLoader + fixed POST request logic
Benjamin Golinvaux <bgo@osimis.io>
parents:
828
diff
changeset
|
96 msg << "Unknown resource! message.GetAnswer() = " << message.GetAnswer() << " message.GetAnswerHeaders() = "; |
|
266e2b0b9abc
better error reporting in DicomStructureSetLoader + fixed POST request logic
Benjamin Golinvaux <bgo@osimis.io>
parents:
828
diff
changeset
|
97 for (const auto& it : message.GetAnswerHeaders()) |
|
266e2b0b9abc
better error reporting in DicomStructureSetLoader + fixed POST request logic
Benjamin Golinvaux <bgo@osimis.io>
parents:
828
diff
changeset
|
98 { |
|
266e2b0b9abc
better error reporting in DicomStructureSetLoader + fixed POST request logic
Benjamin Golinvaux <bgo@osimis.io>
parents:
828
diff
changeset
|
99 msg << "\nkey: \"" << it.first << "\" value: \"" << it.second << "\"\n"; |
|
266e2b0b9abc
better error reporting in DicomStructureSetLoader + fixed POST request logic
Benjamin Golinvaux <bgo@osimis.io>
parents:
828
diff
changeset
|
100 } |
|
266e2b0b9abc
better error reporting in DicomStructureSetLoader + fixed POST request logic
Benjamin Golinvaux <bgo@osimis.io>
parents:
828
diff
changeset
|
101 auto msgStr = msg.str(); |
|
266e2b0b9abc
better error reporting in DicomStructureSetLoader + fixed POST request logic
Benjamin Golinvaux <bgo@osimis.io>
parents:
828
diff
changeset
|
102 LOG(ERROR) << msgStr; |
| 815 | 103 throw Orthanc::OrthancException(Orthanc::ErrorCode_UnknownResource); |
| 104 } | |
| 105 | |
| 106 const std::string instanceId = lookup[0]["ID"].asString(); | |
| 107 | |
| 108 { | |
| 109 std::auto_ptr<OrthancRestApiCommand> command(new OrthancRestApiCommand); | |
| 110 command->SetHttpHeader("Accept-Encoding", "gzip"); | |
| 111 command->SetUri("/instances/" + instanceId + "/tags"); | |
| 112 command->SetPayload(new AddReferencedInstance(loader, instanceId)); | |
| 113 Schedule(command.release()); | |
| 114 } | |
| 115 } | |
| 116 }; | |
| 117 | |
| 118 | |
| 119 class DicomStructureSetLoader::LoadStructure : public LoaderStateMachine::State | |
| 120 { | |
| 121 public: | |
| 122 LoadStructure(DicomStructureSetLoader& that) : | |
| 123 State(that) | |
| 124 { | |
| 125 } | |
| 126 | |
| 127 virtual void Handle(const OrthancRestApiCommand::SuccessMessage& message) | |
| 128 { | |
| 129 DicomStructureSetLoader& loader = GetLoader<DicomStructureSetLoader>(); | |
| 130 | |
| 131 { | |
| 132 OrthancPlugins::FullOrthancDataset dicom(message.GetAnswer()); | |
| 133 loader.content_.reset(new DicomStructureSet(dicom)); | |
| 134 } | |
| 135 | |
| 136 std::set<std::string> instances; | |
| 137 loader.content_->GetReferencedInstances(instances); | |
| 138 | |
| 828 | 139 loader.countReferencedInstances_ = static_cast<unsigned int>(instances.size()); |
| 815 | 140 |
| 141 for (std::set<std::string>::const_iterator | |
| 142 it = instances.begin(); it != instances.end(); ++it) | |
| 143 { | |
| 144 std::auto_ptr<OrthancRestApiCommand> command(new OrthancRestApiCommand); | |
| 145 command->SetUri("/tools/lookup"); | |
| 146 command->SetMethod(Orthanc::HttpMethod_Post); | |
| 147 command->SetBody(*it); | |
| 148 command->SetPayload(new LookupInstance(loader, *it)); | |
|
841
266e2b0b9abc
better error reporting in DicomStructureSetLoader + fixed POST request logic
Benjamin Golinvaux <bgo@osimis.io>
parents:
828
diff
changeset
|
149 //LOG(TRACE) << "About to schedule a /tools/lookup POST request. URI = " << command->GetUri() << " Body size = " << (*it).size() << " Body = " << (*it) << "\n"; |
| 815 | 150 Schedule(command.release()); |
| 151 } | |
| 152 } | |
| 153 }; | |
| 154 | |
| 155 | |
| 156 class DicomStructureSetLoader::Slice : public IExtractedSlice | |
| 157 { | |
| 158 private: | |
| 159 const DicomStructureSet& content_; | |
| 160 uint64_t revision_; | |
| 161 bool isValid_; | |
| 162 | |
| 163 public: | |
| 164 Slice(const DicomStructureSet& content, | |
| 165 uint64_t revision, | |
| 166 const CoordinateSystem3D& cuttingPlane) : | |
| 167 content_(content), | |
| 168 revision_(revision) | |
| 169 { | |
| 170 bool opposite; | |
| 171 | |
| 172 const Vector normal = content.GetNormal(); | |
| 173 isValid_ = ( | |
| 174 GeometryToolbox::IsParallelOrOpposite(opposite, normal, cuttingPlane.GetNormal()) || | |
| 175 GeometryToolbox::IsParallelOrOpposite(opposite, normal, cuttingPlane.GetAxisX()) || | |
| 176 GeometryToolbox::IsParallelOrOpposite(opposite, normal, cuttingPlane.GetAxisY())); | |
| 177 } | |
| 178 | |
| 179 virtual bool IsValid() | |
| 180 { | |
| 181 return isValid_; | |
| 182 } | |
| 183 | |
| 184 virtual uint64_t GetRevision() | |
| 185 { | |
| 186 return revision_; | |
| 187 } | |
| 188 | |
| 189 virtual ISceneLayer* CreateSceneLayer(const ILayerStyleConfigurator* configurator, | |
| 190 const CoordinateSystem3D& cuttingPlane) | |
| 191 { | |
| 192 assert(isValid_); | |
| 193 | |
| 194 std::auto_ptr<PolylineSceneLayer> layer(new PolylineSceneLayer); | |
| 195 layer->SetThickness(2); | |
| 196 | |
| 197 for (size_t i = 0; i < content_.GetStructuresCount(); i++) | |
| 198 { | |
| 199 const Color& color = content_.GetStructureColor(i); | |
| 200 | |
| 201 std::vector< std::vector<DicomStructureSet::PolygonPoint> > polygons; | |
| 202 | |
| 203 if (content_.ProjectStructure(polygons, i, cuttingPlane)) | |
| 204 { | |
| 205 for (size_t j = 0; j < polygons.size(); j++) | |
| 206 { | |
| 207 PolylineSceneLayer::Chain chain; | |
| 208 chain.resize(polygons[j].size()); | |
| 209 | |
| 210 for (size_t k = 0; k < polygons[j].size(); k++) | |
| 211 { | |
| 212 chain[k] = ScenePoint2D(polygons[j][k].first, polygons[j][k].second); | |
| 213 } | |
| 214 | |
| 215 layer->AddChain(chain, true /* closed */, color); | |
| 216 } | |
| 217 } | |
| 218 } | |
| 219 | |
| 220 return layer.release(); | |
| 221 } | |
| 222 }; | |
| 223 | |
| 224 | |
| 225 DicomStructureSetLoader::DicomStructureSetLoader(IOracle& oracle, | |
| 226 IObservable& oracleObservable) : | |
| 227 LoaderStateMachine(oracle, oracleObservable), | |
| 228 revision_(0), | |
| 229 countProcessedInstances_(0), | |
| 230 countReferencedInstances_(0) | |
| 231 { | |
| 232 } | |
| 233 | |
| 234 | |
| 235 void DicomStructureSetLoader::LoadInstance(const std::string& instanceId) | |
| 236 { | |
| 237 Start(); | |
| 238 | |
| 239 instanceId_ = instanceId; | |
| 240 | |
| 241 { | |
| 242 std::auto_ptr<OrthancRestApiCommand> command(new OrthancRestApiCommand); | |
| 243 command->SetHttpHeader("Accept-Encoding", "gzip"); | |
| 244 command->SetUri("/instances/" + instanceId + "/tags?ignore-length=3006-0050"); | |
| 245 command->SetPayload(new LoadStructure(*this)); | |
| 246 Schedule(command.release()); | |
| 247 } | |
| 248 } | |
| 249 | |
| 250 | |
| 251 IVolumeSlicer::IExtractedSlice* DicomStructureSetLoader::ExtractSlice(const CoordinateSystem3D& cuttingPlane) | |
| 252 { | |
| 253 if (content_.get() == NULL) | |
| 254 { | |
| 255 // Geometry is not available yet | |
| 256 return new IVolumeSlicer::InvalidSlice; | |
| 257 } | |
| 258 else | |
| 259 { | |
| 260 return new Slice(*content_, revision_, cuttingPlane); | |
| 261 } | |
| 262 } | |
| 263 } |