Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  532 1674 (N25E56 L=327)   939    108        1 ( 483/1636- 572/1701)  -4.8 E61-E49 N23-N26
      2  385 1499 (N11E57 L=326)   969    170        3 ( 339/1430- 422/1557) -32.5 E61-E50 N07-N14
      3 1248 1397 (N24W05 L= 28)   277     70        0 (1212/1373-1291/1416)   0.0 W02-W07 N23-N25
      4 1393 1108 (N11W20 L= 43)   338    113        0 (1371/1052-1409/1163) -63.4 W17-W21 N07-N13
      5 1831  965 (N14W58 L= 81)   882    137        0 (1785/ 894-1857/1006) -21.8 W52-W62 N09-N15
      6  310  831 (S30E47 L=336)   448     70        0 ( 289/ 788- 334/ 866)  68.2 E49-E44 S33-S28
      7 1084  614 (S25W14 L= 37)   234     53        0 (1054/ 598-1113/ 631)   0.0 W11-W16 S26-S24
      8 1110  456 (S35W22 L= 45)   765    149        2 (1041/ 418-1171/ 498)  -4.8 W15-W27 S35-S32
      9 1545  437 (S24W55 L= 78)   682     93        0 (1500/ 418-1592/ 457)  12.5 W49-W59 S24-S21
     10 1224  262 (S45W45 L= 68)   490     96        0 (1196/ 245-1273/ 289)  18.4 W40-W48 S46-S42
  

    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