Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  603 1756 (N31E53 L=355)  1701    397        3 ( 555/1666- 652/1826) -27.8 E62-E43 N25-N35
      2 1041 1547 (N28E13 L= 35)  1271    370        2 ( 900/1429-1191/1666) -13.0 E26-W00 N23-N32
      3  516 1605 (N18E50 L=358)  1411    384        4 ( 423/1510- 559/1664) -19.4 E61-E43 N13-N21
      4  365 1548 (N12E61 L=347)  1180    192        3 ( 332/1478- 408/1621) -71.6 E64-E57 N07-N16
      5 1474 1498 (N37W17 L= 66)   520     91        1 (1438/1461-1517/1534) -15.9 W14-W21 N35-N38
      6  349 1433 (N05E56 L=352)   752    170        3 ( 293/1394- 405/1471)  36.9 E61-E51 N01-N07
      7  655 1435 (N11E33 L= 15)   257     75        1 ( 638/1402- 669/1469) -36.9 E34-E30 N09-N12
      8  495 1391 (N05E42 L=  6)   266     75        0 ( 484/1346- 508/1429) -68.2 E44-E40 N02-N07
      9  567 1244 (S02E33 L= 16)   553    123        0 ( 519/1224- 624/1271)   0.0 E36-E28 S03-N00
     10  709  570 (S39E09 L= 40)   392     88        0 ( 664/ 543- 759/ 595)   7.1 E12-E04 S39-S37
  

    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