Objective. Temperature plays a crucial role in influencing the spatiotemporal dynamics of the heart. Electrical instabilities due to specific thermal conditions typically lead to early period-doubling bifurcations and beat-to-beat alternans. These pro-arrhythmic phenomena manifest in voltage and calcium traces, resulting in compromised contractile behaviors. In such intricate scenario, dual optical mapping technique was used to uncover unexplored multi-scale and nonlinear couplings, essential for early detection and understanding of cardiac arrhythmia. Approach. We propose a methodological analysis of synchronized voltage-calcium signals for detecting alternans, restitution curves, and spatiotemporal alternans patterns under different thermal conditions, based on integral features calculation. To validate our approach, we conducted a cross-species investigation involving rabbit and guinea pig epicardial ventricular surfaces and human endocardial tissue under pacing-down protocols. Main results. We show that the proposed integral feature, as the area under the curve, could be an easily applicable indicator that may enhance the predictability of the onset and progression of cardiac alternans. Insights into spatiotemporal correlation analysis of characteristic spatial lengths across different heart species were further provided. Significance. Exploring cross-species thermoelectric features contributes to understanding temperature-dependent proarrhythmic regimes and their implications on coupled spatiotemporal voltage-calcium dynamics. The findings provide preliminary insights and potential strategies for enhancing arrhythmia detection and treatment.