Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  525 1726 (N47E52 L=137)  1286    311        3 ( 481/1672- 559/1764) -23.2 E60-E46 N43-N50
      2 1101 1568 (N39W03 L=193)   231     71        2 (1074/1537-1128/1593) -45.0 W01-W05 N36-N40
      3  258 1283 (N15E59 L=131)  4149    719       10 ( 181/1159- 341/1464)  78.7 E68-E52 N06-N26
      4  369 1230 (N12E47 L=143)   361     84        1 ( 351/1196- 396/1262)  63.4 E48-E44 N10-N14
      5  248  854 (S12E56 L=134)  1762    193        4 ( 206/ 786- 300/ 928) -63.4 E61-E51 S15-S07
      6 1761  830 (S08W54 L=244)   398     72        0 (1740/ 796-1777/ 861) -53.1 W51-W55 S09-S05
      7 1099  451 (S35W08 L=198)   503     77        0 (1056/ 432-1147/ 472)   0.0 W05-W11 S36-S33
      8  720  418 (S39E21 L=168)  1151    285        4 ( 600/ 379- 831/ 447)  15.9 E33-E12 S43-S36
  

    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