Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1787 1271 (N14W59 L=313)  1493    249        4 (1699/1181-1851/1341) -50.7 W52-W66 N08-N19
      2  335 1162 (N14E49 L=205)   329     72        1 ( 297/1144- 365/1177)   0.0 E52-E45 N13-N15
      3 1118  868 (S06W05 L=258)   277    125        0 (1057/ 834-1172/ 888) -23.2 W01-W08 S07-S04
      4  440  842 (S05E40 L=214)   279     82        3 ( 423/ 813- 462/ 886) -63.4 E41-E38 S07-S02
      5  208  661 (S18E71 L=182)   794     43        1 ( 184/ 637- 236/ 678)   9.5 E77-E65 S19-S16
      6 1462  641 (S22W28 L=282)   552     76        1 (1434/ 603-1495/ 675) -68.2 W26-W30 S24-S19
      7  595  599 (S21E32 L=222)   261     45        0 ( 571/ 575- 615/ 621) -71.6 E33-E30 S22-S19
      8 1322  523 (S30W19 L=272)   588    119        0 (1277/ 489-1379/ 556)  29.7 W15-W23 S31-S27
      9  574  518 (S27E36 L=218)   326    119        0 ( 549/ 470- 615/ 554)  40.6 E38-E31 S30-S24
     10 1025  519 (S28E02 L=251)   271    101        0 ( 976/ 505-1073/ 545)  18.4 E05-W01 S29-S26
     11  910  472 (S31E11 L=242)   298     65        1 ( 882/ 443- 934/ 506)  51.3 E13-E09 S33-S28
     12 1591  386 (S44W51 L=304)   560     58        0 (1568/ 357-1614/ 413) -24.8 W45-W58 S46-S40
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y