Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  512 1589 (N56E29 L=180)   416     48        2 ( 494/1557- 531/1613)  78.7 E31-E25 N53-N58
      2  331 1192 (N32E46 L=163)   296     63        0 ( 301/1180- 369/1209) -14.0 E49-E41 N30-N33
      3 1053 1122 (N12W04 L=213)   190     42        0 (1024/1108-1080/1138) -45.0 W02-W05 N10-N13
      4 1409 1069 (S00W23 L=232)   333    111        0 (1354/1056-1470/1081) -18.4 W19-W26 S01-N01
      5 1601  815 (S20W30 L=239)   759    296        1 (1501/ 695-1698/ 881)  16.4 W21-W39 S24-S17
      6 1041  778 (S07E05 L=204)   365    140        2 (1003/ 714-1083/ 836) -82.9 E06-E03 S11-S03
      7 1641  727 (S27W32 L=241)   820    222        1 (1602/ 632-1685/ 817)  51.3 W28-W37 S32-S22
      8  963  532 (S20E16 L=192)  1396    427        4 ( 817/ 416-1078/ 647)  21.4 E29-E06 S24-S14
      9 1152  478 (S29E07 L=202)   378    171        0 (1111/ 422-1229/ 526)  20.0 E11-W00 S31-S27
     10  456  407 (S18E59 L=150)  1111     97        2 ( 429/ 354- 483/ 453)  60.3 E62-E54 S21-S14
     11 1330  391 (S41W03 L=211)   434    101        0 (1293/ 368-1375/ 420)   9.5 W00-W06 S42-S39
     12  698  314 (S29E45 L=164)   343     58        2 ( 676/ 288- 721/ 336)  76.0 E47-E41 S30-S26
     13 1092  239 (S46E22 L=187)   351     54        0 (1062/ 225-1125/ 253) -11.3 E24-E19 S47-S44
  

    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