Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1666 1113 (N16W41 L=284)  5642    903       10 (1368/ 785-1858/1428) -29.5 W15-W68 N00-N30
      2  312  923 (S19E47 L=196)   322     85        1 ( 293/ 882- 338/ 949) -52.1 E48-E44 S22-S17
      3 1726  698 (S07W56 L=299)  1191    257        2 (1655/ 669-1811/ 726)  10.6 W48-W65 S10-S04
      4  539  653 (S33E27 L=216)   360    112        2 ( 493/ 633- 599/ 686)  -6.3 E30-E21 S33-S31
      5  907  656 (S25E01 L=243)   469    121        1 ( 844/ 631- 964/ 680)   0.0 E04-W03 S26-S23
      6 1441  583 (S20W36 L=279)   590    153        0 (1372/ 564-1506/ 606)   5.7 W30-W40 S20-S18
      7 1382  349 (S36W43 L=286)   366     44        1 (1362/ 326-1401/ 378) -53.1 W41-W45 S38-S34
      8 1287  309 (S41W39 L=282)   338     50        1 (1261/ 293-1310/ 331)   0.0 W36-W41 S42-S39
  

    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