Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1106 1739 (N39E17 L= 62)   613     78        2 (1062/1714-1149/1762)   0.0 E20-E13 N36-N40
      2 1324 1540 (N31W04 L= 83)  2285    707        8 (1230/1244-1406/1721) -53.1 E06-W15 N14-N42
      3 1519 1634 (N46W18 L= 96)   417     51        0 (1492/1607-1547/1659) -21.8 W14-W20 N44-N47
      4  906 1530 (N19E20 L= 58)   803    211        0 ( 850/1448- 983/1618) -31.6 E26-E13 N15-N23
      5  257 1472 (N03E68 L= 11)  1403    157        3 ( 230/1417- 288/1517) -80.5 E73-E63 N00-N06
      6 1775 1213 (N25W44 L=122)  1534    251        2 (1685/1128-1853/1309) -15.5 W34-W52 N22-N28
      7 1552 1030 (N07W30 L=109)   414    150        1 (1501/ 981-1609/1060)   0.0 W27-W35 N04-N08
      8  251  926 (S30E51 L= 28)   648    166        1 ( 201/ 899- 305/ 946)  11.3 E55-E45 S31-S27
      9  729  451 (S50E03 L= 75)   576    124        1 ( 685/ 422- 797/ 492)  -5.7 E07-W03 S50-S46
     10  863  211 (S65W26 L=105)   457     85        1 ( 835/ 186- 885/ 233)  45.0 W22-W29 S66-S61
  

    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