Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1315 1721 (N36W44 L=101)   663     97        2 (1273/1694-1359/1742) -26.6 W40-W47 N33-N37
      2 1620 1495 (N14W52 L=109)   402     63        0 (1588/1467-1650/1523) -63.4 W50-W54 N11-N15
      3  993 1432 (N28W10 L= 67)   288     78        0 ( 959/1402-1034/1449)   9.5 W06-W12 N27-N29
      4 1094 1265 (N16W10 L= 68)   298    109        0 (1053/1234-1131/1309)  18.4 W07-W13 N14-N17
      5  154 1252 (N39E60 L=357)   785    211        0 ( 117/1234- 192/1268) -18.4 E64-E55 N37-N40
      6  836 1207 (N20E06 L= 51)   463    118        1 ( 767/1184- 892/1224) -35.5 E09-E02 N17-N22
      7  274  966 (N19E50 L=  7)   538    125        0 ( 228/ 943- 321/ 984)  -7.1 E54-E46 N17-N20
      8 1808  742 (S35W47 L=104)  2127    235        1 (1735/ 641-1858/ 844)  39.8 W39-W52 S39-S28
      9 1313  602 (S27W06 L= 63)  1324    331        1 (1282/ 447-1362/ 741)  58.6 W01-W12 S38-S18
     10  974  517 (S23E17 L= 40)   996    285        0 ( 897/ 414-1020/ 643)  36.3 E25-E10 S27-S16
     11  690  467 (S19E36 L= 21)   987    149        3 ( 649/ 409- 736/ 529)  63.4 E39-E33 S22-S14
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y