Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  555 1680 (N20E52 L=131)  1732    311        4 ( 506/1615- 629/1748) -23.2 E61-E47 N15-N23
      2  832 1361 (N07E20 L=163)   232     51        3 ( 804/1344- 859/1383)  33.7 E21-E18 N05-N08
      3 1709 1198 (N25W38 L=221)   382    101        1 (1676/1158-1744/1233) -15.9 W34-W41 N23-N25
      4  855  917 (S16E07 L=177)   693    202        0 ( 765/ 857- 945/ 977)  59.7 E10-E02 S22-S10
      5 1816  937 (N14W54 L=237)   371     60        2 (1793/ 915-1837/ 961)  63.4 W51-W55 N11-N15
      6  206  886 (S35E55 L=129)   661     85        2 ( 182/ 852- 232/ 916)  90.0 E57-E51 S36-S32
      7 1575  707 (S07W42 L=226)  1197    212        6 (1541/ 609-1613/ 812) -67.4 W39-W45 S12-N00
      8 1331  780 (S10W24 L=208)   347     81        1 (1301/ 744-1364/ 807) -21.8 W22-W27 S12-S09
      9 1123  686 (S22W16 L=199)   452    150        0 (1048/ 656-1195/ 724)   5.7 W10-W20 S22-S20
     10  669  578 (S42E10 L=174)   986    247        0 ( 545/ 558- 794/ 608)  19.3 E20-W00 S45-S38
     11 1494  388 (S26W54 L=237)  1471    304        4 (1355/ 360-1615/ 421)  34.5 W44-W62 S31-S20
  

    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