Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1661 1551 (N20W58 L=157)   406     47        0 (1637/1527-1683/1575) -56.3 W56-W60 N18-N22
      2  589 1388 (N31E22 L= 77)   596    228        0 ( 533/1308- 668/1433)  29.1 E26-E16 N26-N32
      3 1838 1042 (S12W57 L=156)   413     62        1 (1808/1009-1860/1065) -45.0 W54-W59 S14-S10
      4 1564 1017 (S09W33 L=132)   354     90        1 (1545/ 965-1586/1058)  68.2 W31-W34 S11-S06
      5  332  872 (N04E48 L= 51)   901    171        0 ( 264/ 834- 385/ 938)  18.4 E54-E42 N02-N07
      6 1115  853 (S10W02 L=101)   685    189        1 (1067/ 769-1168/ 925)  53.1 E01-W06 S14-S06
      7  526  765 (S05E35 L= 64)   411     61        1 ( 507/ 732- 562/ 794)  11.3 E37-E32 S06-S03
      8 1838  710 (S33W60 L=158)   846    119        0 (1806/ 663-1862/ 754) -36.0 W53-W66 S37-S29
      9 1728  670 (S34W46 L=145)   754    159        1 (1679/ 621-1779/ 748)  33.7 W41-W50 S36-S29
     10  532  619 (S14E39 L= 60)   622     90        2 ( 493/ 583- 563/ 670)  38.7 E42-E35 S15-S11
  

    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