Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  544 1653 (N26E45 L=232)   985    111        2 ( 498/1597- 598/1700)  69.4 E47-E42 N21-N29
      2 1170 1642 (N35E01 L=276)   801    122        2 (1105/1608-1248/1678)   6.3 E05-W04 N32-N36
      3 1136 1404 (N18W01 L=278)   686    135        3 (1105/1336-1170/1460) -82.9 E01-W02 N13-N21
      4  575 1375 (N08E33 L=244)   466    100        1 ( 536/1333- 604/1412) -33.7 E36-E30 N05-N10
      5 1370 1280 (N14W17 L=294)   480    142        1 (1326/1218-1417/1336)  69.4 W14-W19 N10-N18
      6 1446 1236 (N13W22 L=299)   429     90        0 (1419/1198-1489/1261)  51.3 W20-W24 N10-N15
      7 1727 1041 (N07W45 L=322)   369     43        0 (1712/1008-1743/1069) -56.3 W43-W46 N05-N09
      8  537  874 (S22E29 L=248)  4114   1205        3 ( 173/ 596- 929/ 988)  20.1 E64-E04 S42-S10
      9 1804  940 (N03W54 L=331)   417     52        1 (1784/ 919-1831/ 963)  56.3 W52-W56 N01-N04
     10 1356  662 (S22W27 L=305)   813    202        1 (1258/ 626-1458/ 698)  33.7 W21-W33 S25-S17
  

    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