Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1212 1757 (N46E15 L=220)  1664    228        2 (1141/1684-1289/1803)  81.9 E18-E08 N38-N52
      2 1543 1436 (N34W21 L=255)  1034    157        2 (1511/1352-1566/1506) -85.2 W17-W22 N28-N40
      3 1733 1352 (N36W39 L=274)  1761    287        1 (1657/1214-1788/1454) -35.5 W31-W46 N29-N41
      4  747  786 (S27E09 L=225)   345    137        0 ( 696/ 749- 802/ 821)  54.5 E12-E06 S30-S23
      5 1196  707 (S18W19 L=254)   654    263        2 (1106/ 673-1287/ 777)  24.4 W12-W24 S20-S14
      6  520  646 (S42E23 L=211)   257     40        1 ( 501/ 617- 539/ 671) -45.0 E25-E21 S43-S40
      7 1706  646 (S05W56 L=291)   404     29        2 (1684/ 626-1730/ 667)  11.3 W53-W58 S06-S03
      8  971  492 (S38W15 L=250)   565     99        1 ( 934/ 452-1027/ 527)  74.1 W12-W17 S40-S33
      9 1118  357 (S41W33 L=267)   309     54        0 (1091/ 339-1148/ 379)  76.0 W30-W34 S42-S38
     10 1312  238 (S40W59 L=293)   448     60        0 (1284/ 219-1345/ 265)   0.0 W53-W63 S41-S39
  

    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