Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  846 1602 (N33E24 L= 14)  1280    332        2 ( 728/1512- 972/1685) -19.2 E35-E12 N28-N37
      2  445 1540 (N22E51 L=347)  1088    287        2 ( 382/1424- 512/1657) -86.6 E53-E49 N13-N30
      3 1249 1577 (N39W06 L= 44)   939    202        2 (1151/1529-1343/1638) -14.9 E02-W13 N36-N42
      4 1034 1508 (N30E08 L= 30)   269     47        1 (1015/1479-1058/1531) -63.4 E09-E06 N27-N31
      5 1430 1502 (N37W21 L= 59)   302     65        2 (1400/1485-1461/1516)  21.8 W18-W23 N35-N38
      6  922 1417 (N22E13 L= 25)   301     70        0 ( 886/1385- 953/1439) -21.8 E15-E10 N20-N23
      7  389 1142 (S03E44 L=354)   357     80        1 ( 353/1116- 426/1160)  26.6 E46-E40 S05-S01
      8 1167  472 (S32W20 L= 59)   512    148        1 (1102/ 451-1241/ 504)   0.0 W15-W26 S33-S31
  

    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