Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1573 1326 (N15W42 L=247)   451     82        3 (1548/1295-1607/1360)  56.3 W39-W44 N12-N16
      2 1112 1315 (N22W10 L=216)   390    126        4 (1053/1291-1171/1340) -33.7 W07-W14 N20-N24
      3  803 1226 (N22E11 L=194)   547    193        1 ( 767/1130- 853/1300)  74.7 E13-E06 N15-N26
      4 1671 1230 (N07W47 L=252)   329     77        1 (1639/1197-1700/1260) -51.3 W44-W48 N04-N09
      5 1561 1072 (N00W35 L=240)   352     79        0 (1524/1046-1593/1100) -36.9 W32-W36 S01-N02
      6  878  988 (N06E09 L=196)   258     60        1 ( 857/ 948- 898/1016) -90.0 E10-E08 N03-N08
      7 1377  816 (S12W18 L=224)   299    146        0 (1364/ 744-1394/ 887)  83.7 W17-W19 S16-S07
      8  189  864 (N06E66 L=140)   537    102        2 ( 172/ 837- 214/ 886)  15.9 E68-E61 N03-N06
      9  467  607 (S12E44 L=161)   283     37        0 ( 447/ 588- 484/ 628)  45.0 E46-E42 S13-S10
  

    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