Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  603 1713 (N48E43 L=345)   558     94        2 ( 576/1683- 624/1744) -45.0 E45-E39 N45-N50
      2  840 1498 (N33E15 L= 13)   725    168        0 ( 771/1444- 907/1542) -26.6 E21-E09 N29-N35
      3  437 1390 (N23E44 L=344)   780    134        1 ( 412/1331- 458/1467)  90.0 E45-E41 N19-N28
      4 1381 1323 (N22W23 L= 51)  1729    364        2 (1217/1232-1544/1420) -26.6 W12-W34 N16-N28
      5  907 1313 (N21E09 L= 20)   909    182        3 ( 860/1222- 955/1399) -65.6 E12-E05 N14-N26
      6  717 1087 (N06E19 L=  9)  1529    402        2 ( 629/ 984- 865/1239) -38.7 E26-E09 N00-N14
      7  570  943 (S04E29 L=359)   274     79        1 ( 530/ 924- 605/ 960)  11.3 E31-E26 S04-S02
      8 1304  256 (S52W32 L= 60)   547    161        0 (1260/ 228-1349/ 281) -19.7 W25-W39 S54-S48
  

    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