Shallow Networks for High-Accuracy Road Object-Detection

06/05/2016 ∙ by Khalid Ashraf, et al. ∙ berkeley college 0

The ability to automatically detect other vehicles on the road is vital to the safety of partially-autonomous and fully-autonomous vehicles. Most of the high-accuracy techniques for this task are based on R-CNN or one of its faster variants. In the research community, much emphasis has been applied to using 3D vision or complex R-CNN variants to achieve higher accuracy. However, are there more straightforward modifications that could deliver higher accuracy? Yes. We show that increasing input image resolution (i.e. upsampling) offers up to 12 percentage-points higher accuracy compared to an off-the-shelf baseline. We also find situations where earlier/shallower layers of CNN provide higher accuracy than later/deeper layers. We further show that shallow models and upsampled images yield competitive accuracy. Our findings contrast with the current trend towards deeper and larger models to achieve high accuracy in domain specific detection tasks.

READ FULL TEXT VIEW PDF
POST COMMENT

Comments

There are no comments yet.

Authors

page 3

page 7

page 8

This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

1 Introduction and Motivation

Advanced driver assistance systems (ADAS) and increasingly autonomous vehicles promise to make transportation safe, efficient and cost effective. Driven by the goal of building a safe transportation system, ADAS has emerged as a leading research direction in recent years  [13, 14, 18, 1]

. Some specific ADAS related machine learning tasks include detection of road boundaries, lane topologies, location of other cars, pedestrians, road signs and obstacles. These detection capabilities form the core of an ADAS technology stack. Other parts of the technology stack include decision making and control systems that take action in a certain road situation based on the input from the perception system. Although this picture works in controlled environments, making this technology effective in changing road situations, emergencies, changing weather etc. remains a significant challenge.

In recent times, deep learning has shown leading accuracy in a number of machine learning challenges. Specifically relevant to ADAS application is the dramatic increase in accuracy of image object classification 

[16, 22, 21, 12] and localization [20, 26, 10, 9, 19] in the last few years. A key advantage of DNN-based approaches is that they do not require hand tuned features for detecting every object but rather learn the representation from the data itself. Deep learning based perception systems promise to play a key role in navigation and safety software stack for ADAS.

R-CNN and its faster variants have become the state of the art in different object detection tasks. In this work, we leverage this method to establish a number of observations related to car detection on the challenging KITTI  [8] dataset. Our main results can be summarized as:

  • Bigger input images lead to higher accuracy. Input image resolution increases the accuracy of car detection using the faster R-CNN network.

  • Shallow models can deliver high accuracy. Convolutional features from shallow or earlier layers of DNNs lead to higher accuracy than features from the deeper layers. This holds true for deep models like VGG16. Surprisingly, even shallow models like AlexNet provide high accuracy on the detection task. Using shallow models that require less memory allow us to use very high input image resolutions. In terms of accuracy, shallow models with high resolution are competitive with deeper models with traditional resolution. This result is surprising given the trend of searching for deeper models for achieving high accuracy on object detection tasks.

The rest of the paper is organized as follows. In Section 2 we review related work, and we provide technical background information in Section 3. We describe our initial experimental setup in Section 4. We get to the crux of our results about large image resolution in Section 5 and shallow models in Section 6. We do additional exploration of R-CNN based configurations in Section 7 We summarize our findings in the context of the related work in Section 8, and we conclude in Section 9.

2 Related Work

2.1 Deep Networks for Object Detection

Deformable parts models (DPM) were the state of the art for image object detection [7] before the emergence of deep convolutional neural nets. The R-CNN method uses selective search for object region proposal [10]. The proposed regions in an image are warped to a fixed size and fed into a classification network called R-CNN. Fast R-CNN was introduced to reuse the shared convolution features for the region proposals  [9]. In Fast R-CNN [9], the inference speed is still dominated by the region proposal in the selective search method. Faster R-CNN [19] proposes object bounding boxes directly from the convolutional features. Inspired by the SPPNet [11] method, Faster R-CNN uses a region proposal network (RPN) to regress proposal boxes to ground truth boxes. The regions proposed by the RPN network is fed into the R-CNN network for classification. The network is trained end to end.111

The Faster R-CNN codebase also offers piecewise training of RPN and classifier branches of the network, but we found this cumbersome, and we use end-to-end training in all of our Faster R-CNN experiments.

