Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  641 1605 (N51E11 L=265)   411     72        0 ( 626/1568- 655/1649)  81.9 E12-E08 N47-N54
      2 1305 1586 (N28W36 L=311)   358     74        0 (1261/1572-1347/1598) -18.4 W32-W38 N27-N29
      3  395 1160 (N27E38 L=237)   479     92        0 ( 356/1119- 425/1199) -56.3 E41-E35 N24-N30
      4  475 1074 (N21E33 L=243)   299     85        0 ( 433/1055- 518/1089)   0.0 E36-E29 N19-N21
      5 1853  791 (S33W56 L=331)  1281    118        1 (1817/ 740-1881/ 849)  45.0 W51-W60 S36-S29
      6 1185  545 (S27E03 L=272)   334    129        0 (1134/ 502-1235/ 583)  15.9 E07-W00 S28-S25
      7 1522  313 (S59W20 L=295)   454     43        0 (1490/ 302-1550/ 327) -24.0 W15-W24 S61-S57
      8  718  306 (S30E44 L=231)   312     51        0 ( 687/ 293- 752/ 320) -11.3 E47-E41 S31-S29
      9 1325  219 (S60E07 L=268)   663     83        2 (1288/ 198-1370/ 243) -49.4 E11-E01 S63-S55
  

    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