Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1552 1735 (N30W65 L= 73)  1300    144        1 (1507/1686-1613/1772)   0.0 W56-W76 N27-N31
      2 1809 1211 (S05W56 L= 65)  1004    161        3 (1783/1124-1834/1279)  63.4 W52-W59 S09-S01
      3 1910 1221 (S08W69 L= 78)  1538    110        5 (1875/1166-1938/1270) -54.5 W64-W76 S11-S04
      4  438 1223 (N27E33 L=336)   277     82        0 ( 417/1190- 472/1263) -51.3 E34-E30 N24-N29
      5  206 1144 (N27E56 L=312)   397     48        0 ( 176/1127- 237/1160)   0.0 E60-E52 N26-N28
      6 1337 1022 (S04W17 L= 26)  1681    502        3 (1188/ 845-1462/1157) -65.8 W10-W21 S16-N05
      7 1602 1089 (S07W36 L= 44)   575    138        2 (1559/1022-1643/1141) -82.9 W33-W38 S10-S02
      8 1702 1071 (S10W43 L= 52)   369     83        3 (1663/1045-1736/1097)   0.0 W39-W46 S11-S09
      9  755 1045 (N10E15 L=354)   514    154        1 ( 677/1024- 832/1070) -26.6 E18-E10 N08-N13
     10 1667  987 (S14W38 L= 47)   700    164        3 (1608/ 939-1724/1039)  11.3 W33-W43 S15-S11
     11 1748  806 (S27W43 L= 52)  1487    272        1 (1668/ 774-1854/ 882) -11.9 W36-W55 S30-S20
     12 1563  738 (S26W27 L= 36)   390    116        0 (1548/ 686-1577/ 802)  82.9 W25-W27 S29-S21
     13  897  733 (S10E14 L=355)   337     58        2 ( 868/ 694- 922/ 760) -76.0 E15-E11 S12-S08
     14  949  714 (S13E11 L=357)   286     67        1 ( 921/ 693- 981/ 739) -33.7 E12-E09 S13-S10
     15  758  467 (S23E31 L=338)  3440    825        4 ( 518/ 247-1022/ 691)  15.3 E64-E08 S33-S15
     16 1574  552 (S39W27 L= 35)  2452    270        6 (1525/ 419-1649/ 652)  51.8 W21-W33 S48-S33
     17  494  586 (S10E43 L=325)   419    138        0 ( 426/ 567- 548/ 620)   0.0 E49-E38 S11-S09
  

    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