Other than the RPN based method, there are several methods proposed for object bounding box prediction. For example, the OverFeat [20] method predicts a single box for localization whereas the Multibox [5, 23] method predicts multiple boxes in a class-agnostic way. The SPP method  [11] uses shared convolutional feature maps for fast object detection.

2.2 Detection on the KITTI dataset

2.2.1 Detection using 2D data

Deep neural networks are the backbone of most high-accuracy approaches to identifying objects such as cars in KITTI and similar datasets. Many such methods have been proposed; we focus this section on the high-accuracy and peer-reviewed results. A high-accuracy method for identifying objects in KITTI dataset is scale dependent pooling (SDP) combined with cascaded region classifiers (CRC) 

[25]. The crux of SDP+CRC lies in selecting a high-resolution CNN layer (e.g. conv3_3 in VGG16 [21]

) or a heavily downsampled CNN layer (e.g. conv5_3), depending on the resolution of each region proposal. By combining features from multiple convolution layers, they were able to achieve very high accuracy on KITTI’s object detection task. Our method introduced in this paper is even simpler in that we use only a single layer for feature extraction.

Another approach is Monocular 3D (Mono3D) which actually uses 2D images, but it aims to identify the pose of objects, with the goal of detecting objects as 3D bounding boxes. Like SDP+CRC, Mono3D is built around a version of R-CNN. There are also a number of anonymous and/or sparsely-explained submissions to the KITTI website’s leaderboard that are reportedly built on top of R-CNN.

2.2.2 Detection using 3D data

The KITTI dataset provides 3D information in the form of stereo images and LIDAR point clouds. Recent results such as 3DVP [24] and 3DOP [3] leverage both 2D and 3D data to achieve higher accuracy relative to comparable 2D baselines.

To build supervised 2D datasets such as ImageNet 

[4] and PASCAL [6], a widely-used approach is to have mechanical turk workers annotate user-generated images and videos from websites such as Flickr or YouTube. However, to our knowledge, there is no 3D equivalent of Flickr or YouTube that receives petabytes per week of user-uploaded 3D imagery. As a result, the overhead in building a 3D dataset currently requires not only data annotation, but also data collection. The cost of data collection includes hours of human labor, and it can also require expensive sensors. The KITTI dataset was released several years ago. However, the Velodyne HDL-64E LIDAR scanner used by the KITTI team still costs 80,000 USD222http://articles.sae.org/13899, which is more than twice the price of the average new car in the United States. With all of this in mind, we think widespread research on 3D object detection will be slow to emerge until (1) there is an internet hub that attracts large quantities of user-generated 3D imagery, and (2) the equivalent of today’s high-end LIDAR sensors become available for tens of US dollars. With this in mind, we focus our efforts on 2D imagery, where anyone with modest resources can collect a custom training set and apply our object detection approach.

3 Preliminaries

3.1 Relationship of conv and pooling strides to activation grid dimensions

The feature dimension of the output of a spatial convolution operation depends on the dimension of its input and the strides used (ignoring the boundary effects). In convolutional neural nets, strides are used to continually reduce the feature dimensions in the convolution layers. Additionally pooling layers are used to reduce dimension by averaging or taking maximum value within a neighborhood. Thus the spatial dimension of the output of the convolution decreases as we go to deeper layers. For example, in VGG16, for a standard input image size 224x224, the features calculated by the first convolution layer is 224x224 which reduce to 14x14 at the output of conv5_3. This dimension is further reduced by the pool5 layer to 7x7. Similarly for AlexNet, the convolution feature dimension reduces from 55x55 in conv1 to 6x6 in pool5.

In some CNN architectures such as AlexNet [16] and VGG16 [21], the first fully-connected layer expects a specific height and width for its input data (e.g. 6x6 for AlexNet). Increasing the height and width of the input image results in a higher-resolution input to the first FC layer. At the CNN architecture level, an easy way around this is to design a CNN architecture that has global average pooling prior to the first FC layer – this approach was popularized in the Network-in-Network (NiN) architecture [17], and it is now used in other architectures such as SqueezeNet [15] and ResNet architectures [12]. However, when using AlexNet or VGG19, the R-CNN authors developed a technique called ROI Pooling

