[out3Plot]: add functionality for assets

[awirb]

No description provided.

[qltysh]

❌ 2 blocking issues (2 total)

Tool	Category	Rule	Count
black	Style	Incorrect formatting, autoformat by running qlty fmt.	2	❌

@qltysh one-click actions:

• Auto-fix formatting (qlty fmt && git push)

[qltysh]

Coverage Impact

⬇️ Merging this pull request will decrease total coverage on master by 0.10%.

Modified Components (1)

Rating	Component	% Diff
	com1DFA	95.2%

Modified Files with Diff Coverage (4)

Rating	File	% Diff	Uncovered Line #s
	avaframe/in1Data/getInput.py	3.3%	1209-1250
	avaframe/com1DFA/com1DFA.py	100.0%
	avaframe/com1DFA/particleTools.py	94.1%	1054
	avaframe/com1DFA/deriveParameterSet.py	100.0%
	Total	41.2%

🤖 Increase coverage with AI coding...

In the `addInfra` branch, add test coverage for this new code:

- `avaframe/com1DFA/particleTools.py` -- Line 1054
- `avaframe/in1Data/getInput.py` -- Line 1209-1250

🚦 See full report on Qlty Cloud »

🛟 Help

• Diff Coverage: Coverage for added or modified lines of code (excludes deleted files). Learn more.

• Total Coverage: Coverage for the whole repository, calculated as the sum of all File Coverage. Learn more.

• File Coverage: Covered Lines divided by Covered Lines plus Missed Lines. (Excludes non-executable lines including blank lines and comments.)

  • Indirect Changes: Changes to File Coverage for files that were not modified in this PR. Learn more.

[fso42]

Wrong PR, sorry

[fso42]

  def createAssetsRasterFromParticleLocations(particlesTimeArrays, dem, uniqueAssets,   
assetsValues):                                                                          
      """create a raster indicating particle trajectories colorcoded with assets        
classes,                                                                                
      highest overrides lower classes"""                                                
                                                                                        
      nTime, nPart = particlesTimeArrays["ID"].shape                                    
      particlesTimeArrays["assetsValue"] = np.full((nTime, nPart), np.nan)              
      particleAssets = np.full((dem["header"]["nrows"], dem["header"]["ncols"]),        
np.nan)                                                                                 
                                                                                        
      # Process classes from low to high so higher classes naturally override lower     
ones                                                                                    
      for l in uniqueAssets:                                                            
          assetCells = np.asarray(assetsValues["value_%d" % l])                         
          inAsset = np.isin(particlesTimeArrays["inCellDEM"], assetCells)               
                                                                                        
          for pId in range(nPart):                                                      
              hitTimes = np.where(inAsset[:, pId])[0]                                   
              if len(hitTimes) == 0:                                                    
                  continue                                                              
                                                                                        
              # Mark particle trajectory up to the last time it hit this asset class    
              mMax = hitTimes[-1]                                                       
              particlesTimeArrays["assetsValue"][:mMax, pId] = l                        
                                                                                        
              indX = particlesTimeArrays["indXDEM"][: mMax + 1, pId].astype(int)        
              indY = particlesTimeArrays["indYDEM"][: mMax + 1, pId].astype(int)        
              # Only overwrite cells not already set to a higher class                  
              particleAssets[indY, indX] = np.where(                                    
                  particleAssets[indY, indX] >= l, particleAssets[indY, indX], l        
              )                                                                         
                                                                                        
      return particleAssets, particlesTimeArrays