Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  848 1621 (N49W03 L= 99)   595     72        1 ( 824/1584- 887/1656)  68.2 W00-W05 N46-N51
      2  871 1346 (N29E03 L= 93)   463    189        3 ( 825/1279- 939/1398)  24.0 E07-W02 N26-N31
      3  422 1163 (N26E39 L= 57)  1242    377        2 ( 300/1087- 531/1222) -25.3 E50-E30 N19-N30
      4  237 1064 (N22E58 L= 38)   828    175        1 ( 197/ 992- 283/1128)  35.0 E62-E52 N18-N25
      5  725  665 (S09E26 L= 70)   237     44        1 ( 699/ 646- 749/ 683) -56.3 E27-E24 S10-S07
      6 1069  482 (S27E09 L= 87)   288     72        0 (1037/ 453-1097/ 514)  31.0 E12-E06 S28-S25
      7  523  433 (S20E50 L= 46)   462     71        2 ( 501/ 406- 555/ 457)   9.5 E52-E46 S22-S18
      8  667  359 (S28E42 L= 54)  1321    238        1 ( 605/ 269- 735/ 449)  79.4 E45-E39 S36-S20
      9 1304  213 (S59W00 L= 96)   522     62        0 (1269/ 196-1338/ 229) -45.0 E03-W04 S63-S56
  

    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