that allows any size input image to be used in concert with AlexNet/VGG FC layers. ROI pooling is quite simple: no matter what the input image size is, ROI Pooling uses max-pooling to reshape the first FC layer’s input to the size that it expects.

3.2 Region Proposal Network in Faster R-CNN

We briefly review how the region proposal network (RPN) in Faster R-CNN generate proposals [19] that will be useful later. RPN starts with convolution layers, which computes a high dimensional, low resolution feature map for the input image. Next, a small network slides through each spatial position in the feature map and generates rectangular region proposals centered around the position. Instead of computing the proposal’s absolute coordinates, the RPN actually computes coordinates relative to a set of pre-selected reference boxes, or anchors. The transformation from an anchor to a proposal is illustrated in Fig. 1.

Figure 1: Transformation from an anchor box (left) to a proposal (right). 4 relative coordinates are regressed by the RPN network to adjust the center position and the shape of the bounding box.

Intuitively, we want the anchors to be spatially close to the ground truth bounding boxes. In an extreme case, if an anchor box is too far away from the ground truth bounding box, learning to transform the anchor to the ground truth will be hopeless. Since anchors are centered at each spatial position on the feature map, and each position on the feature map corresponds to a patch of pixels on the original image, the resolution of the feature map affects the distance from a ground truth bounding box to its nearest anchor. In VGG16, for example, each position in conv5_3 layer spatially corresponds to a patch on the original image, so in the worst case, the nearest anchor to the center of a ground truth bounding box is . As we will see later in this paper, reducing this distance, or relatively, increasing the “anchor density” will significantly increase the localization accuracy, thus improve the detection accuracy.

4 Experimental Setup

4.1 Networks and training configuration

We train faster R-CNN networks built on the VGG16 [21] and AlexNet [16], pretrained on the ImageNet-1k [4] classification dataset. VGG16 has sixteen convolution layers  [21] and AlexNet has only five convolutional layers [16]. Rather than using the convolution features of the last pooled layer (as is done in the original faster R-CNN paper), we use features from convolutional layers that are the bottom layers of the previous pooling layer. For VGG16 it is conv4_3. We resize the window size accordingly. For VGG16 the window size changes from 7x7 to 13x13. We also reduce the feature stride by a factor of two for avoiding one pooling stage. For VGG16 and AlexNet, fully connected layers are used as the R-CNN branch. The weights in these layers are initialized with random gaussian noise. As the standard procedure introduced in faster R-CNN, we randomly sample 128 positive and 128 negative proposals per batch to train the R-CNN layer. For all the experiments, we use initial learning rate of 0.0005, step size 50000 and momentum 0.9. A total of 70K iterations are run during R-CNN training starting from imagenet pre-trained weights for the convolution layers.

4.2 Dataset

We use the KITTI object detection dataset. KITTI object has three categories that is car, pedestrian, bicyclists. The dataset is annotated in three categories based on the occlusion and truncation of the objects. The hard category is heavily occluded and truncated whereas the easy category is relatively clearly visible. There are about 8000 images for both training and testing. The moderate regime is used to rank the competing methods in the benchmark. Our split of the KITTI train and validation sets (each containing half of the images) is the same as  [3]. The evaluation criteria is the same that is prescribed in the KITTI development kit. In KITTI’s evaluation criteria, proposal boxes having overlap with the ground truth or IoU greater than 70% are counted as true detection for cars.

The Faster R-CNN algorithm has been shown to deliver high accuracy on the PASCAL [6] dataset. Adapting that pipeline from PASCAL to KITTI poses a few natural challenges. First, the image sizes in the KITTI dataset is 1242x375 pixels whereas the image sizes in PASCAL dataset is 500 pixels in the longest dimension (many PASCAL images are 500x333 or 333x500). More importantly, the KITTI dataset contains heavily occluded and truncated objects. These objects come in multiple scales. The presence of objects at multiple scales make it difficult to attain high accuracy specially for small objects.

5 Input image resolution

We performed extensive design space search of Faster R-CNN configurations on the KITTI dataset. Our starting point is the VGG16 network that has achieved high accuracy in both image classification [21] and localization [19]. We performed an input image scaling experiment to find its impact on the accuracy. In these experiments, the shorter side of the KITTI images were fixed at 1295 pixels. In the Faster R-CNN codebase, the default off-the-shelf configuration resizes all images to 1000 pixels in the long dimension.

