Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  676 1524 (N18E31 L=246)   439     89        1 ( 638/1500- 728/1544)  21.8 E33-E27 N16-N19
      2 1699 1448 (N33W41 L=317)  1027    219        1 (1630/1349-1759/1523) -42.0 W34-W46 N27-N36
      3  357 1374 (N05E50 L=227)   432    105        1 ( 338/1319- 369/1418) -50.2 E52-E47 N01-N07
      4  150 1169 (S08E67 L=209)   830    151        0 ( 133/1099- 162/1245) -83.7 E68-E65 S12-S03
      5 1845 1156 (N17W55 L=331)   629    108        1 (1794/1131-1898/1174)  18.4 W49-W61 N15-N20
      6 1049 1133 (S00E00 L=276)   226     86        0 (1004/1118-1082/1150)   0.0 E03-W01 N00-N01
      7 1891  974 (N08W62 L=338)  4431    221        9 (1828/ 875-1956/1065)  67.4 W55-W72 N01-N13
      8 1009  717 (S25W04 L=281)   419    160        0 ( 959/ 649-1064/ 778)  68.2 W02-W06 S29-S19
      9 1132  597 (S30W15 L=291)  2063    505        2 ( 989/ 486-1340/ 682)  31.8 W06-W27 S38-S22
     10 1655  450 (S26W62 L=338)  2593    242        3 (1597/ 374-1716/ 545) -21.8 W53-W78 S30-S20
     11 1109  367 (S45W22 L=299)  3001    730        4 ( 831/ 261-1443/ 472)   1.9 E02-W59 S47-S41
     12 1739  430 (S25W75 L=352)   740     54        2 (1712/ 407-1766/ 458)  15.9 W70-W81 S26-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