This paper presents a technique, combining the integral equations (IE) and the Generalized Sheet Transition Conditions (GSTCs) with bianisotropic susceptibility tensors, to compute electromagnetic wave scattering by cylindrical metasurfaces -forming two-dimensional porous cavities -of arbitrary cross sections. Moreover, it applies this technique to two problems -cloaking with circular and rhombic shapes and illusion optics with an elliptic shape -that both validate it, from comparison with specifications used in an exact synthesis of the metasurfaces, and reveal interesting capabilities of such metasurface structures. Particularly, active cylindrical metasurfaces can perfectly cloak and hence eliminate the extinction cross section of various cylindrical shapes, and simple purely passive versions of them, practically more accessible, still perform quite good cloaking and provide remakable extinction cross section reduction. onto the x and y directions of the global system, while leaving the shared coordinate z unchanged. We have then, for instance,