KITTI images have a native resolution of 1295x375, so the default Faster R-CNN behavior is to resize KITTI images to 1000x302. But, is this resizing scheme ideal for obtaining high accuracy? To find out, we doubled the input image height and width to 2000x604. This has the effect of doubling the height and width of the activations (outputs) from all convolutional layers. For example, the conv5_3 activations – which serves as input to both the region proposal network (RPN) and the classification network – double in height and width. With the image upsampled to 2000x604, we see in Table 1 that the KITTI car-detection accuracy increases for easy, medium, and hard by 7.1, 15.5, and 12.6 percentage points, respectively. In a world where half of a percentage point is considered significant, we can say with certainty that the input resolution has a major impact on accuracy.

Can further upsampling of the image lead to further improvements in accuracy? We attempted to perform experiments with upsampling beyond 2000x604, but the volume of activation planes exceeded the 12GB of available memory on an NVIDIA Titan X GPU. In the next section, we consider shallower networks with fewer layers of activation planes, which enables us to move to even higher input resolutions.

AP
Input resolution Easy Medium High
1000x302 80.3 63.0 52.3
2000x604 87.4 78.5 64.9
Table 1: KITTI car detection accuracy using different input image sizes to VGG16. In these experiments, we use conv5_3 features from VGG16.

6 Shallow convolutional models

So far, we have upsampled the input image until we ran out of on-chip GPU memory when training R-CNN models with a VGG16-based feature representation. Based on what we have seen so far, it seems that further upsampling may lead to further gains in accuracy. We need to find a configuration that requires less memory for a given image size, and then we will exploit this extra memory to further upsample the image. One idea would be to decrease the batch size to save memory, but we are already using a batch size of 1, and it’s not clear how to reduce the batch size below 1. Could we reduce the memory footprint by reducing the number of layers in the CNN? Much of the recent literature shows that fewer layers in a CNN leads to lower accuracy (all else held equal). But, our goal is to configure a CNN with fewer layers (and moderately lower accuracy) and then increase the input image resolution (leading to much higher accuracy).

To evaluate this idea, we configure R-CNN to use 2000x604 images (2x height and 2x width compared to our original starting point), using conv4_3 instead of conv5_3 features. In VGG16’s scheme for naming layers, conv4_3 is the 10th layer, and conv5_3 is the 13th layer in the CNN. We expected that the accuracy of R-CNN with conv4_3 would be slightly lower than R-CNN with conv5_3, but as we show in Table 2 that the accuracy is higher with conv4_3 by 5.5, 9.4, and 12.4 percentage-points for easy, medium, and hard detections. Cumulatively, the improvement in accuracy from conv5_3 with 1000x302 images to conv4_3 with 2000x604 images is a whopping 12.6, 24.9, and 25.0 percentage-points for easy, medium, and hard, respectively.

How does further reducing the CNN’s depth affect accuracy? We initially considered using the earlier layers of VGG16 as input to the Region Proposal Network. But, earlier layers in VGG16 have been downsampled less, so their activations have a larger height and width. We found that the off-the-shelf implementation of RPN comes to dominate the end-to-end computation time with very large height and width input grids. Besides depth, one of the differences between VGG16 and AlexNet is that AlexNet downsamples more aggressively in the early layers – for example AlexNet has stride=4 in the conv1 layer (4x downsampling), while the conv1 layer of VGG16 has stride=1 (no downsampling). So, to evaluate this question of how using shallower (10 conv layers) network impacts accuracy, we use AlexNet instead of VGG16. We use conv5 (5th layer) activations as input to the R-CNN region-proposal and classification branches, and we report the results in Table 2. With resolution of 2000x604 for both AlexNet-conv5 and VGG16-conv4_3, the VGG16-based configuration delivers significantly higher accuracy on easy, medium, and hard detections in Table 2. We have additional memory available when running AlexNet with 2000x604 input images, so we now try upsampling the AlexNet input images to 5000x1510. In this configuration, on the easy detections, AlexNet with 5000x1510 input is within 0.5 of a percentage-point of our best VGG16-based result so far. On medium and hard categories, VGG16 conv4_3 with an input resolution of 2000x604 delivers higher accuracy than AlexNet with 2000x604 or 5000x1510 input images.

