Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1344 1745 (N54W35 L=210)   352     30        0 (1328/1722-1362/1766)  29.7 W31-W39 N51-N55
      2 1394 1740 (N53W41 L=216)   514     70        0 (1367/1715-1426/1759)  14.0 W35-W47 N50-N54
      3  975 1729 (N53E05 L=170)   313     33        0 ( 954/1708- 993/1751)  53.1 E07-E03 N50-N54
      4  855 1732 (N53E18 L=157)   335     68        0 ( 826/1717- 887/1745)   9.5 E20-E14 N51-N54
      5  684 1551 (N38E28 L=147)   637    180        1 ( 609/1498- 745/1582) -19.7 E35-E21 N33-N39
      6 1726 1548 (N36W68 L=243)  1051    104        1 (1682/1526-1777/1563)   0.0 W60-W82 N34-N36
      7  851 1492 (N34E13 L=162)   497    141        1 ( 793/1466- 921/1517)  16.7 E17-E07 N31-N35
      8 1057 1460 (N31W02 L=177)   275     46        0 (1027/1445-1091/1477) -14.0 W00-W04 N30-N32
      9  243 1255 (N16E61 L=114)   618    104        1 ( 221/1216- 283/1302) -24.0 E64-E55 N13-N18
     10  370 1277 (N18E48 L=127)   428    133        0 ( 327/1244- 427/1301)  12.5 E51-E42 N15-N19
     11  377 1141 (N10E45 L=130)   409     68        1 ( 350/1113- 415/1171) -18.4 E47-E41 N08-N11
     12 1231  673 (S19W14 L=189)   330    134        0 (1170/ 656-1298/ 693)  12.5 W09-W18 S19-S17
     13  663  487 (S32E27 L=148)  3242    505        5 ( 504/ 345- 833/ 663)  32.7 E51-E12 S44-S19
     14  376  524 (S31E54 L=121)  1335    228        0 ( 345/ 420- 429/ 617)  32.6 E69-E44 S39-S23
     15 1544  530 (S29W40 L=215)   623    142        1 (1484/ 488-1599/ 572)  36.9 W35-W43 S32-S26
     16 1399  417 (S38W31 L=206)   400     63        1 (1370/ 389-1426/ 451)  78.7 W28-W32 S39-S34
     17 1428  323 (S46W39 L=214)   541     72        1 (1398/ 297-1445/ 359)  90.0 W36-W42 S48-S42
     18 1483  299 (S49W49 L=224)   533     33        0 (1462/ 282-1508/ 322) -21.8 W44-W54 S50-S46
  

    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