Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1611 1579 (N27W52 L=314)   931    114        4 (1573/1533-1653/1632)  81.9 W48-W56 N23-N30
      2 1184 1529 (N29W17 L=278)   766    219        0 (1111/1437-1244/1589) -67.4 W13-W19 N21-N33
      3 1660 1401 (N15W48 L=310)   356     82        1 (1633/1356-1687/1437) -68.2 W47-W50 N12-N17
      4 1820 1054 (S07W56 L=317)   374     49        0 (1808/1022-1834/1085) -76.0 W54-W57 S09-S05
      5 1844  912 (S16W59 L=320)   473     75        0 (1822/ 878-1875/ 940)  31.0 W56-W62 S17-S14
      6 1700  627 (S33W47 L=308)   960    138        4 (1642/ 590-1745/ 671)  -9.5 W41-W53 S35-S29
      7 1398  563 (S33W21 L=282)   646    139        1 (1373/ 495-1418/ 626)  83.7 W18-W22 S38-S29
      8 1528  555 (S36W32 L=293)   408    108        0 (1486/ 515-1573/ 597)  35.5 W28-W35 S38-S33
      9 1079  439 (S38E04 L=257)   368     89        0 (1051/ 402-1100/ 482)  45.0 E06-E01 S40-S35
     10  826  421 (S36E23 L=238)   460     90        1 ( 789/ 385- 872/ 457) -63.4 E25-E20 S39-S33
     11  886  238 (S52E28 L=233)   399     59        0 ( 849/ 221- 919/ 254)  -8.1 E31-E24 S53-S50
     12  753  168 (S56E53 L=208)   624     69        0 ( 737/ 133- 771/ 206)  16.9 E66-E43 S59-S52
  

    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