orthanc-stone: OrthancStone/Sources/Scene2D/AnnotationsSceneLayer.cpp comparison

comparison OrthancStone/Sources/Scene2D/AnnotationsSceneLayer.cpp @ 1984:187a261d7ae2

computation of mean and stddev in rectangle probes

author	Sebastien Jodogne <s.jodogne@gmail.com>
date	Mon, 31 Oct 2022 14:42:46 +0100
parents	20fa913272b7
children	bb307007f8e2

comparison

equal deleted inserted replaced

-:20fa913272b7
+:187a261d7ae2
 Segment&  segment2_;
 Segment&  segment3_;
 Segment&  segment4_;
 Text&     label_;
-void UpdateLabel()
+protected:
-{
+virtual void UpdateProbeForLayer(const ISceneLayer& layer) ORTHANC_OVERRIDE
-TextSceneLayer content;
+{
+double x1 = handle1_.GetCenter().GetX();
-const double x1 = handle1_.GetCenter().GetX();
+double y1 = handle1_.GetCenter().GetY();
-const double y1 = handle1_.GetCenter().GetY();
+double x2 = handle2_.GetCenter().GetX();
-const double x2 = handle2_.GetCenter().GetX();
+double y2 = handle2_.GetCenter().GetY();
-const double y2 = handle2_.GetCenter().GetY();
 // Put the label to the right of the right-most handle
 //const double y = std::min(y1, y2);
 const double y = (y1 + y2) / 2.0;
 if (x1 < x2)
 {
-content.SetPosition(x2, y);
+label_.SetPosition(x2, y);
 }
 else
 {
-content.SetPosition(x1, y);
+label_.SetPosition(x1, y);
 }
-content.SetAnchor(BitmapAnchor_CenterLeft);
+std::string text;
-content.SetBorder(10);
-const double area = std::abs(x1 - x2) * std::abs(y1 - y2);
 char buf[32];
-switch (GetUnits())
+if (GetUnits() == Units_Millimeters)
 {
-case Units_Millimeters:
+const double area = std::abs(x1 - x2) * std::abs(y1 - y2);
-sprintf(buf, "%0.2f cm%c%c",
-area / 100.0,
+sprintf(buf, "%0.2f cm%c%c",
-0xc2, 0xb2 /* two bytes corresponding to two power in UTF-8 */);
+area / 100.0,
-break;
+0xc2, 0xb2 /* two bytes corresponding to two power in UTF-8 */);
+text = buf;
-case Units_Pixels:
+}
-// Don't report area (pixel-times-pixel is a strange unit)
-sprintf(buf, "hello");
+if (layer.GetType() == ISceneLayer::Type_FloatTexture)
-break;
+{
+const TextureBaseSceneLayer& texture = dynamic_cast<const TextureBaseSceneLayer&>(layer);
+const AffineTransform2D sceneToTexture = AffineTransform2D::Invert(texture.GetTransform());
+const Orthanc::ImageAccessor& image = texture.GetTexture();
+assert(image.GetFormat() == Orthanc::PixelFormat_Float32);
+sceneToTexture.Apply(x1, y1);
+sceneToTexture.Apply(x2, y2);
+int ix1 = static_cast<int>(std::floor(x1));
+int iy1 = static_cast<int>(std::floor(y1));
+int ix2 = static_cast<int>(std::floor(x2));
+int iy2 = static_cast<int>(std::floor(y2));
+if (ix1 > ix2)
+{
+std::swap(ix1, ix2);
+}
+if (iy1 > iy2)
+{
+std::swap(iy1, iy2);
+}
+LinearAlgebra::OnlineVarianceEstimator estimator;
+for (int y = std::max(0, iy1); y <= std::min(static_cast<int>(image.GetHeight()) - 1, iy2); y++)
+{
+int x = std::max(0, ix1);
-default:
+const float* p = reinterpret_cast<const float*>(image.GetConstRow(y)) + x;
-throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
-}
+for (; x <= std::min(static_cast<int>(image.GetWidth()) - 1, ix2); x++, p++)
+{
-content.SetText(buf);
+estimator.AddSample(*p);
+}
-label_.SetContent(content);
+}
-}
+if (estimator.GetCount() > 0)
-protected:
+{
-virtual void UpdateProbeForLayer(const ISceneLayer& layer) ORTHANC_OVERRIDE
+if (!text.empty())
 {
-// TODO
+text += "\n";
+}
+sprintf(buf, "Mean: %0.1f\nStdDev: %0.1f", estimator.GetMean(), estimator.GetStandardDeviation());
+text += buf;
+}
+}
+label_.SetText(text);
 }
 public:
 RectangleProbeAnnotation(AnnotationsSceneLayer& that,
 Units units,
 segment2_(AddTypedPrimitive<Segment>(new Segment(*this, p2.GetX(), p1.GetY(), p2.GetX(), p2.GetY()))),
 segment3_(AddTypedPrimitive<Segment>(new Segment(*this, p2.GetX(), p2.GetY(), p1.GetX(), p2.GetY()))),
 segment4_(AddTypedPrimitive<Segment>(new Segment(*this, p1.GetX(), p2.GetY(), p1.GetX(), p1.GetY()))),
 label_(AddTypedPrimitive<Text>(new Text(that, *this)))
 {
+TextSceneLayer content;
+content.SetAnchor(BitmapAnchor_CenterLeft);
+content.SetBorder(10);
+content.SetText("?");
+label_.SetContent(content);
 label_.SetColor(COLOR_TEXT);
-UpdateLabel();
 }
 virtual unsigned int GetHandlesCount() const ORTHANC_OVERRIDE
 {
 return 2;
 else
 {
 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
 }
-UpdateLabel();
+TagProbeAsChanged();
-}
-virtual void UpdateProbe(const Scene2D& scene) ORTHANC_OVERRIDE
-{
 }
 virtual void Serialize(Json::Value& target) ORTHANC_OVERRIDE
 {
 target = Json::objectValue;

Mercurial > hg > orthanc-stone

comparison OrthancStone/Sources/Scene2D/AnnotationsSceneLayer.cpp @ 1984:187a261d7ae2