We also conduct a sweep of input image sizes applied to an AlexNet-based R-CNN model that uses conv5 features. We show the results of this sweep in Figure 2. We observe that KITTI car detection accuracy steadily climbs from a baseline resolution of 1000x302 to a plateau at resolution 5000x1510. Beyond this resolution, we have not observed further accuracy improvements at larger sizes such as 6000x1811 or 7000x2133.

Figure 2: Evaluation of how input resolution affects accuracy on KITTI car detection. In this configuration, we use AlexNet conv5 as input features to the R-CNN branches.
AP
CNN Architecture Conv layer name (depth) Input image resolution Easy Medium Hard
VGG16 conv5_3 (13) 2000x604 87.4 78.5 64.9
VGG16 conv4_3 (10) 2000x604 92.9 87.9 77.3
AlexNet conv5 (5) 2000x604 86.7 71.6 56.1
AlexNet conv5 (5) 5000x1510 92.4 82.5 68.2
Table 2: Impact of CNN depth on accuracy. Conventional wisdom would suggest that deeper representations would produce higher accuracy, but we find otherwise. AP numbers are for car detection on the KITTI dataset.

7 Further improvements

7.1 Context windows

Chen et al. [3] proposed context windows as a way to include information from adjacent pixels of a proposed bounding box. A context window is a bounding box that is scaled up from the original bounding box proposal of the RPN network. In the experiments with context window, in addition to the original R-CNN branch, an extra R-CNN branch is added that trains on the features extracted from context window. The original R-CNN features and the context R-CNN features are concatenated before classification. We add a context branch in addition to the usual R-CNN branch with a spatial bounding box scaling of 1.5. When applying the context window to an AlexNet-based R-CNN configuration, we find that the accuracy of all the categories improve as shown in Table 4. The improvement is significant in small image sizes and provides diminishing returns as we scale up the input image size.

7.2 Optimal Anchor-box shape

(a) Distribution of bounding box shapes of cars in the KITTI dataset
(b) 9 Anchor box shapes (white crosses) selected by K-Means
Figure 3: Bounding box shape distribution of the car object in the KITTI dataset is plotted on the left. 9 anchor shapes computed by K-Means are plotted in the right figure as white crosses.

In [19], default anchor shapes are arbitrarily chosen by reshaping a square box by 3 scales and 3 aspect ratios. But is there a better way to choose anchor shapes for input images and target objects? Intuitively, we want the anchors to have similar shapes with the ground truth bounding boxes. The shape of a bounding box can be characterized by its width and its height . The width and height distribution of the car object in the KITTI training data set is plotted in Fig. 3(a). The problem of choosing the ”most similar” anchor shapes can be formulated as the following: given a set of ground truth bounding box shape observations , find anchors such that the sum of the distance (in the shape space of ) between each ground truth box to its nearest anchor is minimized. This problem can be effectively solved by K-means. The optimal anchor shapes are plotted in Fig. 3(b). These anchors are optimized specifically for the car category, but the idea of optimizing anchors by considering ground truth bounding box statistics can be generalized to multi-category object detection as well.

We tested the AlexNet-based Faster R-CNN’s detection accuracy with different anchor box selection schemes, and the result is shown in Table 3. We fix the number of anchors to be 9. From the bounding box shapes of cars in the training data set, we used K-Means to select 9 optimal anchor boxes. As comparison, we used a set of default anchors with 3 scales and 3 aspect ratios as in [19]. When image width is , using the K-Means selected anchors improves AP significantly comparing with the default anchors. When the image width is , we could see that AP with default shapes are slightly better for easy and moderate category, but using K-Means selected anchors still improves AP for the hard category by 4 percentage points. As we further scale the image width to , the performance gain saturates, we still observe some improvement by using K-Means selected anchors. We have not yet combined context windows with our anchor box improvements, but it is possible that this combination will yield a further improvement in accuracy.

