Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1475 1723 (N31W58 L=156)   396     44        0 (1453/1704-1504/1738)   0.0 W54-W61 N29-N32
      2 1216 1456 (N24W24 L=122)   229     35        1 (1189/1437-1244/1470) -14.0 W21-W25 N22-N24
      3 1660 1381 (N07W50 L=148)   840    188        0 (1599/1297-1704/1450) -84.8 W47-W52 N01-N12
      4  711 1188 (N22E14 L= 84)  1003    212        1 ( 614/1163- 819/1218) -21.0 E20-E07 N19-N25
      5 1835 1086 (S14W55 L=153)   338     39        1 (1822/1059-1855/1109)  26.6 W53-W57 S15-S12
      6 1618  932 (S16W33 L=131)   846    315        0 (1525/ 829-1706/1049)  22.4 W24-W41 S19-S11
      7 1634  766 (S26W31 L=129)  1844    556        2 (1510/ 562-1744/1011)  38.7 W19-W44 S36-S14
      8  189  911 (N17E62 L= 36)  2526    277        7 (  90/ 864- 324/ 984) -16.7 E78-E48 N12-N22
      9 1320  972 (S05W15 L=113)   183     76        0 (1279/ 960-1355/ 983) -26.6 W12-W16 S06-S04
     10  521  830 (N06E34 L= 64)   546    145        1 ( 497/ 772- 556/ 889)  74.1 E36-E31 N02-N09
     11 1819  747 (S34W50 L=148)   521    129        2 (1784/ 705-1860/ 815)  24.0 W45-W55 S35-S30
     12  378  764 (N05E47 L= 51)   384     72        0 ( 342/ 734- 419/ 795) -51.3 E49-E44 N02-N07
     13  794  701 (S07E21 L= 77)   318    109        0 ( 749/ 673- 842/ 735)   8.1 E24-E17 S08-S06
     14  376  602 (S04E53 L= 45)   887    180        1 ( 322/ 557- 430/ 675)  19.7 E59-E45 S06-N00
     15 1465  411 (S46W13 L=111)   524     53        1 (1431/ 384-1499/ 433) -45.0 W10-W15 S48-S42
     16 1153  315 (S42E14 L= 84)   380     84        1 (1108/ 290-1193/ 342)   7.1 E18-E10 S43-S40
     17 1025  224 (S45E29 L= 69)   381     38        2 ( 998/ 200-1051/ 249) -11.3 E31-E26 S47-S43
     18  826  166 (S44E51 L= 48)   677     72        0 ( 787/ 152- 865/ 180)   0.0 E55-E45 S45-S42
  

    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