Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  600 1534 (N30E37 L=211)  1517    378        2 ( 533/1408- 703/1627) -31.8 E47-E25 N22-N36
      2 1583 1517 (N38W42 L=290)   949    132        1 (1530/1463-1641/1564) -36.9 W37-W46 N34-N41
      3 1503  925 (S01W32 L=279)   615    155        1 (1431/ 906-1580/ 956)  10.3 W26-W37 S02-N01
      4 1825  927 (N01W60 L=308)   534     40        0 (1804/ 900-1844/ 956)  90.0 W57-W62 N00-N03
      5 1751  851 (S04W53 L=301)   555     35        0 (1727/ 818-1775/ 878) -71.6 W50-W55 S05-S02
      6  400  832 (S16E41 L=206)   344     71        1 ( 370/ 812- 429/ 862) -26.6 E44-E39 S16-S13
      7 1827  774 (S08W64 L=312)   501     80        0 (1811/ 739-1846/ 811)  78.7 W61-W66 S10-S05
      8 1435  721 (S14W30 L=277)   961    277        2 (1328/ 654-1533/ 779) -19.4 W21-W38 S17-S10
      9 1243  686 (S17W17 L=265)   455     90        2 (1195/ 657-1283/ 717)  21.8 W13-W19 S19-S15
     10  646  655 (S25E23 L=225)   277     99        0 ( 598/ 645- 700/ 663)   7.1 E26-E18 S25-S23
     11 1456  532 (S26W36 L=284)   296    106        0 (1425/ 490-1481/ 585)  81.9 W34-W37 S29-S22
  

    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