Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1127 1783 (N50E16 L= 50)  3787    333        3 (1045/1644-1208/1913) -32.5 E36-E03 N39-N62
      2  381 1530 (N14E58 L=  7)  1870    246        5 ( 331/1446- 466/1625) -56.3 E66-E50 N08-N20
      3  593 1225 (N00E30 L= 36)  1209    248        3 ( 528/1128- 673/1327)  69.0 E32-E27 S06-N07
      4  796 1001 (S08E13 L= 53)   496    152        1 ( 738/ 954- 865/1040)  53.1 E15-E09 S12-S04
      5 1766  791 (N01W55 L=121)  3141    300        3 (1701/ 669-1820/ 937) -37.4 W47-W64 S05-N08
      6 1421  866 (S03W27 L= 92)   592    109        0 (1382/ 823-1461/ 913)  68.2 W25-W28 S05-N00
      7  425  689 (S35E35 L= 30)   307    126        1 ( 383/ 630- 450/ 716) -33.7 E39-E33 S39-S33
      8 1741  641 (S08W59 L=125)   642     79        1 (1709/ 602-1769/ 671) -18.4 W55-W64 S09-S06
      9 1172  175 (S54W50 L=116)   462     55        0 (1153/ 144-1193/ 205) -31.0 W46-W56 S56-S50
  

    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