目录
任务简介
任务解决思路与经验收获
具体步骤
总结
我是政胤 期待你的关注
大家好 我是政胤 今天教大家一个比较刑的爬虫案例
Python爬虫 爬取下载美国科研网站数据
制作不易 点个免费的关注 支持一下吧
任务简介
目标网站:https://app.powerbigov.us/view” />
目标网站是漂亮国的科研数据,是PowerBI实现的网页数据,无法使用Ctrl+C复制内容,因此,求助于我们进行爬取。
任务解决思路与经验收获
首先任务可拆解为两个部分:一是从网站爬取数据到本地,二是解析数据输出CSV文件:
爬取数据部分:
解析网页,找到数据异步加载的实际请求地址与参数
书写爬虫代码获取全部数据
解析数据部分
这是本次任务的主要难点所在,难点在于:在返回的数据list中,元素不是固定的个数,只有与上一行不同的数值,而具体哪一列不同、哪一列相同,是使用一个“R”值表示,正常解决思路是要通过JS逆向,找出解析R关系的函数,完成解析。但是,由于网页的JS非常复杂,且许多函数名都是简写,阅读十分困难,一直没有逆向成功
在解决该问题上,先是手工查询总结关系,完成了第一个版本,没想到后续在写写这篇分享文章时突然思路打开,改变了请求数据方式,绕过了分析R关系的步骤:
方案一:按正常请求,使用R关系解析数据
下载完整的数据后分析,所需要的2009至2013年的数据中,R关系一共有124种,最小值0、最大值4083,通过人工查询这124种关系,制作成字典,完成解析。总结出的关系如下图(手工查询了5个小时,累啊):
方案二:以时间换空间,每次仅请求一行数据,绕过解析R关系的难题
在复盘时,突然头脑开窍,请求到的数据第一行一定是完整的,要是每次只请求一行数据,那就不可能存在与上一行相同的情形了,这种情况下就能绕过解析R关系这一难题。测试后方案可行,只是需要考虑以下问题:开启多线程加速以缩时间,但即使开启多线程,也只能按12个年份开启12个线程,而行数最多的年份约2万行,爬虫需要运行约5至6个小时
断点续爬,避免程序异常中断后,需要从头开始;
具体步骤
目标网站分析
第一步当然是对目标网站进行分析,找到数据正确的请求地址,这点很容易,打开Chrome的开发者模式,向下拖动滚动条,看到新出现的请求,就是真实的地址了,直接上结果:
然后看一下POST请求的参数
请求的参数
再看一下Headers,意外发现,竟然没有反爬!没有反爬!没有反爬!好吧,漂亮国的网站就是大气。
- 分析小结:
# 完整参数就略过,关键参数以下三项:# 1.筛选年份的参数,示例:param['queries'][0]['Query']['Commands'][0]['SemanticQueryDataShapeCommand']['Query']['Where'][0]['Condition']['In']['Values'][0][0]['Literal']['Value'] = '2009L'# 2.请求下一批数据(请求首批数据时无需传入该参数),示例:param['queries'][0]['Query']['Commands'][0]['SemanticQueryDataShapeCommand']['Binding']['DataReduction']['Primary']['Window']['RestartTokens'] = [["'Augusto E Caballero-Robles'","'Physician'","'159984'","'Daiichi Sankyo, Inc.'","'CC0131'","'Basking Ridge'","'NJ'","'Compensation for Bona Fide Services'","2009L","'4753'"]]# 注:以上"RestartTokens"的值在前一批数据的response中,为上一批数据的返回字典值,示例res['results'][0]['result']['data']['dsr']['DS'][0]['RT']# 3.请求页面的行数(浏览器访问默认是500行/页,但爬虫访问的话...你懂的),示例:param['queries'][0]['Query']['Commands'][0]['SemanticQueryDataShapeCommand']['Binding']['DataReduction']['Primary']['Window']['Count'] = 500# 参数还有很多,例如排序的参数ordby,各种筛选项等
请求数据的网址:
https://wabi-us-gov-virginia-api.analysis.usgovcloudapi.net/public/reports/querydata" />
在获取全部数据上,则使用一个while True死循环,每次请求返回值中如有"RT"关键字,则修改POST参数,发起下一个请求,直至返回值中没有"RT"关键字,代表全部数据爬取结束(详见代码)
请求过程中出现的异常需要捕获,根据异常类型决定下一步操作,由于该网站没有反爬,只有超时或连接错误的异常,因此,只需要重启发起请求即可,因此可以不考虑断点续爬
以上步骤,详细代码如下:
网站没有反爬,找到正确的路径和参数后,在爬虫代码实现上相对简单,直接发起post请求即可,代码中通过PageSpider类实现(详细代码附后)
在断点续传上,通过流程解决,把每行数据存储到TXT文件中,文件名记录年份以及行数,先读取已爬取的记录,找到最后一次请求结果,然后发起后续请求。
"""爬取页面数据的爬虫"""import pathlib as plimport requestsimport jsonimport timeimport threadingimport urllib3def get_cost_time(start: time.time, end: time.time = None):"""计算间隔时长的方法:param start: 起始时间:param end: 结束时间,默认为空,按最新时间计算:return: 时分秒格式"""if not end:end = time.time()cost = end - startdays = int(cost / 86400)hours = int(cost % 86400 / 3600)mins = int(cost % 3600 / 60)secs = round(cost % 60, 4)text = ''if days:text = f'{text}{days}天'if hours:text = f'{text}{hours}小时'if mins:text = f'{text}{mins}分钟'if secs:text = f'{text}{secs}秒'return textclass PageSpider:def __init__(self, year: int, nrows: int = 500, timeout: int = 30):"""初始化爬虫的参数:param year: 下载数据的年份,默认空,不筛选年份,取得全量数据:param nrows: 每次请求获取的数据行数,默认500,最大30000(服务器自动限制,超过无效):param timeout: 超时等待时长"""self.year = year if year else 'all'self.timeout = timeout# 请求数据的地址self.url = 'https://wabi-us-gov-virginia-api.analysis.usgovcloudapi.net/public/reports/querydata" /> dict:"""制作请求体中的参数:param key: 下一页的关键字RestartTokens,默认空,第一次请求时无需传入该参数:return: dict"""params = self.params.copy()if key:params['queries'][0]['Query']['Commands'][0]['SemanticQueryDataShapeCommand']['Binding']['DataReduction']['Primary']['Window']['RestartTokens'] = keyreturn paramsdef crawl_pages(self, idx: int = 1, key: list = None):"""爬取页面并输出TXT文件的方法,:param idx: 爬取的索引值,默认为1,在每行爬取时,代表行数:param key: 下一页的关键字RestartTokens,默认空,第一次请求时无需传入该参数:return: None"""start = time.time()while True:# 创建死循环爬取直至结束try:res = requests.post(url=self.url, headers=self.headers, json=self.make_params(key),timeout=self.timeout)except (requests.exceptions.ConnectTimeout,requests.exceptions.ConnectionError,urllib3.exceptions.ConnectionError,urllib3.exceptions.ConnectTimeoutError):# 捕获超时异常 或 连接异常print(f'{self.year}_part{idx}: timeout, wait 5 seconds retry')time.sleep(5)# 休息5秒后再次请求continue# 跳过后续步骤except Exception as e:# 其他异常,打印一下异常信息print(f'{self.year}_part{idx} Error: {e}')time.sleep(5)# 休息5秒后再次请求continue# 跳过后续步骤if res.status_code == 200:with open(f'./tmp/{self.year}_part{idx}.txt', 'w', encoding='utf-8') as fout:fout.write(res.text)if idx % 100 == 0:print(f'{self.year}的第{idx}行数据写入完成,已用时: {get_cost_time(start)}')key = json.loads(res.text)['results'][0]['result']['data']['dsr']['DS'][0].get('RT', None)if not key:# 如果没有RT值,说明已经全部爬取完毕了,打印一下信息退出print(f'{self.year} completed max_idx is {idx}')returnidx += 1else:# 打印一下信息重新请求print(f'{self.year}_part{idx} not 200,check please', res.text)continuedef mul_crawl(year: int, nrows: int = 2):"""多线程爬取的方法,注按行爬取:param year: 需要爬取的年份:param nrows: 每份爬取的行数,若每次仅爬取1行数据,nrows参数需要为2,才会有下一行,否则都是第一行"""# 定义爬虫对象spider = PageSpider(year, nrows=nrows)# 获取爬取对象已爬取的idx,key和是否完成爬取的信号singleidx, key, single = spider.get_idx_and_rt()if not single:print(f'{year}年的共{idx}行数据已经全部下载,无需爬取')returnprint(f'{year}年的爬虫任务启动, 从{idx+1}行开始爬取')spider.crawl_pages(idx+1, key)# 特别注意,已经爬取了idx行,重启时,下一行需要+1,否则重启后,会覆盖一行数据if __name__ == '__main__':pools = []for y in range(2009, 2021):pool = threading.Thread(target=mul_crawl, args=(y, 2), name=f'{y}_thread'# 按行爬取,nrows参数需要为2)pool.start()pools.append(pool)for pool in pools:pool.join()print('任务全部完成')
解析数据
方案一
代码运行示例:
以时间换空间,每次仅请求一行,绕过R关系解析
解析数据困难的部分就是找出R关系规律,这部分是使用手工查询来解决的,直接上代码吧:
class ParseData: """ 解析数据的对象 """ def __init__(self, file_path: pl.Path = None): """ 初始化对象 :param file_path: TXT数据存放的路径,默认自身目录下的tmp文件夹 """ self.file_path = pl.Path('./tmp') if not file_path else file_path self.files = list(self.file_path.glob('2*.txt')) self.cols_dict = None self.colname_dict = { 'D0': 'License Type', 'D1': 'License Number', 'D2': 'Manufacturer Full Name', 'D3': 'Manufacturer ID', 'D4': 'City', 'D5': 'State', 'D6': 'Full Name', 'D7': 'Payment Category', 'D8': 'Covered Recipient ID' } self.colname_dict_T = {v: k for k, v in self.colname_dict.items()} def make_excels(self): """ 将每个数据文件单独转换为excel数据表用于分析每份数据 :return: """ for file in self.files: with open(file, 'r') as fin: res = json.loads(fin.read()) dfx = pd.DataFrame(res['results'][0]['result']['data']['dsr']['DS'][0]['PH'][0]['DM0']) dfx['filename'] = file.stem dfx[['year', 'part']] = dfx['filename'].str.split('_', expand=True) dfx['C_count'] = dfx['C'].map(len) writer = pd.ExcelWriter(self.file_path / f'{file.stem}.xlsx') dfx.to_excel(writer, sheet_name='data') for k, v in res['results'][0]['result']['data']['dsr']['DS'][0]['ValueDicts'].items(): dfx = pd.Series(v).to_frame() dfx.to_excel(writer, sheet_name=k) writer.save() print('所有数据均已转为Excel') def make_single_excel(self): """ 将所有数据生成一份excel文件,不包含字典 :return: """ # 合并成整个文件 df = pd.DataFrame() for file in self.files: with open(file, 'r') as fin: res = json.loads(fin.read()) dfx = pd.DataFrame(res['results'][0]['result']['data']['dsr']['DS'][0]['PH'][0]['DM0']) dfx['filename'] = file.stem dfx[['year', 'part']] = dfx['filename'].str.split('_', expand=True) dfx['C_count'] = dfx['C'].map(len) df = pd.concat([df, dfx]) return df def get_cols_dict(self): """ 读取列关系的字典 :return: """ # 读取列字典表 self.cols_dict = pd.read_excel(self.file_path.parent / 'cols_dict.xlsx') self.cols_dict.set_index('R', inplace=True) self.cols_dict = self.cols_dict.dropna() self.cols_dict.drop(columns=['C_count', ], inplace=True) self.cols_dict.columns = [col.split(':')[-1] for col in self.cols_dict.columns] self.cols_dict = self.cols_dict.astype('int') def make_dataframe(self, filename): """ 读取TXT文件,转换成dataframe :param filename: 需要转换的文件 :return: """ with open(filename, 'r') as fin: res = json.loads(fin.read()) df0 = pd.DataFrame(res['results'][0]['result']['data']['dsr']['DS'][0]['PH'][0]['DM0']) df0['R'] = df0['R'].fillna(0) df0['R'] = df0['R'].map(int) values_dict = res['results'][0]['result']['data']['dsr']['DS'][0]['ValueDicts'] dfx = [] for idx in df0.index: row_value = df0.loc[idx, 'C'].copy() cols = self.cols_dict.loc[int(df0.loc[idx, 'R'])].to_dict() row = {} for col in ['License Type', 'License Number', 'Manufacturer Full Name', 'Manufacturer ID', 'City', 'State', 'Full Name', 'Payment Category', 'Disclosure Year', 'Covered Recipient ID', 'Amount of Payment', 'Number of Events Reflected']: v = cols.get(col) if v: value = row_value.pop(0) if col in self.colname_dict.values(): if not isinstance(value, str): value_list = values_dict.get(self.colname_dict_T.get(col), []) value = value_list[value] row[col] = value else: row[col] = None row['R'] = int(df0.loc[idx, 'R']) dfx.append(row) dfx = pd.DataFrame(dfx) dfx = dfx.fillna(method='ffill') dfx[['Disclosure Year', 'Number of Events Reflected']] = dfx[ ['Disclosure Year', 'Number of Events Reflected']].astype('int') dfx = dfx[['Covered Recipient ID', 'Full Name', 'License Type', 'License Number', 'Manufacturer ID', 'Manufacturer Full Name', 'City', 'State', 'Payment Category', 'Amount of Payment', 'Number of Events Reflected', 'Disclosure Year', 'R']] return dfx def parse_data(self, out_name: str = None): """ 解析合并数据 :param out_name: 输出的文件名 :return: """ df = pd.DataFrame() for n, f in enumerate(self.files): dfx = self.make_dataframe(f) df = pd.concat([df, dfx]) print(f'完成第{n + 1}个文件,剩余{len(self.files) - n - 1}个,共{len(self.files)}个') df.drop(columns='R').to_csv(self.file_path / f'{out_name}.csv', index=False) return df
方案二
使用方案二处理数据时,在进行数据后验后发现,还有两个细节问题需要解决:
一是返回值中出现了新的关键字“Ø”,经手工验证才知道代表输出的行中,存在本身就是空值的情况,遍历数据后,发现只有出现3个不同值(60, 128, 2048),因此,手工制作了col_dict(详见代码)。\class ParseDatav2: """ 解析数据的对象第二版,将按行爬取的的json文件,转换成dataframe,增量写入csv文件, 因每次请求一行,首行数据不存在与上一行相同情形,因此,除个别本身无数据情况,绝大多数均为完整的12列数据, """ def __init__(self): """ 初始化 """ # 初始化一行的dataframe, self.row = pd.DataFrame([ 'Covered Recipient ID', 'Full Name', 'License Type', 'License Number', 'Manufacturer ID', 'Manufacturer Full Name', 'City', 'State', 'Payment Category', 'Amount of Payment', 'Number of Events Reflected', 'Disclosure Year' ]).set_index(0) self.row[0] = None self.row = self.row.T self.row['idx'] = None # 根据 Ø 值的不同选择不同的列,目前仅三种不同的Ø值,注0为默认值,指包含所有列 self.col_dict = { # 完整的12列 0: ['License Type', 'License Number', 'Manufacturer Full Name', 'Manufacturer ID', 'City', 'State', 'Full Name', 'Payment Category', 'Disclosure Year', 'Covered Recipient ID', 'Amount of Payment', 'Number of Events Reflected'], # 有4列是空值,分别是 'Manufacturer Full Name', 'Manufacturer ID', 'City', 'State' 60: ['License Type', 'License Number', 'Full Name', 'Payment Category', 'Disclosure Year', 'Covered Recipient ID', 'Amount of Payment', 'Number of Events Reflected'], # 有1列是空值,是 'Payment Category' 128: ['License Type', 'License Number', 'Manufacturer Full Name', 'Manufacturer ID', 'City', 'State', 'Full Name', 'Disclosure Year', 'Covered Recipient ID', 'Amount of Payment', 'Number of Events Reflected'], # 有1列是空值,是 'Number of Events Reflected' 2048: ['License Type', 'License Number', 'Manufacturer Full Name', 'Manufacturer ID', 'City', 'State', 'Full Name', 'Payment Category', 'Disclosure Year', 'Covered Recipient ID', 'Amount of Payment'], } # 列名转换字典 self.colname_dict = { 'License Type': 'D0', 'License Number': 'D1', 'Manufacturer Full Name': 'D2', 'Manufacturer ID': 'D3', 'City': 'D4', 'State': 'D5', 'Full Name': 'D6', 'Payment Category': 'D7', 'Covered Recipient ID': 'D8' } # 储存爬取的json文件的路径 self.data_path = pl.Path('./tmp') # 获取json文件的迭代器 self.files = self.data_path.glob('*.txt') # 初始化输出文件的名称及路径 self.file_name = self.data_path.parent / 'data.csv' def create_csv(self): """ 先输出一个CSV文件头用于增量写入数据 :return: """ self.row.drop(0, axis=0).to_csv(self.file_name, index=False) def parse_data(self, filename: pl.Path): """ 读取按1行数据请求获取的json文件,一行数据 :param filename: json文件的路径 :return: None """ row = self.row.copy() # 复制一行dataframe用于后续修改 res = PageSpider.read_json(filename) # 获取数据中的valuedicts valuedicts = res['results'][0]['result']['data']['dsr']['DS'][0]['ValueDicts'] # 获取数据中每行的数据 row_values = res['results'][0]['result']['data']['dsr']['DS'][0]['PH'][0]['DM0'][0]['C'] # 获取数据中的'Ø'值(若有),该值代表输出的行中,存在空白部分,用于确定数据列 cols = ic(self.col_dict.get( res['results'][0]['result']['data']['dsr']['DS'][0]['PH'][0]['DM0'][0].get('Ø', 0) )) # 遍历每行数据,修改row这个dataframe的值 for col, value in zip(cols, row_values): ic(col, value) colname = self.colname_dict.get(col) # colname转换,D0~D8 if colname: # 如果非空,则需要转换值 value = valuedicts.get(self.colname_dict.get(col))[0] # 修改dataframe数据 row.loc[0, col] = value # 写入索引值 row['idx'] = int(filename.stem.split('_')[-1].replace('part', '')) return row def run(self): """ 运行写入程序 """ self.create_csv() for idx, filename in enumerate(self.files): row = self.parse_data(filename) row.to_csv(self.file_name, mode='a', header=None, index=False) print(f'第{idx + 1}个文件{filename.stem}写入表格成功') print('全部文件写入完成')
二是每行数据请求,nrows需要设置为2,而最后一行数据无法通过该方式获取,因此,需要从最后一个返回的json数据中解析出最后一行数据(详见LastRow类)
class LastRow: """ 获取并写入最后一行数据的类 由于每次请求一行数据的方式,存在缺陷,无法获取到最后一行数据, 本方法是对最后一个能够获取的json(倒数第二行)进行解析,取得最后一行数据, 本方法存在缺陷,即默认最后一行“Amount of Payment”列值一定与倒数第二行不同, 目前2009年至2020年共12年的数据中,均满足上述条件,没有出错。 除本方法外,还可以通过逆转排序请求的方式,获取最后一行数据 """ def __init__(self): """ 初始化 """ self.file_path = pl.Path('./tmp') # 存储爬取json数据的路径 self.files_df = pd.DataFrame() # 初始化最后一份请求的dataframe # 列名对应的字典 self.colname_dict = { 'D0': 'License Type', 'D1': 'License Number', 'D2': 'Manufacturer Full Name', 'D3': 'Manufacturer ID', 'D4': 'City', 'D5': 'State', 'D6': 'Full Name', 'D7': 'Payment Category', 'year': 'Disclosure Year', 'D8': 'Covered Recipient ID', 'M0': 'Amount of Payment', 'M1': 'Number of Events Reflected' } self.data = pd.DataFrame() # 初始化最后一行数据data def get_last_file(self): """ 遍历文件夹,取得最后一份请求的dataframe """ self.files_df = pd.DataFrame(list(self.file_path.glob('*.txt')), columns=['filename']) self.files_df[['year', 'idx']] = self.files_df['filename'].map(lambda x: x.stem).str.split('_', expand=True) self.files_df['idx'] = self.files_df['idx'].str.replace('part', '') self.files_df['idx'] = self.files_df['idx'].astype(int) self.files_df.sort_values(by=['year', 'idx'], inplace=True) self.files_df = self.files_df.drop_duplicates('year', keep='last') def get_last_row(self, ser: pd.Series) -> pd.DataFrame: """ 解析文件,获取最后一行的数据 :param ser: 一行文件信息的series """ # 读取json数据 res = PageSpider.read_json(ser['filename']) # 获取values_dict values_dict = res['results'][0]['result']['data']['dsr']['DS'][0]['ValueDicts'] # 获取文件中的第一行数据 row_values = res['results'][0]['result']['data']['dsr']['DS'][0]['PH'][0]['DM0'][0]['C'] # 获取文件中的下一行数据,因文件是倒数第二行的数据,因此下一行即为最后一行 next_row_values = res['results'][0]['result']['data']['dsr']['DS'][0]['PH'][0]['DM0'][1]['C'] # 初始化Series row = pd.Series() # 解析数据填充series for k, col in self.colname_dict.items(): value = row_values.pop(0) if k.startswith('D'): # 如果K值是D开头 values = values_dict[k] if len(values) == 2: value = next_row_values.pop(0) value = values[-1] elif k == 'year': pass else: if next_row_values: value = next_row_values.pop(0) row[col] = value row['idx'] = ser['idx'] + 1 row = row.to_frame().T return row def run(self): """ 运行获取最后一行数据的方法 """ self.get_last_file() for i in self.files_df.index: self.data = pd.concat([self.data, self.get_last_row(self.files_df.loc[i])]) self.data = self.data[[ 'Covered Recipient ID', 'Full Name', 'License Type', 'License Number', 'Manufacturer ID', 'Manufacturer Full Name', 'City', 'State', 'Payment Category', 'Amount of Payment', 'Number of Events Reflected', 'Disclosure Year', 'idx' ]] filename = self.file_path.parent / 'data.csv' self.data.to_csv(filename, mode='a', index=False, header=None) return self.data
结果展示
2010年
2015年
2020年尾部
延伸思考
如果将上述方案一与方案二结合,整理出所有不同R关系的行样例,使用方案二爬取少量的部分示例,然后推导出完整的R关系字典,再使用方案一的方法进行爬取解析,将大大节约时间。该方式在数据量远远超过当前数量时,可以考虑使用。
总结
完成整个项目过程中历经了:暗爽(不到1小时就完成了爬虫部分功能)->迷茫(JS逆向失败,无法总结R关系规律)->焦虑与烦躁(担心无法完成任务,手工查询规则5个多小时)->开窍(复盘过程中突然发现新思路)一系列过程。最终结果还是较为顺利的完成了整个任务,而最大的感触还是思路的开拓:一条路走不通时,也许换个方向就能解决问题(注:count参数500一开始就使用了,只是一直在增加请求的行数,而一直没有想到减少请求的行数这么一个小小的改变,就能带来巨大的突破)。
我是政胤 期待你的关注