Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  764 1767 (N31E45 L=256)   822    200        0 ( 723/1679- 815/1837) -60.3 E49-E40 N23-N37
      2 1008 1695 (N33E24 L=276)   521    123        1 ( 990/1629-1026/1749) -45.0 E27-E19 N28-N37
      3 1069 1352 (N13E07 L=294)  1075    334        3 (1044/1210-1093/1515) -57.1 E13-E02 N05-N22
      4  326 1122 (S18E48 L=252)   474    197        0 ( 277/1082- 382/1191)  41.2 E52-E44 S20-S13
      5 1558  523 (S16W51 L=351)   371     69        0 (1533/ 496-1579/ 547)  -9.5 W47-W53 S18-S15
      6 1647  472 (S16W62 L=  3)   807    111        3 (1597/ 451-1699/ 493)  29.1 W57-W67 S17-S12
      7 1382  442 (S27W44 L=344)  1198    257        1 (1273/ 418-1504/ 480)  37.6 W36-W50 S30-S20
      8 1574  394 (S22W63 L=  3)  1457    127        5 (1522/ 353-1623/ 428)  26.6 W56-W69 S25-S19
      9  922  257 (S56W29 L=329)   621     83        2 ( 897/ 220- 945/ 293) -24.4 W23-W34 S59-S53
     10  747  240 (S67W17 L=317)  1210    158        4 ( 704/ 190- 790/ 288) -38.7 W09-W28 S72-S60
  

    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