A lot of researches have been conducted to estimate the volatility smile effect shown in the option market. This paper proposes a method to approximate an implied volatility function, given noisy real market option data. To construct an implied volatility function, we use Gaussian Processes (GPs). Their output values are implied volatilities while moneyness values (the ratios of strike price to underlying asset price) and time to maturities are as their input values. To show the performances of our proposed method, we conduct experimental simulations with Korean Equity-Linked Warrant (ELW) market data as well as toy data.