Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1045 1605 (N36W11 L=151)   581    121        0 ( 983/1577-1104/1624) -26.6 W06-W14 N32-N37
      2  864 1549 (N35E03 L=137)   602    168        1 ( 785/1498- 939/1611)  22.6 E09-W03 N32-N38
      3 1827  943 (S15W56 L=196)   650    157        1 (1793/ 881-1858/1024)  58.0 W52-W59 S19-S10
      4 1568  919 (S13W33 L=173)   312     71        0 (1528/ 892-1591/ 951)  31.0 W29-W34 S14-S11
      5 1556  791 (S21W31 L=171)   978    260        1 (1462/ 714-1648/ 894)  32.7 W24-W38 S24-S15
      6  572  689 (S14E33 L=106)   511     98        2 ( 525/ 648- 607/ 735)  23.2 E37-E30 S16-S11
      7 1042  642 (S23E04 L=136)   629    275        0 ( 988/ 532-1133/ 730)  45.0 E09-W02 S29-S18
      8  276  600 (S15E63 L= 77)  2270    170        3 ( 229/ 532- 327/ 661)  76.0 E69-E58 S19-S11
      9  272  519 (S20E69 L= 70)   801     99        0 ( 241/ 496- 300/ 546) -38.7 E73-E64 S21-S17
     10  744  437 (S32E30 L=110)   375    117        0 ( 703/ 401- 782/ 480)  29.1 E34-E25 S34-S29
     11  631  286 (S41E49 L= 91)   497     68        2 ( 603/ 258- 658/ 316)  51.3 E51-E46 S43-S38
  

    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