Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1682 1438 (N21W49 L=101)  2155    348        6 (1572/1350-1804/1530) -29.1 W41-W60 N15-N27
      2  454 1331 (N19E37 L= 15)   312    140        0 ( 418/1266- 473/1374) -54.5 E40-E35 N14-N21
      3  290 1255 (N15E51 L=  1)   400     98        0 ( 262/1215- 325/1305)  45.0 E53-E48 N12-N17
      4 1814  907 (S11W56 L=108)   996    167        2 (1771/ 840-1859/ 972)  48.8 W52-W60 S15-S07
      5 1712  678 (S25W50 L=102)   617    115        0 (1683/ 633-1741/ 741)  63.4 W48-W51 S28-S21
      6  579  618 (S26E33 L= 19)   231     34        0 ( 559/ 594- 600/ 642) -76.0 E33-E31 S27-S23
      7 1260  581 (S31W15 L= 67)   430    145        0 (1226/ 514-1317/ 622)  49.4 W12-W18 S35-S28
      8  503  557 (S29E41 L= 11)   742    157        0 ( 433/ 526- 578/ 586)   8.1 E47-E33 S30-S27
      9  885  325 (S48E16 L= 36)   606     91        0 ( 842/ 298- 933/ 357) -39.8 E19-E12 S50-S45
     10 1012  321 (S49E05 L= 47)   365     64        1 ( 978/ 306-1051/ 349)   9.5 E07-E01 S50-S46
     11  686  286 (S50E38 L= 14)   822    173        1 ( 614/ 240- 741/ 347)  19.2 E51-E28 S54-S45
  

    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