Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  617 1756 (N32E50 L=286)  1674    281        4 ( 556/1674- 662/1810) -35.5 E57-E40 N26-N36
      2  495 1708 (N27E57 L=278)   488     91        2 ( 473/1664- 530/1743) -71.6 E60-E54 N23-N29
      3  546 1676 (N25E50 L=286)  1135    209        3 ( 468/1614- 601/1716)  60.9 E53-E45 N19-N28
      4  581 1480 (N12E37 L=298)  1438    270        4 ( 504/1363- 665/1597)  83.3 E40-E35 N04-N21
      5  728 1458 (N13E27 L=309)   483    109        1 ( 702/1407- 757/1501) -80.5 E28-E25 N10-N16
      6 1545 1421 (N28W26 L=  1)   249     73        0 (1527/1383-1571/1454) -45.0 W24-W28 N26-N30
      7  923 1341 (N10E12 L=323)   546    145        2 ( 877/1278- 972/1413)  90.0 E13-E09 N05-N14
      8 1839 1221 (N24W54 L= 30)   574     63        1 (1811/1185-1864/1260) -45.0 W51-W57 N21-N26
      9  649 1229 (S02E26 L=309)   252     87        0 ( 613/1193- 674/1254) -21.8 E28-E23 S03-N00
     10 1544 1217 (N15W28 L=  3)   444    118        1 (1495/1168-1573/1249) -18.4 W24-W30 N12-N16
     11 1369 1201 (N10W17 L=352)   280     65        1 (1350/1163-1392/1231) -56.3 W15-W17 N08-N11
     12  200 1065 (S17E60 L=275)  1692    227        3 ( 160/ 969- 259/1153)  65.6 E65-E54 S22-S11
     13 1576 1051 (N06W33 L=  8)   329     69        2 (1540/1028-1611/1083)  11.3 W30-W35 N05-N08
     14 1774 1011 (N09W49 L= 25)   558    134        0 (1755/ 943-1797/1063)  90.0 W47-W50 N04-N12
     15  568  803 (S28E25 L=310)  1589    538        3 ( 380/ 682- 685/ 999)  46.2 E42-E19 S38-S14
     16 1818  881 (N03W57 L= 32)   502     50        1 (1796/ 848-1834/ 911)  90.0 W54-W58 N01-N05
     17 1780  668 (S09W62 L= 37)   501     80        0 (1744/ 642-1812/ 698)  38.7 W57-W64 S11-S06
     18 1585  486 (S25W54 L= 29)  2063    345        2 (1528/ 329-1645/ 634) -51.7 W48-W64 S34-S15
  

    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