Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1400 1509 (N34W29 L= 89)   496     98        0 (1371/1457-1423/1552) -90.0 W27-W30 N30-N37
      2 1628 1487 (N31W48 L=108)   775    107        1 (1592/1439-1664/1534) -63.4 W45-W50 N28-N34
      3  296 1480 (N32E60 L=  0)   880    220        2 ( 265/1445- 339/1515) -21.8 E64-E54 N29-N33
      4  544 1394 (N27E34 L= 26)   251     87        1 ( 505/1380- 582/1409)   0.0 E37-E31 N26-N27
      5 1884 1182 (N10W69 L=129)   893    145        1 (1856/1137-1915/1228) -40.6 W64-W74 N07-N13
      6 1757 1150 (N09W53 L=113)   471    133        1 (1722/1111-1800/1194) -40.6 W49-W56 N06-N12
      7 1069  961 (S00W03 L= 63)   285     88        0 (1023/ 941-1121/ 984)  21.8 W00-W05 S01-N01
      8 1821  822 (S11W59 L=119)   712    173        0 (1799/ 788-1860/ 876) -40.6 W56-W64 S13-S07
      9  453  785 (S11E39 L= 21)   803    254        0 ( 347/ 728- 546/ 855) -25.0 E47-E31 S14-S07
     10  344  625 (S22E51 L=  9)  1428    433        5 ( 263/ 534- 403/ 713)  35.8 E63-E44 S29-S16
     11  992  577 (S25E03 L= 57)   245     47        0 ( 962/ 563-1022/ 593)  26.6 E04-W00 S25-S23
  

    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