Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1288 1777 (N41W52 L=169)   412     71        0 (1272/1737-1303/1813)  21.8 W47-W57 N38-N43
      2  175 1133 (N32E65 L= 53)  2638    383        7 ( 122/ 997- 209/1327)  68.2 E74-E59 N23-N44
      3  368 1042 (N23E42 L= 75)  1210    316        3 ( 254/ 931- 453/1121) -43.0 E53-E36 N15-N29
      4  203  922 (N19E62 L= 56)  1361    289        3 ( 174/ 789- 234/1037)  66.8 E64-E56 N11-N25
      5  259  952 (N20E54 L= 63)   625    116        0 ( 244/ 890- 278/1010)  81.9 E56-E52 N16-N23
      6 1759  934 (S20W45 L=162)   700    160        0 (1722/ 873-1801/1009)  33.7 W40-W50 S22-S16
      7 1494  893 (S14W24 L=141)   899    251        1 (1382/ 845-1595/ 948)  -7.1 W15-W31 S15-S12
      8 1727  622 (S40W40 L=157)  1497    295        2 (1620/ 541-1827/ 726)   9.1 W27-W52 S42-S35
      9  866  393 (S26E29 L= 88)  3256    903        2 ( 692/ 252-1151/ 628)  13.1 E46-E03 S34-S18
     10 1214  512 (S28E03 L=114)   616    240        0 (1144/ 442-1248/ 623)  38.7 E09-W01 S31-S22
  

    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