Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1260 1611 (N41W24 L=112)   336     44        0 (1233/1591-1291/1633) -45.0 W21-W25 N39-N43
      2 1327 1580 (N39W28 L=117)   386     69        2 (1306/1551-1354/1609)  53.1 W25-W30 N36-N40
      3 1763 1387 (N21W59 L=148)   671    134        2 (1733/1350-1792/1425)  33.7 W55-W63 N19-N23
      4 1673 1383 (N22W50 L=139)   525     85        1 (1651/1344-1694/1418)  68.2 W48-W52 N19-N24
      5  890 1358 (N26E07 L= 82)   208     39        0 ( 871/1333- 911/1383)  45.0 E08-E05 N24-N27
      6 1761  681 (S23W53 L=142)  3171    474       11 (1693/ 519-1822/ 819)  78.7 W47-W58 S34-S14
      7  931  660 (S18E08 L= 81)   220     95        0 ( 892/ 645- 976/ 685) -31.0 E10-E05 S18-S15
      8  463  594 (S21E42 L= 47)   362    124        0 ( 429/ 565- 490/ 636)  33.7 E45-E39 S22-S17
      9 1298  577 (S26W16 L=105)   638    110        2 (1249/ 539-1354/ 606)  23.2 W12-W20 S28-S23
     10 1052  375 (S39E02 L= 87)  1062    195        2 ( 952/ 331-1116/ 431)  26.6 E10-W03 S42-S35
  

    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