Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1354 1808 (N55W08 L=194)   772    105        0 (1305/1788-1400/1827)  26.6 W03-W13 N51-N58
      2 1360 1724 (N47W10 L=195)   493     63        1 (1339/1694-1383/1756)  74.1 W07-W11 N43-N50
      3  501 1463 (N11E42 L=144)   402    110        0 ( 483/1420- 519/1521) -63.4 E43-E39 N08-N14
      4 1775 1276 (N25W46 L=232)   873    161        1 (1754/1209-1789/1359)  84.3 W44-W47 N20-N30
      5 1121 1270 (N09W01 L=187)   257     61        1 (1080/1243-1156/1292)  14.0 E01-W03 N07-N10
      6 1221  875 (S12W14 L=200)  1238    410        3 (1053/ 709-1322/1040)  68.6 W07-W17 S24-N00
      7  310  878 (S26E47 L=139)   374     55        1 ( 282/ 855- 336/ 903)  36.9 E48-E44 S27-S24
      8  628  842 (S24E22 L=164)   728    220        0 ( 515/ 787- 721/ 903)  39.8 E29-E15 S29-S19
      9 1641  816 (S06W43 L=228)   248    112        0 (1595/ 799-1683/ 827)   8.1 W39-W46 S07-S04
     10  910  752 (S25E02 L=184)   458     70        3 ( 865/ 711- 941/ 797)  71.6 E03-W00 S27-S21
     11 1536  495 (S26W47 L=233)  5081    941        6 (1153/ 413-1838/ 617)  24.8 W23-W77 S40-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