Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1030 1550 (N32W08 L=158)   507    138        0 ( 974/1491-1082/1600) -56.3 W04-W10 N27-N36
      2 1872 1316 (N07W69 L=219)   884     59        2 (1849/1285-1905/1354) -33.7 W66-W74 N04-N09
      3  433 1271 (N22E37 L=113)   746    142        1 ( 373/1238- 499/1296)   5.7 E41-E31 N20-N23
      4  383 1213 (N19E41 L=109)   410     85        1 ( 340/1188- 428/1232)   7.1 E44-E36 N17-N19
      5 1599  897 (S15W35 L=185)  1370    343        2 (1558/ 743-1631/1050)  87.0 W33-W37 S24-S05
      6 1886  915 (S17W64 L=214)  3782    263        9 (1817/ 827-1936/1007) -40.2 W55-W74 S22-S11
      7 1485  770 (S21W27 L=176)   267    110        0 (1447/ 733-1512/ 821)  45.0 W23-W28 S23-S18
      8 1652  551 (S38W44 L=194)  1748    285        1 (1536/ 472-1762/ 638)  18.4 W32-W54 S42-S33
      9  598  571 (S22E35 L=115)   266     61        1 ( 584/ 538- 615/ 608)  38.7 E37-E32 S24-S20
     10  315  570 (S19E60 L= 90)   506     48        3 ( 295/ 546- 339/ 597) -90.0 E62-E57 S20-S16
     11  284  520 (S21E68 L= 82)   548    111        0 ( 251/ 480- 324/ 569) -74.1 E70-E64 S24-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