Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1035 1597 (N35W10 L=151)   531    160        1 ( 966/1570-1113/1616) -21.8 W05-W15 N32-N37
      2  889 1571 (N36E01 L=140)   302     80        1 ( 863/1535- 923/1608)  59.0 E03-W02 N33-N38
      3  782 1508 (N33E09 L=132)   281     73        1 ( 746/1483- 820/1529)  11.3 E12-E06 N31-N34
      4 1809  954 (S15W54 L=195)  1069    184        3 (1775/ 879-1865/1037)  41.6 W50-W60 S19-S10
      5 1617  991 (S10W37 L=178)   261     81        2 (1599/ 953-1639/1024)  63.4 W35-W39 S11-S07
      6 1561  807 (S20W32 L=173)   714    227        1 (1480/ 731-1639/ 895)  33.7 W25-W38 S23-S15
      7 1742  717 (S29W49 L=190)   340     63        1 (1724/ 684-1762/ 741) -63.4 W46-W50 S30-S26
      8  558  688 (S14E35 L=107)   487    106        3 ( 516/ 641- 587/ 733)  45.0 E37-E31 S16-S11
      9 1019  633 (S23E06 L=135)   471    234        0 ( 969/ 529-1080/ 718)  45.0 E11-W00 S29-S18
     10  326  567 (S18E59 L= 82)  2490    378       10 ( 234/ 444- 453/ 663) -60.6 E66-E52 S27-S11
     11  737  439 (S32E30 L=111)   391    106        1 ( 704/ 401- 773/ 488)  29.1 E34-E25 S34-S29
     12  619  291 (S41E50 L= 91)   433     51        2 ( 595/ 261- 639/ 316)  76.0 E52-E47 S43-S38
  

    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