Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  912 1906 (N48E52 L=220)   601     74        0 ( 891/1875- 938/1925) -36.9 E56-E47 N44-N50
      2  786 1832 (N37E50 L=222)   783     99        1 ( 765/1782- 805/1887) -39.8 E57-E44 N32-N42
      3 1092 1638 (N32E16 L=256)   643    134        2 (1061/1576-1116/1706) -63.4 E18-E13 N27-N37
      4 1351 1529 (N33W05 L=276)  1172    346        1 (1274/1372-1421/1659) -46.8 E04-W12 N24-N40
      5 1476 1172 (N15W22 L=293)   905    232        3 (1420/1078-1515/1281) -51.3 W17-W25 N10-N21
      6  487 1155 (S13E36 L=236)   491    107        1 ( 466/1097- 511/1212) -56.3 E38-E33 S15-S09
      7  740 1013 (S14E16 L=256)   188    113        0 ( 712/ 968- 780/1048)  71.6 E16-E14 S17-S11
      8 1449  990 (N04W25 L=297)   352    112        0 (1409/ 949-1481/1035)  71.6 W23-W25 N01-N07
      9  986  415 (S42W20 L=292)   610    206        1 ( 886/ 396-1080/ 435)  18.4 W12-W27 S45-S39
     10 1306  407 (S31W41 L=313)   641    131        1 (1248/ 390-1372/ 432)  10.3 W35-W46 S33-S28
     11 1557  299 (S27W74 L=346)   812    116        0 (1509/ 282-1611/ 315)  26.6 W67-W80 S29-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