Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  472 1729 (N25E63 L=148)   748    158        2 ( 443/1693- 487/1772) -63.4 E66-E61 N22-N28
      2  295 1580 (N12E73 L=138)  1852    234        1 ( 260/1542- 326/1620)  56.3 E79-E70 N09-N15
      3  241 1397 (N00E68 L=144)   897    152        1 ( 206/1367- 272/1426)  56.3 E72-E63 S01-N02
      4  946  910 (S12E02 L=210)   684    216        1 ( 872/ 836-1063/ 967)  54.0 E04-W04 S18-S07
      5 1763  883 (N08W52 L=263)   474     90        1 (1739/ 847-1791/ 922) -38.7 W49-W55 N06-N10
      6 1294  795 (S10W21 L=233)   394    125        0 (1260/ 731-1320/ 859) -83.7 W19-W22 S14-S05
      7  482  672 (S38E28 L=183)  1000    273        1 ( 379/ 651- 595/ 692)  21.3 E37-E19 S41-S34
      8  851  583 (S34W01 L=213)   890    221        0 ( 730/ 535- 948/ 660)  -6.7 E08-W09 S35-S32
      9 1319  513 (S25W33 L=245)  1862    483        1 (1089/ 444-1486/ 572)   9.2 W17-W48 S28-S21
     10 1695  453 (S17W67 L=278)  1188    253        1 (1644/ 428-1739/ 483)   3.6 W59-W75 S18-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