Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  425 1491 (N31E50 L=179)   358     59        1 ( 403/1467- 454/1519) -23.2 E53-E46 N29-N32
      2  529 1190 (N12E33 L=196)   233     86        0 ( 512/1148- 540/1233) -68.2 E35-E32 N09-N14
      3 1766  951 (S01W53 L=282)   373     64        0 (1746/ 927-1793/ 982) -26.6 W50-W56 S02-N00
      4  428  707 (S19E42 L=188)   339     85        0 ( 387/ 692- 469/ 721)   8.1 E45-E38 S20-S17
      5 1699  652 (S20W52 L=281)   384     42        1 (1681/ 629-1720/ 676)  71.6 W50-W53 S21-S18
      6  426  613 (S25E44 L=185)   515    137        0 ( 411/ 550- 462/ 666)  53.1 E46-E39 S29-S21
      7 1381  584 (S24W26 L=255)   282    103        0 (1349/ 540-1412/ 628)  63.4 W24-W27 S27-S21
      8  799  517 (S30E15 L=214)   268     61        1 ( 762/ 502- 830/ 536)  11.3 E17-E12 S31-S28
      9 1513  456 (S34W41 L=270)   329     16        1 (1493/ 438-1533/ 477)  18.4 W39-W43 S35-S32
     10 1135  368 (S41W11 L=240)   514     45        1 (1105/ 343-1169/ 392)  38.7 W08-W13 S43-S39
     11  883  327 (S46E10 L=219)   887    187        1 ( 802/ 284- 941/ 363) -28.6 E17-E05 S49-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