DNN architecture Input resolution Anchor shape selection scheme AP
Easy Moderate Hard
AlexNet 1242x375 Default shape 70.37 54.44 46.33
AlexNet 1242x375 K-Means 76.12 59.29 47.38
AlexNet 2500x755 Default shape 84.12 72.14 58.29
AlexNet 2500x755 K-Means 83.27 71.43 62.42
AlexNet 5000x1510 Default shape 91.33 84.52 69.90
AlexNet 5000x1510 K-Means 91.44 85.98 70.04
Table 3: Impact of anchor box shape on accuracy. These results use AlexNet conv5 features.
AP
Source CNN Architecture Feature layer (Depth) Input resolution Context window Easy Medium Hard
SDP+CRC [25] VGG16 conv3_3, 4_3, 5_3 (7,10,13) multiple no 90.3 83.5 71.1
Mono3D [2] VGG16 conv5_3 (13) not reported yes 92.3 88.7 79.0
ours VGG16 conv5_3 (13) 1000X302 no 80.25 62.96 52.3
ours VGG16 conv5_3 (13) 2000X604 no 87.35 78.49 64.93
ours VGG16 conv4_3 (10) 2000X604 no 92.9 87.9 77.3
ours VGG16 conv4_3 (10) 5000x1510 no out of memory
ours AlexNet conv5 (5) 1000X302 no 67.5 49.44 38.9
ours AlexNet conv5 (5) 2000X604 no 86.7 71.6 58.1
ours AlexNet conv5 (5) 5000x1510 no 92.4 82.5 68.2
ours AlexNet conv5 (5) 1000X302 yes 71.58 51.13 40.9
ours AlexNet conv5 (5) 2000X604 yes 86.98 74.32 60.83
ours AlexNet conv5 (5) 5000x1510 yes 94.7 84.8 68.3
Table 4: Summary of results on KITTI [8] car detection. All of our results are based on Faster R-CNN. To our knowledge, all of the related work discussed in this table also uses a version of R-CNN.

8 Discussion

We show the precision-recall curve for our KITTI car detection in Fig. 4(a). We used AlexNet-conv5 with a context window and input image resolution of 5000x1510, as in the final row of Table 4. We observe that the precision of the easy category is very high even at very high recall. However, the precision in the hard category suffers at high recall. Improving the precision of results on the hard category with our method will be a target of future work. We show a few examples of success and failure modes in the hard category in Fig 5. In Fig 5(a), the model successfully predict a highly occluded car while in Fig 5(b) the predicted bounding box encompass two cars that are adjacent to each other. In Fig 5(c), the predicted bounding box enclose a visible car but completely misses the car that is truncated. The precision-recall curve using the conv4_3 features of VGG16 is shown in Fig.4(b). The precision at high recall for the hard category improves significantly. The inference time using AlexNet conv5 layer with input image size of 5000x1510 and VGG16 conv4_3 layer with input image size of 2000x604 is 0.34s and 0.6s respectively. Inference times for other published high accuracy methods on the KITTI dataset are 3s for 3DOP [3] and 0.4s for SDP+CRC [25].

(a) Precision-recall curve for the AlexNet network with context window. The input image size is 5000. The conv5 features are used in this experiment.
(b) Precision-recall curve for the VGG16 network without context window. The input image size is 2000. The conv4_3 features are used in this experiment.
Figure 4: Precision vs. recall curve for KITTI’s car detection system.
(a) Success.
(b) Predicting one box for two adjacent cars.
(c) Failure
Figure 5: Examples of success(a) and failure modes(b,c) in the hard category using AlexNet with a context window and input resolution of 5000x1510 pixels. The green boxes show the ground truth label of heavily truncated car. In (a), the model successfully predict a highly occluded car. In (b), the predicted bounding box encompass two cars that are adjacent to each other. In (c), the predicted bounding box enclose a visible car but completely misses the car that is truncated.

9 Conclusions

In summary, we have shown that shallow networks perform well in achieving high accuracy on detecting cars in the road. We have shown that input image resolution has a large impact on the accuracy of car detection using the faster R-CNN network. For very deep models, shallow layers can (surprisingly) provide higher accuracy than the later convolutional layers deep in the network. Shallow models like AlexNet can achieve high accuracy when the input image is upsampled. In addition, we have used an anchor box selection method and context window to further enhance car detection accuracy. We believe that our findings will inspire the research community to evaluate shallow models for achieving high accuracy on object detection tasks.

Acknowledgments

Khalid Ashraf was supported by the National Science Foundation under Award number 125127. We thank Kostadin Ilov for many help with computational hardware. Thanks to Ross Girshick for comments on some initial results. Thanks to Fan Yang and Kaustav Kundu for clarification on their results.

References