Numpy uses inverse indexing order (k, j, i).
In the example above, va[138, 277, 163] should give you the results you expect.
See further examples in script repository, such as https://www.slicer.org/wiki/Documentation/Nightly/ScriptRepository#Get_markup_fiducial_RAS_coordinates_from_volume_voxel_coordinates