Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1233 1795 (N58E03 L=326)   373     54        1 (1200/1780-1263/1811)  38.7 E05-W00 N55-N60
      2  684 1750 (N39E48 L=281)   389     36        3 ( 661/1730- 712/1773)   0.0 E50-E44 N37-N41
      3 1479 1710 (N57W28 L=357)   393     45        0 (1462/1685-1497/1736)  50.2 W24-W31 N53-N59
      4 1644 1650 (N56W52 L= 21)   839     75        2 (1605/1619-1670/1671) -31.0 W46-W62 N52-N58
      5  895 1599 (N33E22 L=307)   363     89        2 ( 856/1568- 931/1633) -23.2 E24-E17 N30-N34
      6 1383  781 (S09W27 L=356)   358    120        0 (1322/ 759-1439/ 803)  20.6 W22-W30 S10-S06
      7  449  697 (S32E35 L=294)   389     48        1 ( 423/ 672- 476/ 724)  53.1 E37-E32 S33-S29
      8  991  536 (S32W08 L=337)  2576    658        3 ( 792/ 432-1251/ 687) -11.0 E08-W28 S37-S24
  

    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