Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  702 1731 (N31E44 L=352)  1200    224        0 ( 660/1675- 760/1790) -18.4 E51-E36 N28-N35
      2  651 1590 (N21E40 L=357)  1714    277        4 ( 539/1521- 723/1635)  -6.7 E49-E32 N17-N23
      3 1190 1486 (N28E02 L= 34)  1365    404        1 (1025/1354-1366/1617) -14.0 E16-W12 N23-N33
      4  482 1436 (N08E45 L=351)   545    149        2 ( 445/1396- 521/1480) -26.6 E49-E41 N05-N09
      5 1638 1413 (N37W32 L= 68)   290     64        0 (1612/1388-1666/1440) -18.4 W29-W35 N35-N37
      6  820 1371 (N11E21 L= 16)   301     79        0 ( 793/1328- 847/1409) -38.7 E23-E17 N09-N13
      7 1751 1291 (N32W43 L= 80)   514    100        1 (1710/1250-1780/1314) -15.9 W39-W46 N30-N33
      8  706  931 (S16E16 L= 20)   318     97        1 ( 674/ 882- 741/ 985) -39.8 E19-E13 S18-S13
      9  865  498 (S39W05 L= 41)   404     88        1 ( 823/ 475- 913/ 527)   8.1 W01-W08 S40-S37
  

    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